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CHAPTER  1 


introduction 

1.1  Problem  Definition  and  Issues 

A broadcast  channel  is  a communications  channel  where  a signal  generated 
by  one  transmitter  can  be  received  by  many  receivers.  Some  examples  of  these 
channels  are:  a satellite  channel  where  a satellite  acts  as  a transponder; 
a ground  radio  network  where  a number  of  terminals  have  access  to  common  freq- 
uencies; and  coaxial  cables  where  a number  of  terminals  use  a single  cable  for 
communicating  by  time  multiplexing.  When  a number  of  spatially  Isolated  Inde- 
pendent sources  have  access  to  a common  broadcast  channel,  a problem  usually 
arises  In  allocating  the  channel  capacity  to  the  various  users.  This  Is  so 
because  the  only  means  these  Independent  users  have  of  communicating  their 
requests  for  channel  capacity  Is  through  the  channel  Itself. 

There  are  three  main  Issues  that  are  Involved  In  multl-accesslng  a 
comann  channel:  a)  the  percentage  of  the  capacity  used  In  accessing,  b)  the 
time  It  takes  to  access  the  channel,  and  c)  the  stability  of  the  multl-accesslng 
system,  l.e.,  how  likely  Is  It  for  the  system  to  be  In  a state  where  many 
sources  are  attempting  to  access  the  channel  and  very  few  are  actually 
succeeding.  These  three  properties  (throughput,  delay,  and  stability)  will 
be  defined  more  precisely  later  on;  however.  It  should  be  clear  at  this  point 
that  a desirable  multl-accesslng  system  Is  one  that  Is  stable,  has  high 
throughput  and  low  accessing  delay.  As  we  will  see  In  Section  1.3,  when  the 
nuiber  of  sources  Is  large  and  the  message  lengths  short,  present  multl- 
accesslng  schemes  suffer  In  at  least  one  of  these  attributes.  In  this  thesis. 
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therefore,  we  will  be  mainly  concerned  with  the  development  of  stable  multi- 
accessing  techniques  that  have  high  throughput  and  low  delay. 

More  specifically,  we  will  be  concerned  with  the  time  domain  multl- 
accesslng  of  a common  broadcast  channel  by  having  the  sources  transmit 
their  data  In  packet  form.  A packet  Is  a block  of  fixed  length  digital  data 
that  contains  the  Information  to  be  transmitted,  along  with  the  source  and 
destination  addresses  and  any  other  overhead  Information  that  might  be 
necessary,  such  as  error  correcting  and  error  detecting  bits.  We  will 
assume  that  If  more  than  one  transmitter  transmits  simultaneously,  then,  they 
will  Interfere  with  each  other,  and  all  the  packets  will  be  received 
Incorrectly.  If,  on  the  other  hand,  no  packet  collisions  occur,  then  we 
assume  that  error  free  transmission  results. 

It  Is  Important  for  each  transmitter  (also  to  be  referred  to  as  a 
source)  to  be  able  to  determine  whether  there  are  zero,  one,  or  more  than 
one  packets  In  the  channel  at  any  one  time.  More  than  one  packet  corresponds 
to  a collision,  l.e.,  a channel  contention.  This  channel  state  Information 
My  be  obtained  directly  by  listening  to  the  channel  or  by  some  other  means 
such  as  a central  observer  along  with  an  auxiliary  feedback  channel.  How 
the  transmitters  obtain  this  Infonsatlon  Is  not  Important  for  our  work.  Once 
it  Is  determined  that  a collision  occurred,  then  the  sources  must  take  action 
to  resolve  this  conflict.  The  resolution  of  this  conflict  Is  the  heo*’t  of  the 
multi-aecesslng  problem,  and  It  Is  here  that  we  will  focus  our  attention  by 
developing  and  analyzing  a new  class  of  conflct  resolving  algorltlns.  What 
makes  this  problem  Interesting  Is  that,  when  an  Initial  collision  occurs, 
each  of  the  contending  sources  knows  that  Its  packet  collided;  however.  It 
does  not  know  the  identities  or  the  number  of  the  other  contending  sources. 
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The  nultlple  access  system  to  which  our  work  Is  directly  applicable 
Is  the  so  called  packet  switching  broadcast  network  where  the  packets  might 
contain  data  from  such  sources  as  computers,  teletype  terminals,  or  vocoders. 
The  results  of  this  thesis,  however,  are  equally  applicable  to  the  various 
dynamic  reservation  multi-accessing  systems  In  which  there  are  two  channels, 
a data  channel  and  a reservation  channel.  In  this  system  the  reservation 
channel  Is  used  to  make  reservation  requests  for  the  data  channel.  Here  the 
channel  contention  difficulties  arise  In  the  reservation  channel  and  again 
one  must  resolve  the  conflicts  between  simultaneous  requests  on  this  channel. 

Our  results,  however,  need  not  be  restricted  to  conminlcatlon  systems. 
They  may  be  extended  to  more  general  systems  In  which  a central  facility  Is 
accessible  by  a number  of  Independent  users.  If  the  number  of  users  that 
the  facility  can  service  slmiltaneously  Is  less  than  the  maximum  that  can 
place  demands  upon  It,  then  contentions  will  arise  that  might  be  solvable 
with  the  techniques  developed  here. 

From  the  preceding  discussion,  we  see  that  a multi-access  system  may 
be  decomposed  Into  three  major  components:  the  sources  and  messages,  the 
channel  and  the  multl-accesslng  protocol.  These  are  discussed  In  more  detail 
In  the  following  sections.  In  Section  1.2  we  present  the  channel  and  source 
models  that  will  be  used  la  this  thesis  and  In  Section  1.3  we  consider  the 
multl-scceaslng  algorithms.  In  Section  1.4  we  present  an  outline  of  the 
thesis,  the  analysis,  and  the  main  results.  Finally  In  Section  1.5  we  present 
the  history  and  the  work  conducted  by  others  In  this  field. 
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1.2  Channel  and  Source  Models 


In  Chls  seccion  we  will  present  Che  channel  and  source  models  Chat 
will  be  used  In  this  thesis. 

1.2.1  Channel  Hodel 

The  channel  Is  assumed  to  be  slotted,  that  Is,  the  channel  Cine  Is 
divided  Into  equal  segments  called  slots.  The  length  of  each  slot  equals 
Chat  of  a packet,  and  It  Is  assuned  chat  a source,  which  Is  synchronized  to 
the  channel  time,  transmits  a packet  within  only  one  slot.  Furthermore, 

Che  channel  Is  such  Chat  Che  sources  can  determine  whether  there  are  zero, 
one,  or  multiple  packets  In  any  one  slot.  Multiple  packets  per  slot 
correspond  to  a collision  and  under  such  circumstances  no  one  gets  through. 
If,  however,  a packet  does  not  collide  %d.th  other  packets,  then  It  Is  assumed 
Chat  the  S/N  la  high  enough  or  enough  forward  error  correction  Is  applied  so 
chat  the  packet  Is  successfully  transmitted.  This  last  asstmiptlon  Is  made  so 
as  to  allow  us  to  focus  on  the  multi-access  properties  of  the  channel. 

1.2.2  Source  Models 

Ve  will  consider  two  source  models  In  this  thesis.  They  will  be 
designated  as  the  Poisson  and  the  finite  source  models. 

1.  Poisson  Source  Model 

The  Poisson  source  model  asstnes  the  existence  of  an  Infinite  number  of 
Independent  sources  that  collectively  generate  k packets  per  slot,  where  k 
Is  a Poisson  random  variable  %rlth  constant  mean  X.  A source  can  have  at 
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nose  tvo  packets,  one  that  has  undergone  a collision  and  Is  In  the  process 
of  being  retransmitted  and  one  that  may  have  arrived  since  the  collision 
of  the  first  packet  occurred.  If  a second  packet  arrives.  It  Is  not  trans- 
mitted until  after  the  first  Is  successfully  transmitted.  This  last  assump- 
tion will  be  discussed  In  more  detail  In  Section  1.3  ^Aere  the  algorithms 
are  considered. 

11.  Finite  Source  Model 

Here  we  assume  that  there  are  2 Independent  sources.  This  model  Is 
similar  to  the  Poisson  In  that  a source  can  have  at  moat  one  packet  In  the 
process  of  being  transmitted  or  retransmitted  and  at  most  one  waiting  to  be 
processed.  If  a source  has  at  most  one  packet  then  the  probability  that  It 
will  receive  a packet  In  the  next  round  trip  Interval  Is  constant  and  It 
is  given  by  p.  It  can  be  shown  that  the  Poisson  source  model  Is  the 
limiting  ease  of  the  finite  source  model.  That  is.  If  we  let  p2**  « constant 
and  let  It***  then  the  finite  model  approaches  the  Poisson  model. 
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1.3  Conflict  Reaolvlng  Algorithms 


Here  we  will  consider  the  third  component  of  the  multiple  access  system, 
the  conflict  resolving  algorithm.  This  section  is  organized  Into  four 
subsections  as  follows.  In  Subsection  1.3.1,  we  first  consider  the  TEHA  and 
Aloha  protocols,  the  two  algorithms  presently  In  use,  and  then  present 
a third  alternative,  the  tree  algorithm.  In  Subsections  1.3.2  and  1.3.3 
we  expand  upon  the  tree  algorithm  and  In  1.3.4  we  present  an  Information 
theoretic  discussion  of  the  multiple  access  channel. 

1.3.1  TDMA.  Aloha  and  the  Tree  Algorithm 

There  are  three  basic  techniques  for  accessing  a ccmnunlcatlons 
channel  In  the  time  domain,  the  TDMA,  the  Aloha  and  the  tree  algorithm.  In 
the  TDMA  scheme,  contention  Is  avoided  by  allocating  a portion  of  the  channel 
to  each  of  the  sources.  Although  TDMA  Is  effective  In  situations  where  the 
number  of  sources  Is  small  and  the  message  lengths  are  long,  it  suffers  from 
low  throughput  and  large  delays  when  the  nuiri>er  of  sources  is  large  and  the 
duty  cycle  Is  short. 

In  the  Aloha  algorithm,  when  a source  has  a new  packet.  It  transmits 
It,  and  then  listens  to  the  channel  to  determine  whether  or  not  the  packet 
collided  with  packets  from  other  sources.  If  a packet  collision  Is  detected, 
then  the  source  retransmits  the  packet  at  a randosdy  selected  time.  The 
retransmission  takes  place  at  a randosdy  selected  time  so  that  conflicting 
packets  will  not  surely  collide  again.  It  has  been  shorn  (see  Section  1.5) 
that  when  the  sources  satisfy  the  Poisson  source  model,  then  the  maximum 
throughput  for  the  Aloha  system  Is  1/e.  However,  a multi-access  Aloha  type 
system  Is  unstable,  and  It  eventually  overflows.  Therefore,  although  the 
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delay  and  throughput  properties  might  be  satisfactory  In  the  short  term 
they  are  quite  poor  when  observed  over  a long  Interval  of  time. 

There  are  numerous  of  dynamic  reservation  schemes  In  the  literature 
(see  Section  1.5).  Most  of  these  schemes  use  two  channels;  a reservation 
and  a data  channel.  Before  data  can  be  transmitted,  the  data  channel  Is 
reserved  by  TOMA  In  Aloha  techniques  on  the  reservation  channel.  It  follows 
then,  that  the  disadvantages  of  the  TIMA  and  Aloha  algorithms  apply  to  the 
reservation  channel. 

The  tree  algorithm  that  we  are  about  to  Introduce  has  (when  used  In 
conjunction  with  Poisson  source  model)  a maximum  average  throughput  of 
.430  packets/slot.  Is  stable  In  that  all  the  moments  of  the  delay  are  finite 
If  the  arrival  rate  Is  less  than  .430  packets/slot  and  It  has  good  delay 
properties.  (The  results  of  the  analysis  are  presented  In  more  detail  In 
Section  1.4.)  Below  In  Table  1.3. 1.1  we  present  a qualitative  comparison 
of  the  three  algorithms  when  the  number  of  sources  Is  large  and  the 
message  lengths  sre  short. 


TDMA 

Aloha 

Tree  Algr. 

Delay 

Poor 

Fair 

Good 

Throughput 

Poor 

Good 

Good 

Stsblllty 

Good 

Poor 

Good 

Table  1.3. 1.1.  A Qualitative  Comparison  of  TIMA,  Aloha  and  Tree  Algorithm 
when  the  number  of  sources  Is  large  and  the  message  lengths 
are  short. 
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It  Is  easiest  to  Introduce  the  tree  algorithm  by  an  example.  However, 
before  we  do  this,  we  need  to  make  the  following  definitions.  Figure  1.3. 1.1 
should  be  helpful 

Root  node  - the  initial  node  of  a tree;  in  Fig.  1.3. 1.1  It  Is  n^^. 

Depth  of  a node  - correponds  to  the  tier  at  which  the  node  Is  found. 

The  root  node  Is  at  depth  zero. 

Degree  of  a node  - Is  the  number  of  branches  that  emanate  from  a node. 

Subtree  - the  subtree  whose  root  node  Is  n^^^ . 

Note  that  in  a binary  tree,  j corresponds  to  the  particular 
node  of  depth  1,  and  that  there  are  2^  nodes  of  depth  1. 

Symmetry  - a tree  Is  symmetric  if  all  nodes  of  equal  depth  have  equal 
degrees . 

Example;  The  Binary  Tree  Algorithm 

Let  there  be  16  sources  {Sq,  S^^,  ...,  S^^^}  and  let  each  correspond  to 
a leaf  of  a 16-leaf  binary  tree  as  shown  In  Fig.  1.3. 1.1.  The  tree  may  be 
considered  as  an  addressing  procedure,  in  which  each  source  has  a 4-blt 
address  depending  on  Its  location  on  the  tree.  In  Fig.  1.3. 1.1  we  also 
present  the  slotted  satellite  time.  For  convenience,  we  assume  that  the 
round  trip  delay  Is  zero;  the  effect  of  a nonzero  round  trip  delay  is 
considered  in  Section  1.3.3.  Note  that  the  slots  are  paired  and  that  a slot 
pair  Is  designated  by  There  Is  a relationship  between  the  subscripts  of 

the  nodes  and  those  of  the  slot  pairs;  as  we  will  see,  T^^  transmits  Its 
packets  In  . 

Now  assume  that  no  collisions  have  occurred  until  the  beginning  of  SL^^, 
when  sources  S^,  82,  S^,  Sg,  and  each  has  a packet  to  transmit.  Then 
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SL22  I SLl^  I 


I SLqo  I SL|o  I SL20  I SL^l  I 

y 


Figure  1.3. 1.1  An  Exaaple  of  the  Binary  Tree  Algorlttaai 
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beginning  with  where  the  first  contention  arises,  the  tree  algorithm 

takes  the  following  steps  In  the  designated  slot  pairs. 

SLqq  Sources  In  transmit  their  packets  In  the  first  slot  of 

SLqq,  and  the  sources  In  transmit  theirs  In  the  second 
slot.  This  results  In  two  collisions,  one  among  S^,  $2,  and 
S^  and  the  other  between  Sg  and  Since  there  was  at 

least  one  collision  In  any  new  packets  that  arrive  are 

not  transmitted  until  that  contention  Is  resolved. 

SL^q  Since  there  was  a collision  In  the  sources  In  T^q  are 

divided  In  half  and  the  packets  In  T2Q  and  T22  are  transmitted 
In  the  first  and  second  slots  of  respectively.  This 

results  In  a collision  between  S^  and  S2  and  to  a successful 
transmission  by  S^. 

SL2g  Since  there  was  a collision  In  T2Q,  T^g  and  T^^  transmit 

their  packets  In  the  first  and  second  slots  of  SL2g, 
respectively.  This  results  In  two  successful  transmissions 
by  Sg  and  $2. 

SL^^  Since  there  was  a collision  In  ^22  and  T22  transmit 

their  packets  Is  succession.  This  results  In  a collision 
between  Sg  and  S^g  In  the  first  slot  and  no  transmission 
In  the  second. 
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SL22  Since  there  was  a collision  in  T22>  and  transmit. 

This  results  is  two  successful  transmissions  by  Sg  and 

SL^  Since  all  the  sources  were  Involved  in  some  transmission  in 

uU 

which  no  collisions  occurred,  we  know  that  the  original 
contention  has  been  resolved.  Any  new  packets  that  may  have 
arrived  to  during  this  conflict  resolution  interval,  are 

transmitted  in  the  first  slot  of  SLgg,  and  packets  that  arrived 
to  are  transmitted  in  the  second  slot.  The  process 
continues  on,  as  described  above. 

Note  that  in  this  example  we  used  10  slots  to  transmit  5 packets.  Next  we 
will  state  and  discuss  the  binary  tree  algorithm. 


The  Binary  Tree  Algorithm 

Let  each  source  correspond  to  a leaf  in  a binary  tree.  If  the  number 
of  sources  is  infinite,  then  the  tree  extends  to  infinity.  The  slots  are 
paired  into  odd  and  even  slots,  and  until  a collision  occurs  the  sources  in 
T^g  transmit  their  packets  in  the  odd  slots,  whereas,  the  sources  in  T^^ 
transmit  theirs  in  the  even  slots. 

How  let  T^^  and  T^2  variables  and  assume  that  no  collision  occurred 

up  to  the  beginning  of  the  present  pair  of  slots.  Thus,  the  binary  tree  algor- 
ithms is  as  follows: 


1. 

2. 


■ ■'10=  ^2 
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T^2  transmits  in  the  first  slot  of  the  present  pair  of  slots,  and  T^2 
transmits  in  the  second. 
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3.  If  any  collisions  occur  In  the  preceding  step,  then 

a.  Until  they  are  resolved,  no  new  packets  are  transmitted. 

b.  Resolve  the  first  collision  before  resolving  the  second. 

A collision  In  (1  ■ 1,2)  Is  resolved  by  dividing  In  two  halves 

(say  A and  B),  letting  - A,  T^2  ~ ^ then  repeating  steps  2 and  3. 

This  algorithm  Is  equivalent  to  the  following  tree  search.  Beginning 
with  the  root  node  and  at  each  succeeding  node,  one  asks  whether  there  are 
zero,  one,  or  more  than  one  packet  In  each  of  the  two  emanating  branches. 

If  the  answer  for  any  of  the  branches  Is  more  than  one,  then  proceed  to  the 
base  nodes  of  those  branches  and  repeat  the  question.  This  continues  until 
all  the  leaves  are  separated  Into  sets  such  that  each  set  contains  at  most 
one  packet.  Note  that  the  number  of  slots  used  In  a conflict  resolution 
Interval  equals  twice  the  number  of  nodes  visited. 

The  tree  search  may  be  carried  out  In  one  of  two  ways,  serially  or  In 
parallel.  In  the  serial  search,  two  branches  transmit  their  packets  In  two 
consecutive  slots,  and  the  results  of  those  two  transmissions  are  resolved 
before  another  two  subtrees  are  allowed  to  transmit.  In  the  parallel  search, 
all  the  branches  at  some  depth,  whose  parents  have  had  a collision,  transmit 
their  packets  In  consecutive  slots.  It  should  be  clear  that  the  number  of 
slots  needed  to  process  any  particular  sec  of  sources  Is  the  same  for  both 
schemes. 

The  question  arises  as  to  whether  some  other  tree  besides  the  binary 
tree  might  not  be  more  efficient.  In  that  It  requires  fewer  slots  to  resolve 
a conflict.  The  determination  of  the  optimum  tree  Is  one  of  the  major  problems 
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that  Is  solved  in  this  thesis.  More  specifically,  we  will  develop  the 
optimum  dynamic  tree  algorithm.  In  which  the  tree  that  Is  used  to  process 
a conflict  Is  that  which  minimizes  the  average  number  of  slots  used,  given 
the  history  of  the  transmission  process;  more  will  be  said  about  this  In 
Section  1.4. 

The  control  for  the  tree  algorithm  may  be  centralized  or  distributed. 

If  a central  control  Is  used,  then  an  observer  observes  the  transmission 
process  and  notifies,  via  a feedback  channel,  the  sources  that  may  transmit 
In  the  next  slot.  If,  on  the  other  hand,  a distributed  control  Is  used, 
then  all  the  sources  must  observe  the  channel  and  each  must  execute  the 
algorithm  by  Itself. 

We  will  present  an  Interesting  variation  to  the  above  algorithm,  before 
concluding  this  section.  Step  3a  of  the  algorithm  may  be  changed  somewhat, 
so  that  some  packets  that  arrive  after  a collision  occurs  may  be  transmitted 
before  that  Initial  contention  Is  resolved.  In  the  example  above,  for 
Instance,  packets,  that  arrive  to  T22  before  could  be  transmitted  In 

and  they  would  be  treated  similarly  to  the  way  they  would  have  been,  had  they 
arrived  In  This  variation  to  the  binary  tree  algorithm  appears  to  be 

more  efficient  then  the  original  algorithm  Itself.  However,  It  Is  more 
complicated  and  we  have  not  analyzed  It. 

1.3.2  Deterministic  and  Random  Source  Addressing 

In  the  tree  algorithm,  as  It  la  given  in  the  preceding  section,  each 
source  is  preaaslgned  detemlnlstlcally  an  address,  l.e.,  a position  on  the 
tree.  A variation  to  the  deterministic  addressing  which  has  certain  Imple- 


mentation  advantages.  Is  the  random  address  assignment  scheme.  Here,  the 
sources  are  not  preassigned  addresses  on  the  tree.  However,  when  a collision 
occurs,  a contending  source  conducts  the  tree  search  by  Independently  and 
with  equal  probability  deciding  to  take  the  upper  or  the  lower  branch  that 
emanate  from  each  node,  beginning  with  the  root  node.  In  general.  If  any 
collisions  occur  at  nodes  at  one  tier,  the  contending  sources  move  to  the  next 
tier  by  randomly  deciding  which  branch  to  take.  This  process  continues  until 
no  collisions  occur.  As  can  be  seen,  the  objective  of  this  search  Is  the  same 
as  that  of  the  deterministic  addressing  scheme;  It  Is  to  divide  the  sources 
Into  sets  such  that  each  set  contains  at  most  one  active  source. 

Note  that  In  the  random  addresssing  scheme,  the  tree  is  Infinite  whether 
the  number  of  sources  Is  finite  or  Infinite.  Therefore,  It  can  be  shown  that 
when  the  number  of  sources  Is  finite,  random  addressing  Is  slightly  Inferior 
to  deterministic  addressing  In  terms  of  delay  and  throughput.  However,  when 
the  number  of  sources  Is  Infinite,  such  as  In  the  Poisson  model,  then  the  two 
schemes  are  Identical.  To  see  this,  note  that  in  random  addressing,  the 
sequential  random  choice  decisions  that  each  contending  source  makes  are 
equivalent  to  allowing  each  of  the  contending  sources  to  choose  an  Infinite 
dimensional  address,  and  then  to  execute  the  tree  algorithm  as  It  Is  given  In 
the  preceding  section.  The  only  difference,  therefore,  between  the  two 
schemes  Is  the  way  the  addresses  are  chosen.  Note,  however,  that  the  statistics 
of  the  addresses  of  the  contending  sources  are  the  same  under  both  schemes. 

This  Is  so  because  In  one  scheme,  the  contending  sources  randomly  choose  them 
and  In  the  other  the  addresses  are  preassigned  but  the  sources  are  chosen 
randomly.  Since  the  address  statistics  are  Identical  and  since  the  algorithm 
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Is  Che  same  In  both  cases,  we  conclude  chat  in  the  Poisson  source  model  the 
throughput,  delay  and  stability  properties  are  the  same  for  both  addressing 
procedures . 

The  Implementation  advantage  of  the  random  addressing  scheme  Is  due  to 
the  fact  that  the  addresses  are  not  preassigned.  Therefore,  sources  may  enter 
or  leave  the  system  with  much  greater  ease.  Note  Chat  the  tree  search  for 
the  random  addressing  scheme  may  be  carried  out  serially  or  In  parallel,  and 
chat  Che  tree  need  not  be  binary. 

1.3.3  Sample  Transmission  Process  and  Definitions 

Here,  we  will  define  some  of  the  quantities  Chat  we  will  be  using, 
and  we  will  consider  the  effects  of  the  round  trip  delay  by  presenting  a 
sample  transmission  process  for  the  serial  search  binary  tree  (SSBT)  algorithm. 

In  the  execution  of  the  SSBT  algorithm,  two  branches  transmit  their 
packets  In  a pair  of  consecutive  slots;  following  this,  no  action  Is  taken 
until  the  results  of  these  two  transmissions  are  received,  when  two  more 
branches  are  allowed  to  transmit.  An  example  of  such  a transmission  process. 

In  which  the  round  trip  delay  equals  four  slots.  Is  illustrated  In  Fig. 

1.3. 3.1.  In  this  example,  since  one  algorithm  uses  only  1/3  of  the  channel 
capacity,  one  nay  either  divide  the  sources  Into  three  groups  and  process 
each  group  Independently  on  1/3  of  the  channel  by  a tree  algorithm  or  use 
one  tree  algorithm  on  1/3  of  the  channel  to  reserve  the  ocher  2/3 's  In  a 
dynamic  reservation  scheme. 

In  any  case,  we  will  focus  our  attention  on  one  algorithm.  Therefore, 
we  %rlll  assume  chat  our  channel  consists  only  of  those  slots  that  are  used 
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arrival  slots 


Figure  1.3. 3.1  A Seaple  Transnlsalon  Proceas 
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by  that  algorithm.  In  Fig.  1.3. 3.1  the  channel  that  we  are  Interested  In 
Is  that  composed  of  SL^,  SL2»  SL^,  ...,  etc. 

At  this  point,  we  will  make  several  definitions. 

Algorithm  Step  - It  consists  of  the  transmissions  taken  In  a pair  of 
slots,  the  observation  of  the  results  of  those  transmissions,  and 
the  decision  as  to  what  action  to  take  In  the  next  pair  of  slots. 

The  time  span  of  a step,  therefore,  equals  the  round  trip  delay 
plus  the  length  of  two  slots.  In  Fig.  1.3. 3.1,  a step  equals  all 
the  actions  taken  from  the  beginning  of  SL^  to  the  beginning  of  SL^. 

Epoch  - an  Interval  of  conflict  resolution.  If  a conflict  exists;  other- 
wise, It  Is  a pair  of  slots.  In  Fig.  1.3. 3.1,  [SL^,  SL2,  SL^],  [SL^l, 
(SL^,  SLg]  are  three  consecutive  epochs. 

Ij  - the  length  of  the  J'th  epoch  In  algorithmic  steps. 

hj  - the  number  of  slots  used  In  the  J'th  epoch.  This  equals  21^. 

A - the  average  of  the  total  nusd>er  of  packats  arriving  In  one  slot. 

Vij  - the  average  number  of  packets  arriving  In  the  J'th  epoch.  In  the 

text.  It  Is  shown  that  for  the  Poisson  source  model  ■ Ahj. 

5 - packet  delay;  l.e.,  the  time  spent  In  the  system  by  a packet. 

Normally  we  will  express  the  delay  In  terms  of  algorithm  steps. 

average  delay  - the  delay  of  a randomly  chosen  packet. 

average  throughput  - the  fraction  of  the  slots,  over  a very  long  Interval, 
that  contain  exactly  one  packet  each. 

stability  - the  systam  Is  K'th  order  stable  If  the  K'th  moment  of  the 
delay  Is  finite. 
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In  terms  of  the  preceding  definitions.  In  packet  switching  networks, 
packets  that  arrive  In  one  epoch  are  transmitted  In  the  following  epoch. 

Note,  however,  that  since  only  a portion  of  the  channel  Is  used  by  the  algor- 
ithm, only  those  packets  arriving  in  slots  corresponding  to  those  used  by 
the  algorithm  In  one  epoch  will  be  processed  by  the  same  algorithm  In  the 
following  epoch.  In  the  example  of  Fig.  1.3. 3.1,  the  packets  that  arrive  In 
slot  pairs  SLj^,  SL^,  and  SL^  are  transmitted  In  SL^.  Those  that  arrive  In 

SL,'  are  transmitted  to  SL,  and  those  that  arrive  In  SL'  and  SLI  are  transmitted 
4 5 5 6 

In  SL^. 

An  Interesting  and  a very  useful  property  of  the  transmission  process 
Is  that  under  the  finite,  as  well  as  the  Polsslon  source  model,  Is  an 

embedded  Harkov  chain.  That  Is,  can  be  considered  to  be  the  state  of  a 
Markov  chain  after  the  j'th  transition.  To  see  this,  first  note  that 
given  Vj , Is  Independent  of  the  transmission  process  up  to  the  end  of  the  J ' th 
epoch.  (Vj  Is  the  number  of  packets  arriving  In  the  j'th  epoch.)  This 
observation  coupled  with  the  following  probabilities  (which  are  developed 
In  the  text)  prove  the  Markovian  property. 

V -(2li.) 

(2Xip  J e ^ 

p(Vj  ^ V I “o***^* 

• l^j) 

p(a  packet  arrives  to  a particular  source  In  the  J'th  epoch|l^,£.^_^, .. .) 

*1 

■ 1 - (1-p)  for  the  finite  model. 
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1.3.4  An  Information  Theoretic  Approach  to  Multl-Accesslng 

As  we  will  see,  the  naxlsrum  average  throughput  of  the  optimum  dynamic 
tree  algorithm  la  .430  packets/slot  when  the  sources  satisfy  the  Poisson 
model;  whereas,  the  maximum  throughput  for  the  Aloha  system  Is  1/e  * .368. 

Since  the  performance  of  these  two  algorithms  Is  different,  the  question 
arises  as  to  whether  there  are  any  other  algorithms  which  are  better,  or, 
more  Interestingly,  what  Is  the  maximum  possible  throughput  to  a particular 
multiple  access  system  tinder  all  possible  algorithms.  We  have  attempted  to 
answer  these  questions  but  were  unsuccessful.  However,  our  approach  to  this 
problem  offers  some  Insight  Into  the  multl-accesslng  problem  and  since  It  might 
be  the  basis  for  further  research,  we  outline  It  here. 

We  have  pointed  out  In  the  preceding  section  that  the  basic  problem 
In  a multiple  access  system  Is  the  resolution  of  the  conflicts.  So,  let  us 
examine  this  a little  more  carefully.  Let  us  assume  that  we  have  a set  of 
Independent  sources,  and  v of  them  are  active,  l.e. , have  packets  to  transmit. 
Furthermore,  let  the  probability  measure  on  V be  p(v). 

Now,  note  that  the  conflict  Is  resolved  Iff  the  sources  are  subdivided 
Into  sets  such  that  each  set  contains  at  isost  one  active  source.  Therefore, 
this  partitioning  of  the  sources  must  be  the  objective  of  any  conflict 
resolving  algorithm,  whether  It  be  the  Aloha,  the  tree,  or  any  other  that  might 
be  proposed  In  the  future.  In  other  words,  the  execution  of  the  algorithm 
must  supply  enough  Information  so  that  one  can  partition  the  sources  Into 
sets  such  that  each  set  contains  at  most  one  active  source.  Let  (source) 
be  the  minimum  average  Information  required  to  do  this. 
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Next,  by  observing  the  contents  of  a slot,  we  lear’'i  whether  there  are 

zero,  one,  or  more  Chan  one  active  sources  In  the  set  that  la  transmitting  In 

chat  slot.  Let  H (trans.)  be  the  maximum  average  Information  (maximized 

over  all  partitions)  that  can  be  obtained  from  any  one  slot;  certainly, 

(trans. ) ^ log2(3).  Therefore,  the  minimum  average  number  of  slots 

required  to  resolve  Che  conflict  under  any  algorithm  must  be  equal  to  or 

greater  than  H . (source)/H  (trans.). 
mxn  max 

We  did  not  proceed  beyond  this  formulation  because  we  were  unable  Co 

obtain  good  bounds  to  H . (source). 

nun 


1.4  Thesis  Outline  and  Results 


The  thesis  Is  organized  Into  four  chapters.  Chapter  1 Is  this  Intro- 
duction. The  analysis  Is  carried  out  In  Chapters  2,  3 and  4 and  In  the 
appendices  that  accompany  these  chapters.  Here  we  will  consider  the 
objectives  and  results  of  each  of  these  chapters. 

In  Chapter  2,  we  analyze  the  serial  search  binary  tree  algorlthn  when 
It  Is  used  In  conjunction  with  the  Poisson  source  model.  There  are  two 
main  results  In  this  chapter.  The  first  concerns  the  average  delay  vs. 
average  throughput  trade-off;  we  obtain  upper  and  lower  bounds  to  the 
average  delay,  as  a function  of  X.  These  bounds  are  given  In  Eqs.  (2. 2. 4. 2) 
and  (2. 2. 4. 6)  and  they  are  Illustrated  In  Fig.  2. 2. 4.1.  It  Is  shown  that 


these  results  may  be  Interpreted  as  average  delay  vs.  average  throughput, 
and  we  show  that  the  maximum  average  throughput  Is  .347  packets/slot. 
Furthermore,  we  show  that  It  Is  possible  to  obtain  a throughput  of  up  to 
.430  packets/slot,  but  only  for  a limited  time. 

The  second  main  result  of  Chapter  2 concerns  the  stability  of  the  binary 
tree  algorithm.  Here,  we  prove  that  If  X < 1/3  packets/slot,  then  all  the 
moments  of  the  delay  are  finite.  However,  we  observe  that  this  Is  an  overly 
conservative  result,  and  point  out  that  Indications  are  that  all  the  moments 
of  the  delay  are  finite  for  X < .347  packets/slot . 

In  Chapter  3,  we  determine  and  analyze  the  optimum  dynamic  tree  algorithm 
and  examine  a suboptimum  algorithm  that  has  certain  Implementation  advantages. 
This  algorithm  Is  called  optimum  dynamic,  because  the  tree  Is  allowed  to  vary 
from  epoch  to  epoch  optimally  depending  on  the  traffic.  The  source  model 
that  Is  assuzMd  In  this  chaptar  Is  the  Poisson. 


More  specifically,  we  show  that  the  tree  which  minimizes  the  expected 
number  of  slots  needed  to  process  V packets,  where  V Is  a Poisson  random 
variable.  Is  binary  everywhere  except  for  the  root  node  whose  degree  g^ 
depends  on  w and  Is  given  by  Eq.  (3. 1.0. 2). 

Ue  point  out  that  y « hX;  therefore,  by  observing  the  number  of  slots 
In  the  preceding  epoch,  we  can  determine  y and  hence  from  Eq.  (3. 1.0. 2) 
determine  the  optimum  tree  to  be  used  In  the  next  epoch.  In  order  to 
simplify  the  analysis,  we  restricted  the  degree  of  the  root  node  to  be  2r. 
Thus,  the  root  node  corresponds  to  r algorithm  steps  and  each  of  the  other 
nodes  corresponds  to  one  step.  For  this  dynamic  algorithm  we  have  obtained 
upper  and  lower  bounds  to  the  average  delay  as  a function  of  A.  These  results 
are  given  In  Eqs.  (3.3.2.13)  and  (3.3.2.14)  and  displayed  In  Fig.  3. 3. 2. 2. 

We  also  show  that  the  maximum  average  throughput  Is  .430  packets/slot,  and 
we  prove  that  all  the  moments  of  the  delay  are  finite  for  X < .430  packets/ 
slot. 

In  Chapter  3,  we  also  consider  the  more  easily  Implemented  algorithm  In 

which  the  root  node  degree  Is  restricted  to  be  2 (K  > 0) , and  all  other  nodes 

* 

are  binary.  Subject  to  the  above  constraints,  first,  we  determine  K , the  K 
which  minimizes  E{hj^^|yj .K};  this  Is  given  In  Eq.  (3. 1.0. 3).  Next  we 
determine  upper  and  lower  bounds  to  the  average  delay  vs.  X;  these  are 
given  in  Eqs.  (3.4.2.13)  and  (3.4.2.14)  and  displayed  In  Fig.  3. 4. 2.1.  The 
maximum  average  throughput  Is  shown  to  be  greater  than  .420  but  less  than 
.430  packets/slot. 
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In  Chapter  4 we  consider  the  static  and  dynamic  tree  algorithms  when 
they  are  used  In  conjunction  with  the  finite  source  model.  For  the  static 
binary  algorithm  we  obtain  an  average  delay  upper  bound  vs.  average  through- 
put lower  bound  curve.  This  Is  shown  In  Fig.  4. 2. 2. 2 for  64  sources.  The 
maximum  throughput  for  the  64  source  model  Is  .507  packets/slot. 

The  dynamic  tree  was  restricted  to  be  binary  everywhere  except  for 
the  root  node  whose  degree  was  restricted  to  be  8q*2  . This  algorithm 
was  optimized,  as  before,  over  K and  the  results  are  given  by  Eqs.  (4. 3. 1.1) 
and  (4. 3. 1.6).  An  Interesting  observation,  on  the  optimum  dynamic  algorithm. 
Is  that,  under  low  traffic.  It  Is  Identical  to  the  binary  tree  algorithm. 
However,  as  the  traffic  Increases  this  algorithm  adaptively  changes  to  a tree 
that  has  only  one  node  with  2 branches,  which  Is  recognized  to  be  the  TDMA 
protocol.  (Note  2 equals  the  number  of  sources.) 

The  delay- throughput  characteristics  of  this  optimum  dynamic  tree  are 
determined,  and  are  Illustrated  In  Fig.  4.3.2.S  for  2 * 64.  The  maxlmtan 

average  throughput  for  this  algorithm  Is  one  packet/slot.  Chapter  4 Is 
concluded  with  a theorem  proving  that  the  average  delay  of  the  optimum  dynamic 
tree  Is  less  than  or  equal  to  that  of  the  TDHA  protocol. 
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1.5  History 

The  types  of  multiple  access  problems  that  have  been  considered  by 
others  fall  Into  two  main  categories,  circuit  switching  and  packet  switching. 

A third  category  contains  the  various  dynamic  reservation  techniques,  but 
this  can  be  considered  to  be  a hybrid  of  the  first  two  classes. 

With  circuit  switching,  the  channel  Is  partitioned  and  allocated  (leased) 
before  hand.  The  sharing  of  the  channel  this  way  has  been  accomplished 
either  by  FDMA  (Frequency  Domain  Multiple  Access),  or  by  TDMA  (Time  Domain 
Multiple  Access)  or  by  a combination  of  FDMA  and  TDMA.  These  techniques 
have  been  the  only  ones  that  were  In  use  up  to  about  1970  [1]  and  they  were 
quite  effective  for  the  communication  needs  of  that  time.  If,  however,  there 
Is  a large  number  of  sources  that  do  not  require  continuous  full  use  of  the 
channel  or  If  the  data  Is  bursty,  l.e.,  the  ratio  of  peak  data  rate  to 
average  data  rate  Is  high,  then  circuit  switching  can  lead  to  long  transmission 
delays  and  Inefficient  use  of  the  channel. 

In  recent  years  attention  has  shifted  to  packet  switching  forms  of 
multiple  access.  With  this  form  of  accessing,  as  was  pointed  out  earlier, 
the  transmitter  formats  the  data  Into  a packet  of  constant  length  along 
with  source  and  destination  addresses  and  error  detecting  bits.  The  source 
then  transmits  the  packet  to  all  the  receivers  Including  the  desired  one. 

If  a packet  la  destroyed  In  transmission,  the  originating  source  learns 
about  It  either  through  a feedback  channel  or  by  listening  to  Its  own 
transmission.  When  a source  determines  that  Its  packet  has  been  destroyed 
It  retransmits  that  packet  at  a randomly  selected  time.  Most  of  the  wrk 
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that  has  been  done  to  date  on  packet  switching  broadcast  systems  assumes 
that  the  channel  is  noise  free  and  chat  all  Che  packets  involved  in  a 
collision  are  destroyed. 

There  have  been  two  approaches  to  packet  switching  broadcast  systems; 
pure  Aloha  and  slotted  Aloha.  In  pure  Aloha,  packets  are  transmitted  or 
retransmitted  asynchronously;  in  slotted  Aloha  the  sources  are  synchronized 
so  that  packets  are  transmitted  in  phase.  In  dynamic  reservations  the 
sources  first  dynamically  reserve  the  channel  capacity  via  an  Aloha  or  a 
TCHA  channel  and  then  transmit  their  data. 

Pure  or  classical  Aloha  has  been  studied  by  Abramson  [1].  He  assumed  chat 
1)  Che  starting  times  of  the  packets  that  are  offered  to  the  channel  for  the 
first  time  from  all  the  sources  comprise  a Poisson  point  process,  and  2)  Che 
starting  times  of  all  the  packets  (new  plus  retransmitted  packets)  comprise 
another  Poisson  point  process.  Given  these  two  assumptions  which  taken  together 
imply  an  equilibrium  condition,  he  proves  that  the  capacity  of  the  system  is 
or  approximately  .184. 

Roberts  [2]  pointed  out  that  considerable  Improvanent  can  be  made  in  the 
capacity  of  the  Aloha  channel  by  synchronizing  the  sources  so  that  all  the 
packets  arrive  at  the  channel  in  phase.  He  proposed  that  the  channel  time  be 
divided  into  slots  and  sources  be  allowed  to  use  at  most  one  slot  per  packet. 

By  doing  this  he  showed  that  capacity  of  the  slotted  system  for  the  Poisson 
source  model  is  or  twice  that  of  pure  Aloha. 

Metcalf  [3]  considered  Aloha  systems  with  blocking  (that  is  a source  may 
not  generate  a second  packet  until  the  first  is  successfully  transmitted)  and 
exasdned  several  random  retransmission  policies.  More  Important  he  showed 
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that  Aloha  channels  may  be  unstable  In  that  the  number  of  blocked  sources  can 
become  very  large.  He  also  proposed  control  techniques  by  varying  the  retrans- 
mission probabilities. 

Klelnrock  and  Lam  [4],  [5]  quantified  and  extended  Metcalf's  results. 

They  modeled  the  slotted  Aloha  system  with  blocking  by  a Markov  Chain  whose  state 
corresponded  to  the  number  of  blocked  sources.  They  showed  that  the  slotted 
Aloha  channel  with  an  Infinite  but  Independent  population  (l.e. , the  Poisson 
source  model)  Is  unstable.  They  further  showed  that  If  the  number  of  sources  Is 
finite  then  the  Aloha  system  may  be  bistable  In  that  It  Is  possible  to  have  two 
stable  operating  points  - one  with  a small  number  and  one  with  a large  number  of 
blocked  sources.  Random  perturbation  In  the  channel  traffic  will  cause  the 
system  to  vacillate  between  these  two  stable  points.  They  also  proposed  con- 
trolling the  system  through  the  retransmission  probabilities  and  derived  an 
optimum  control  policy  based  on  exact  knowledge  of  the  state  of  the  channel. 
Carlelal  and  Heilman  [6]  also  modeled  the  Aloha  system  as  a Markov  chain  and 
examined  Its  bistable  behavior. 

Several  dynamic  reservation  schemes  have  been  proposed.  Basically,  these 
techniques  use  either  slotted  Aloha  or  TDMA  to  make  the  reservations.  Two 
protocols  that  use  Aloha  techniques  to  reserve  the  channel  are  Reservation- 
Aloha  Introduced  by  Crowther  et.al.  [7]  and  Interleaved  Reservatlon-Aloha  sug- 
gested by  Roberts  [8].  In  Reservatlon-Aloha  the  channel  slots  are  grouped  Into 
frames  that  are  at  least  a round-trip  delay  long.  A source  that  has  success- 
fully used  a particular  slot  In  one  frasM  has  access  to  that  slot  In  the  fol- 
lowing frsme.  If  a slot  Is  unused  then  It  Is  up  for  grabs  and  any  one  can  con- 
tend for  It  by  random  Aloha  techniques.  In  the  Interleaved  Reservatlon-Aloha 
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the  channel  Is  divided  Into  two  states.  Reservation  and  Aloha.  On  the  Reser- 
vation state  the  sources  attempt  to  reserve  the  Aloha  state  through  slotted 
Aloha  techniques.  As  traffic  Increases  the  percentage  of  Che  channel  being  In 
Che  Aloha  state  Increases  allowing  for  greater  utilization  of  the  channel.  In 
the  Interleaved  Reservation  Aloha  as  Che  ratio  of  message  length  to  reservation 
packet  length  Increases  the  maximum  throughput  approaches  lOOZ.  Roberts  also 
showed  (In  an  example  where  Che  average  message  length  was  27  tines  Chat  of 
Che  reservation  packet)  that  Che  transmission  delay  for  this  scheme  Is  better 
Chan  slotted  Aloha  at  large  throughput  but  not  as  good  for  throughput  less  chan 
20Z.  Crowther,  et.al.  do  not  present  an  analysis  of  their  scheme  but  one  would 
expect  a behavior  similar  to  Robert's  algorithm.  Stability  Is  not  considered  In 
either  paper  but  as  long  as  random  accessing  techniques  are  used  for  Che  reserva- 
tion, stability  Is  an  Issue. 

Binder  [9]  offers  a dynamic  reservation  scheme  chat  uses  TDMA  techniques. 
Essentially  what  he  proposes  Is  that  the  slots  be  grouped  Into  frames  and  each 
source  be  allocated  a slot  to  which  It  has  first  priority.  If  a slot  is  not 
used  by  Its  owner  then  it  Is  available  to  the  ocher  sources  on  a round-robin 
basis.  Unfortunately,  here  again  there  Is  no  analysis.  Limited  simulations 
Indicate  Chat  this  protocol  Is  better  chan  slotted  Aloha  at  high  traffic  and 
worse  at  low  traffic. 

In  the  work  chat  has  been  discussed  up  to  now  the  authors  assumed  chat  If 
more  than  one  packet  Is  transadcted  simultaneously  chan  a collision  occurs  and 
all  Che  packets  are  destroyed.  Roberts  [2]  points  out,  however,  that  this  need 
not  be  Che  case.  FM  receivers  will  crack  the  strongest  of  many  signals  as  long 
as  Che  next  strongest  la  down  by  1.5  to  3 dB.  He  examines  the  problem  where 
one  receiver  Is  surrounded  by  an  equal  density  population  of  equal  power  trans- 
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mltters  and  shows  that  It  Is  possible  for  the  capacity  to  Increase  to  .60  with 
FM  capture.  His  example  Is  for  ground  radio  systems,  but  one  should  be  able  to 
take  advantage  of  FM  capture  In  satellite  systems  where  the  sources  are  essen- 
tially equidistant  from  the  satellite  by  regulating  the  power  of  the  various 
transmitters. 

Another  technique  that  offers  considerable  Improvement  to  Aloha  systems 
that  have  a small  maximum  time  delay  Is  CSMA  (Carrier  Sense  Multiple  Access). 

In  CSMA  a terminal  with  a packet  to  transmit  first  listens  for  the  carrier  of 
other  users  to  determine  whether  the  channel  Is  busy.  If  It  Is  busy  then  the 
source  refrains  from  transmitting,  whereas  If  the  channel  Is  empty  of  other 
carriers  the  source  will  transmit  Its  packet  with  some  probability.  A colli- 
sion will  occur  If  the  time  between  transmissions  of  two  packets  Is  less  than 
the  time  delay  between  the  corresponding  sources.  In  case  of  a collision  the 
contending  sources  retransmit  at  randomly  selected  times.  Klelnrock  and 
Tobagl  [10],  [11]  have  examined  CSMA  and  have  shown  that  If  a)  there  are  no 
hidden  terminals,  l.e.,  every  terminal  can  hear  every  other  terminal,  b)  the 
time  to  detect  the  carrier  Is  zero,  and  c)  the  distance  between  terminals  Is 
not  large  than  considerable  Improvement  can  be  made.  For  example.  If  the  ratio 
of  the  maximum  time  delay  between  terminals  to  packet  length  la  .01  then  the 
channel  capacity  la  .86.  CSMA  can  be  quite  effective  In  ground  packet  switching 
channels,  but  It  is  ineffective  In  satellite  systems. 
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CHAPTER  2 


STATIC  BINARY  TREE  ALGORITIM  WITH  POISSON  SOURCE  MODEL 

2.1  Introduction 

In  this  chapter  we  will  exaalne  the  delay,  throughput  and  stability 
properties  of  the  static  binary  tree  algorithm  when  It  Is  used  In  conjunction 
with  the  Poisson  source  model.  We  will  restrict  the  analysis  to  the  serial 
search  and  to  the  deterministic  address  assignment.  As  was  pointed  out  In 
Chapter  1,  under  the  Poisson  assumption,  Che  random  and  deterministic  source 
address  assignment  schemes  are  Identical  In  terms  of  delay,  throughput  and 
stability.  However,  the  parallel  and  serial  search  algorithms  have  the  same 
throughput,  but  the  average  delay  of  the  parallel  search  scheme  Is  less  than 
that  of  the  serial. 

The  average  packet  delay  Is  defined  to  be  the  delay  that  a randomly 
selected  packet  undergoes,  and  It  Is  the  topic  of  Section  2.2.  There,  we 
obtain  upper  and  lower  bounds  to  the  mean  packet  delay  as  a function  of  the 
arrival  rate  X.  These  bounds  are  given  by  Eqs.  (2. 2. 4. 2)  and  (2. 2. 4. 6)  and 
are  displayed  In  Fig.  2. 2. 4.1. 

The  average  throughput  la  defined  to  be  the  fraction  of  slots  over  a 
very  long  Interval  of  time  that  contain  exactly  one  packet  each.  It  Is 
considered  In  Section  2.3,  where  we  argue  that  if  the  average  delay  Is 
finite,  then  the  average  arrival  rate  equals  the  average  throughput.  In 
that  section  we  also  show  that  the  maximum  average  throughput  is  .347 
packets/slot  and  that  It  Is  possible  to  attain  a throughput  of  up  to  .430 
packets/slot  but  only  for  a limited  time. 


-38- 


The  system  Is  considered  to  be  k'th  order  stable  If  the  first  k 
moments  of  the  delay  are  finite.  Stability  Is  considered  in  Section  2.4. 
There,  first,  we  show  that  If  X < .347  packets/slot,  then  the  system  is 
first  order  stable.  Secondly,  we  prove  that  If  X < 1/3  then  all  the  moments 
of  the  delay  are  finite.  This  result,  however,  seems  to  be  overly  conserva- 
tive, and  Indications  are  that  all  the  moments  of  the  delay  are  finite  If 
X < .347  packets/slot.  Accompanying  the  main  text  Is  Appendix  A2  where 
some  crucial  results  of  this  chapter  are  developed.  A detailed  outline  of 
Appendix  A2  Is  given  in  the  following  section  after  several  definitions  are 
made. 
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2.2  Average  Delay 


The  objective  of  this  section  Is  to  develop  a characterization  of  the 
average  packet  delay  In  a multi-access  communication  system  that  uses  the 
static  binary  tree  protocol.  More  specifically,  we  will  determine  upper  and 
lower  bounds  to  the  average  delay,  e{5},  that  are  functions  of  the  packet 
arrival  rate  X.  (See  Eqs.  (2. 2.4. 2)  and  (2. 2.4. 6)  and  Fig.  2. 2. 4.1.) 

The  unit  of  measure  for  the  delay  Is  the  time  between  two  successive 
steps  of  the  algorithm.  This  Is  constant  and  It  normally  equals  the  round 
trip  delay  plus  the  transmission  time  of  two  packets.  The  delay,  therefore, 
that  a packet  undergoes  Is  directly  proportional  to  the  number  of  steps  that 
are  executed  by  the  algorithm  from  the  time  when  a packet  arrives  to  the  time 
when  It  Is  received  correctly. 

The  expected  delay  Is  defined  to  be  the  delay  that  a randomly  chosen 
packet  will  undergo.  Note  that  Inherent  In  the  definition  of  E{5}  Is  the 
concept  of  random  Incidence,  this  will  play  a central  role  In  the  analysis. 

As  has  been  discussed  previously,  packets  that  arrive  In  one  epoch  (an 
epoch  Is  an  Interval  of  conflict  resolution  If  a conflict  exists  or  simply  a 
pair  of  slots  If  there  Is  no  conflict)  wait  until  that  epoch  ends  and  are 
transmitted  In  the  following  epoch.  With  this  In  mind,  several  definitions 
are  presented  below.  All  lengths  are  In  algorithmic  steps. 

■ the  epoch  In  which  a randomly  chosen  packet  arrives. 

€2  " Che  epoch  following  e^. 

y^  " the  length  of  This  Is  a random  Incidence  random  variable. 

• a random  variable  that  equals  the  length  Cj.  The  probability  that 
Ij  ■ b Is  the  steady  state  probability  that  an  epoch  of  length 


L occurs 
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■ upper  bound  to  E{A} 

■ lower  bound  to  E{A} 

^ 2 

■ upper  bound  to  E{A  } 

2 2 

■ lower  bound  to  E{A  } 

hj  ■ the  number  of  slots  used  by  the  algorithm  In  . 

Note  that  h^  * 21^  for  the  binary  tree  algorithm, 
dj  * the  time  spent  by  a packet  In  e^. 

Vj  ■■  the  number  of  packets  that  arrive  In  . 

X > the  average  number  of  packets  that  arrive  In  any  one  slot, 

pj  • the  average  number  of  packets  that  arrive  In  given  A^ . Note 

that  Uj  - E{Vj1Aj}  - 2XAj  - Xh^ . 

6 * the  packet  delay,  or  the  time  spent  In  the  system  by  a packet.  This 

quantity  equals  the  sum  of  d^  and  d^. 

■ lower  bound  E{6} 

■ upper  bound  E{6} 

As  we  proceed,  we  will  further  expand  upon  some  of  the  relationships  of 
the  above  quantities. 

Next  we  will  present  an  outline  of  the  analysis  that  follows.  In 
Section  2.2.1  we  develop  the  relstlonshlp  between  p^,  A^  and  X.  In  Section 

2.2.2  upper  and  lower  bounds  to  E{6}  are  developed  that  are  functions  of 

2 2 
E{A}  and  E{A  }.  In  Section  2.2.3  upper  and  lower  bounds  to  E{A}  and  E{A  } 

are  developed  that  are  functions  of  X.  Finally,  In  Section  2.2.4  we  coid>lne 

the  results  of  Sections  2.2.2  and  2.2.3  to  obtain  upper  and  lower  bounds  to 

E{6}  that  are  functions  of  X. 
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In  Appendix  A2  we  prove  some  of  the  more  tedious  but  still  very  crucial 

results  of  this  chapter.  This  appendix  Is  organized  Into  eight  sections.  In 

Section  A2.1,  we  derive  an  expression  for  and  In  A2.2  some 

Important  properties  of  are  developed.  In  Section  A2.3  an 

2 

expression  for  Is  derived  and  Its  properties  are  developed  In 

A2.4.  In  Section  A2.5  It  la  proved  that  ^ . In  A2.6 

E{dj_^^|yj}  Is  derived  and  In  A2.7  we  develop  the  relationship  between 
E{dj_|_j^|Vj|}  and  E{ | } . Finally  In  A2.8  we  obtain  upper  bounds  to 

E{e  ^ ^ ^|yj}/38. 


2.2.1  Relationship  Among  y^,  A^  and  X. 


As  defined  above,  X Is  the  average  number  of  packets  that  arrive  In 
any  one  slot,  y^  Is  the  average  number  of  packets  that  arrive  In  £j  and  A^  Is 
the  nuBber  of  steps  executed  by  the  algorithm  In  . For  the  binary  tree 
which  Is  being  considered  here,  A^  also  equals  the  number  of  nodes  visited  by 
the  algorithm.  It  follows  from  the  above  definitions  that  the  length  of 
in  seconds  Is  given  by  A, (T  4*  2t  ) where  T is  the  round  trip  delay  and  T 

j ^ 8 If  8 

Is  the  length  of  one  slot. 

Another  relationship  which  will  be  proved  here  Is  the  following: 


(2. 2. 1.1) 


( 


To  see  this,  first  note  that  the  binary  tree  algorithm  uses  two  slots  for 
each  step:  therefore 


"j  - “j 
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(2. 2. 1.2) 


4 


Next  note  that  the  total  number  of  packets  that  arrive  in  equals  the  sum 
of  the  packets  that  arrive  In  each  of  the  h^  slots,  and  since  the  expectation 
of  a sum  equals  the  sum  of  the  expectations.  It  follows  that 


hjX 


Equation  (2. 2. 1.1)  follows  from  Eqs.  (2. 2. 1.2)  and  (2. 2. 1.3). 


(2. 2. 1.3) 


2 

2.2.2  Characterization  of  E{ delay}  In  terms  of  E{&}  and  E{A  } 

In  this  section  we  will  develop  upper  and  lower  bounds  to  E{6}  that 
are  functions  of  E{i}  and  E{2.  }.  More  specifically  if  we  let  6^  and  5^^  be 
the  upper  and  lower  bounds  to  E{6},  respectively,  then  we  will  prove  that 


6 


u 


1.05 


EU^} 

eIIF 


+ .321 


(2. 2. 2.1) 


and 


Now  we  will  begin  with  the  analysis.  The  total  delay  that  a randomly 
chosen  test  packet  will  undergo  can  be  decomposed  Into  two  parts;  d^,  the  time 
spent  In  (the  epoch  In  which  it  arrived)  and  d^t  the  time  spent  in  £2 
(the  epoch  in  which  It  Is  transmitted). 

6 - d^  + d2  (2. 2.2.3) 
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and 


E{5}  ■ E{dj^}  + E{d2}  (2. 2. 2. 4) 

Next  we  will  derive  expressions  for  E{dj^}  and  E{d2}>  We  begin  with  E{dj^}. 

If  we  let  be  the  length  of  the  epoch  that  the  test  packet  entered, 
then  is  a random  incidence  random  variable.  It  is  well  kno%m  [Ref.  12, 
p.  149]  that  the  densltY  of  y^^  expressed  in  terms  of  the  densltY  of  as 

follows, 

■ E(ij)  ■ (2.2.2.5) 

Since  the  slot  in  which  the  test  packet  enters  can  occur  with  equal  proba- 
bllltY  anYwhere  in  we  have 

E{di}  • j RCyj^}  (2. 2. 2. 6) 

From  Eqa.  (2.2.2. 5)  and  (2. 2. 2. 6)  follows 


E{dj^} 


yi  Pi^<yi> 
“TTiT 


or 


E{dj^} 


1 11^ 

2 eTIT 


(2. 2. 2. 7) 
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The  expectation  of  62  can  be  written  as  follows: 

E{<i,)  • J E(d,|y  } p (y  ) 

El  '1 

(2. 2. 2. 8) 

and  substituting  Eq.  (2. 2. 2. 5)  Into  Eq.  (2. 2. 2. 8),  we  have 

(2. 2. 2. 9) 

Combining  Eqs.  (2. 2. 2. 4),  (2. 2. 2. 7)  and  (2. 2. 2. 9),  we  have 

the  following 

result. 

" 2 “eTIT  eTeT  1 V^‘^2l*'i^Pa-  ^*'1> 

1 ^ 

(2.2.2.10) 

In  Appendix  A2.6,  the  following  expression  for  E{d2|f'2^ 

derived. 

00 

E{d2|y)  - 1 + I 0(y/2^)(l  +-|D(y/2^)) 

(2.2.2.11) 

J-1 

where 

0(y)  - 1^-  ^ 

1-e  ^ 

(2.2.2.12) 

flO 

D(y)  - 1 2h(u/2h 

(2.2.1.13) 

1-0 

C(W)  - 1 - e"P  - ye'P 

(2.2.2.13a) 
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1 

and 


y - 2X1^  (2.2.2.14) 

Note  that  the  subscript  of  y has  been  dropped  for  convenience. 

Equations  (2.2.2.10)  through  (2.2.2.14)  is  as  far  as  we  can  go  in 
obtaining  an  exact  closed  form  expression  for  E{6}.  Since  this  expression 
of  E{6},  as  it  stands,  is  quite  complex,  we  will  turn  our  attention  to  deriving 
upper  and  lower  bounds. 

In  Appendix  A2.1,  the  following  expression  for  £{12!^^  is  derived. 

- 1 + 2D(y/2)  (2.2.2.15) 

and  in  Appendix  A2.7  it  is  proved  that 

E{d2|y}  < .55  E{42lw>  + -321  for  y ^ 0 (2.2.2.16) 

But  because  of  Eq.  (2.2.2.14)  it  follows  that 

E{d2|fi}  < .55  E{Jl2|ti>  + -321  (2.2.2.17) 

Nov  conblnlng  Eqs.  (2e2e2.10)  and  (2e2e2.17)  ve  have 

2 

1 i EO)  EOJ  ^ + -321  (2.2.2.18) 

^1 
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It  can  easily  be  shown  Chat 


*■1 


-i  E{£i^}  + j E{l^h 


(2.2.2.19) 


But  since 


2 

£ we  have  from  Eqs.  (2.2.2.18)  and  (2.2.2.19) 


E{6}  < 1.05 


E{£2} 

“eTIT 


+ .321 


(2.2.2.20) 


This  Is  Che  desired  upper  bound.  Next  we  will  derive  an  equivalent  lower 
bound  to  E(($}. 

In  Appendix  A2.7  It  Is  also  shown  that 

E{d2|Uj^}  > Y £{^2!^^  for  y > 0 (2.2.2.21) 


Here  again,  from  Eqs.  (2. 2. 1.1)  and  (2.2.2.21)  we  have 


E{d2l£i}  for  1 1 


(2.2.2.22) 


Substituting  this  Into  Eq.  (2.2.2.10)  and  then  sunning  we  have 


E{6}  > j 


E{£^} 

■w 


£{£^£2} 

“ETir 


(2.2.2.23) 
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In  Appendix  A2.5  It  is  proved  that 


> E{«'j^}E{A2} 


(2.2.2.24) 


Substituting  this  Into  Eq.  (2.2.2.23)  we  have  the  following  lower  bound. 


(2.2.2.25) 


Another  lower  bound  to  E{5}  la 

E{6}  > E{4} 


(2.2.2.26) 


This  follows  froa  Eq.  (2.2.2.25)  by  noting  that  E{2.^}  ^ E^{E}. 

This  concludes  the  first  part  of  our  analysis.  In  summary,  the  main 

results  up  to  this  point  are  the  derivations  of  the  upper  and  lower  bounds 

(Eqs.  (2.2.2.20),  (2.2.2.25)  and  (2.2.2.26)]  to  E{5}  that  depend  only  on 
2 

E{£}  and  E{£  }.  Our  goal,  however,  as  stated  above  Is  the  characterization 

of  E{5}  In  terms  of  X.  Therefore,  In  the  following  section,  we  will 

2 

determine  upper  and  lower  bounds  to  E{£}  and  E{£  } as  functions  of  X. 

2 

2.2.3  Characterization  of  E(t}  and  E{t  } In  terms  of  X 

In  Appendices  A2.1  and  A2.3  we  derive  expressions  for  and 

2 

where  Is  the  length  of  and  Wj  Is  the  average  nuiri>er  of 

packets  that  arrived  In  e^.  More  precisely,  it  Is  shown  that. 
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(2. 2. 3.1) 


- 1 + 2D(y/2) 


E{4j^l|y}  - 1 + 2D(y/2)  + [2D(y/2)]^  + F(y) 


(2. 2. 3. 2) 


where 


D(y)  - I 2H(y/2^) 
1-0 


(2. 2. 3. 3) 


F(y)  - I 2^[2D(y/2S(l  - 5(y/2^))  + C^(y/2S]  (2. 2. 3.4) 

1-1 


C(y)  - 1 - e’^  - ye‘^ 


(2.2. 3.5) 


Note  that  the  above  expressions  are  Independent  of  X.  The  arrival  rate  X 
enters  the  analysis  through  the  following  relationship. 


2X1. 


(2. 2. 3. 6) 


The  quantities  E{tj^j^|y}  and  i/ E{£j_j_j^|y}  are  plotted  In  Figs.  A2.2.1  and  A2.4.1. 

Equations  (2. 2. 3.1)  through  (2. 2. 3. 6)  present  the  first  and  second 

moments  of  the  length  of  one  epoch  conditioned  on  the  length  of  the  previous 

epoch  with  X being  a parameter.  In  this  section,  we  will  derive  upper  and 

2 

lower  bounds  to  E{1}  and  E{£  },  the  steady  state  first  and  second  moments  of 
Ij,  In  terms  of  X and  the  above  conditional  moments. 
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This  section  Is  divided  Into  four  subsections.  In  Subsectlon-1,  first 

we  derive  an  upper  bound  to  E{Jl}  based  on  a general  set  of  conditions,  next 

we  demonstrate  how  these  conditions  relate  to  the  arrival  rate  and  to  the 

system  parameters  and  finally  compute  the  upper  bound  to  for  various 

values  of  X.  In  Subsections  11,  111,  and  Iv,  we  develop  a lower  bound  to 

2 

E{f>}  and  upper  and  lower  bounds  to  E{)1  } respectively.  The  procedure  of  the 
analysis  In  these  subsections  Is  the  same  as  that  of  Subsection  1.  That  Is, 
first,  we  develop  a bound  based  on  a general  set  of  conditions,  next  we 
relate  these  conditions  to  the  system  parameters  and  finally  compute  these 
bounds  for  various  X's. 


1.  Upper  Bound  to  E{A} 

Here  we  will  develop  an  upper  bound  to  E{£}  that  Is  a function  of  X. 
We  begin  by  proving  the  following  theorem. 

Theorem  2. 2. 3.1;  Let  be  a positive  Integer  corresponding  to  the  state 

of  a Markov  chain  after  the  J'th  transition.  Also  assume  that  for  some 

constants  b and  a , 0 < a <1 
u — u 


< “u(ij-l)  + b for  £j  > 1 


(2. 2. 3. 7) 


lim  E{£J  < 


b - o 


r - 1 - ot, 


^ = T 


(2. 2. 3. 8) 


Proof:  Multiply  both  sides  of  Eq.  (2. 2. 3. 7)  by  p(i^j)  and  then  sun  over 
to  obtain 


(2. 2. 3. 9) 


Since  0 ^ a < 1,  Eq.  (2. 2. 3. 9)  Is  solved  recursively  to  obtain  the  following 
steady  state  solution. 


b - a 

11m  E{E  } < -Ji 

3 1 * Qt 

J-W  XI 


(2.2.3.10) 


It  can  be  shown  that  the  above  bound  Is  the  optimum  bound  over  a class  of 
bounds.  That  Is,  If  Is  nondecreasing,  and  It  Is  upperbounded  as 

shown  In  Fig.  2. 2. 3.1  then  the  tightest  upper  bound  over  all  a occurs  at  a - 1. 
Note  that  the  upper  bound  to  at  o ■ 1 Is  given  by  Eq.  (2. 2. 3. 7). 

Next  we  will  relate  the  above  theorem  to  the  system  parameters. 

By  applying  Property  3C  (Appendix  A2.2)  It  can  be  shown  that  for  any  Mq  ^ 0 


EtAj+llu)  < 1.44(y-yQ)  + ^ 1 ^0  (2.2.3.11) 


Since  y * 2XEj  and  ^ 1>  It  follows  that  mlny^  > 2X.  Therefore,  conparlng 
Eqs.  (2. 2. 3. 7)  and  (2.2.3.11)  and  noting  that  y ■ 2X1.^  tre  have  that 


a • 2.88X 
u 


(2.2.3.12) 


b - E{fj^j|y-  2X} 


4 


or  from  Eq.  (2.2.2.15) 


b - 1 + 2D(X) 


(2.2.3.13) 


Finally  substituting  Eqs.  (2.2.3.12)  and  (2.2.3.13)  Into  Eq.  (2. 2. 3. 8)  we  have 
the  following  desired  expression. 


£^a) 


2D(X)  E{alu  - 2X}  - 2.88X 
1 - 2.88X  “ 1 - 2.88X 


(2.2.3.14) 


An  Interesting  result,  that  follows  from  Eq.  (2.2.3.13)  and  the  fact  that 

a < 1,  Is  that  If  X < .347  then  I < <*>.  The  converse  of  this  statement 
u u 

follows  from  the  work  of  the  next  subsection. 

11.  Lower  Bound  to  E{E} 

Here  we  will  develop  a lower  bound  to  E{1}  that  Is  a function  of  X.  A 
by-product  to  the  work  of  this  section  Is  a lower  bound  to  the  arrival  rate 
at  which  E{1}  ■*  <».  Ue  begin  with  the  following  theorem. 


Theorem  2. 2. 3. 2;  Let  be  the  state  of  a Markov  chain  after  the  J'th 

transition,  and  let  f (E.)  be  a positive,  convex  and  nondecreasing  lower 

c 3 

bound  to  Then  if 

Ij  - fg(Ej)  (2.2.3.15) 

A * 

has  only  one  solution  at  Ej  • 1 and  fg(^j)  ^ fo*'  E^  > E then 
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(2.2.3.16) 


llni  E{«,  } > I* 

If,  on  the  other  hand,  Eq.  (2.2.3.15)  has  two  solutions,  then  there  exists  an 
Initial  state  to  the  Markov  chain  such  that 

11a  E{£  } - • (2.2.3.17) 

J-»«o  J 


Proof ; 


E{4j}  - I (2.2.3.18) 


- I (2.2.3.19) 

*3-1 

> I (2.2.3.20) 

*3-1 

> f (E{1,  ,}]  (2.2.3.21) 

— c 3-1 

Equation  (2.2.3.20)  follows  because  ^^(Ij)  (2.2.3.21) 

follows  froa  the  convexity  of  f and  Janssen's  Inequality. 

c 

Applying  Eq.  (2.2.3.21)  to  • • • • >1]^  recursively  and  using 

the  fact  that  f^  la  nondecraaslng  yields 

E{i,)  > fj{t-}  (2.2.3.22) 

J — c o 
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where 


(2.2.3.23) 


J times 


Now  assume  that  at  only  one  point.  A , and  that  f^(A^)  < for 

A 

Aj  > A . Then  f^(Aj)  Is  as  shown  in  Fig.  2. 2. 3. 2 and  since  is  convex  it 
follows  that 


11m  fJ{A-}  - f {A*}  - A* 
c 0 c 


(2.2.3.24) 


Now  If  Aj  - f(Aj)  has  two  solutions  A^^  and  A^  (A^^  > A^),  then  f (A^ ) and  A^ 

are  as  shown  In  Fig.  2.2. 3. 3.  Here,  as  can  be  seen  from  the  figure.  If 
* 

Aq  > A^  then 


11m  f^{A}  - « 

j— 


(2.2.3.25) 


Equations  (2.2.3.22),  (2.2.3.24)  and  (2.2.3.25)  together  conclude  the  proof. 

QED 

Our  next  step  will  be  to  relate  this  theorem  to  our  system,  and  to  obtain 
a lower  bound  to  E{A}  as  a function  of  X. 

It  follows  from  Property  4 (Appendix  A2.2)  that  fj^(W)  (given  below)  is 
a convex  Increasing  postlve  lower  bound  to  E{Aj^^|p}. 


f,(U)  - 1 + 


2D(v/2) 


for  V < 8 


2.88(^  - 4)  + 2D(4)  y < 8 


(2.2.3.26) 
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J*0 


«i 


Figure  2. 2. 3. 2 Illustrating  the  Conditions  for  Convergence  of  the 

Lower  Bound  to 


since  u * 2X1,  we  have  that 


f^(l)  - fj^(2Xl) 


(2.2.3.27) 


Therefore,  an  equivalent  statement  to  Theorem  2. 2. 3. 2 Is  the  following. 

* 

Let  y be  defined  by 


(2.2.3.28) 


then  If  u Is  unique  and 


fj^(y)  < ^ for  y > y* 


(2.2.3.29) 


we  have  that 


EU}>^ 


(2.2.3.30) 


otherwise 


E{1}  - » 


(2.2.3.31) 


It  follows  from  Eq.  (2.2.3.26),  that  Eq.  (2.2.3.28)  has  a unique 
solution  satisfying  Eq.  (2.2.3.29)  for  y ^ 0,  Iff 


> 1-**U  - 


(2.2.3.32) 


since  Eq.  (2.2.3.32)  holds  for  all  u ^ 0 iff  0 ^ X ^ .347,  we  conclude  that 
if  X > .347  then 

\(X)  - « (2.2.3.33) 


and  if  0 < X < .347  then  l^(\)  is  given  by  Eq.  (2.2.3.30). 
have  the  following  lower  bound  to  E{il}, 


In  sunnnary,  we 


* 

IL. 

2X 


A,(X) 


X < .347 
X > .347 


(2.2.3.34) 


ill.  Upper  Bound  to  E{ii.^} 

2 

Here  we  will  derive  an  upper  bound  to  E{i  } that  is  a function  of  X. 
He  begin  with  the  following  theorem. 

Theorem  2.2.3.3t  Let  Ej  be  the  state  of  a Markov  chain  after  the  j'th 
transition  and  assume  that 


E{Ej^l|E  } < (or  + c)^  for  E > 1 


(2.2.3.35) 


where 


0 < a <1 
— u 


(2.2.3.36) 


and 


11m  E{E  } - E{E} 


59- 


then 


2a  cE{£}  + c 
11a  E{iL/}  < 7- 


J ' - 


(2.2.3.37) 


1 - a 


Proof:  Expand  the  right  side  of  Eq.  (2.2.3.35),  multiply  by  p(Ej)  and  then 
sum  over  all  to  obtain 


< a^E{Aj}  + 2a^cE{aj}  + 


(2.2.3.38) 


2 2 

Equation  (2.2.3.38)  Is  a difference  unequallty  In  E{2.j  } where  2a^cE{£.^}  + c 
Is  the  driving  function.  Since  from  the  hypothesis  E{£.j}  ■ E{1!,}  and 
0 ^ < 1,  it  can  be  solved  recursively  to  obtain  the  following  steady 

state  solution. 


E{Z^}  < — ^ 


2a  cE{l}  + c' 


(2.2.3.39) 


1-0 


Next  we  will  relate  this  theorem  to  our  system.  In  Appendix  A2.4  we 
proved  that  E{f.^|lj}  (1  ***  1.44  y)^.  Since,  as  can  be  seen  from  Fig.  A2.4.1, 


•Sy  '^E{tJ^i|y}  - 1.44 


this  bound  can  be  generalized  as  In  the  following  expression 


E{Ej^l|y}  < 1.44(y-yQ)  + / EUj^Jy^}  for  y ^ y^  (2.2.3.40) 


where  y^  Is  an  arbitrary  parameter  yg  ^ 0 
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Now  since  M - 2X£j  and  ^ 1.  it  follows  that  min  Uq  “ 2X.  Therefore, 
comparing  Eqs.  (2.2.3.35)  and  (2.2.3.40)  we  have  that 


c - / E{aJ^Ju-2X  - 


(2.2.3.41) 


and 


a - 2.88X  (2.2.3.42) 

u 

Finally  substituting  Eq.  (2.2.3.42)  Into  Eq.  (2.2.3.39)  we  have 

E{£^}  < (X)  (2.2.3.43) 


where 


j^2.^v  5.76cX  E{E}  c^ 

“ 1 - (2.88X)^ 


(2.2.3.44) 


and  c Is  given  by  Eq.  (2.2.3.41). 

~2 

This  Is  the  desired  result.  Note  that  as  X .347,  as  well  as, 
and  Ig  approach  Infinity. 

~~2 

Next  we  will  turn  our  attention  to  the  derivation  of  t^(X).  This  la 
performed  In  the  following  subsection. 
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2 

iv.  Lower  Bound  to  £{&  } 


2 

Here  we  will  develop  a lower  bound  to  E{il  } as  a function  of  X. 
Since  the  work  here  Is  similar  to  that  of  the  previous  subsection,  the 
details  will  be  omitted. 


Theorem  2. 2. 3.4;  Let  be  as  given  In  the  statement  of  Theorem  2. 2. 3. 3 
except  for 


(2.2.3.45) 


then. 


11m  E{£^} 


E{i}  + cl 


(2.2.3.46) 


Proof:  This  proof  Is  Identical  to  the  proof  of  Theorem  2. 2. 3. 3 with  the 
Inequalities  reversed. 

QED 

Finally  In  relating  Theorem  2. 2. 3. 4 to  the  system  parameters,  first 
substitute  y * 2X£  Into  Eq.  (2.2.3.45)  to  obtain 


, o.y 

E{£^|y}  > + Cj) 


(2.2.3.47) 


Comparing  this  to  the  results  of  Appendix  A2.4  where  It  Is  shown  that 


E{£^ly}  > (1.44y  + .25)^ 


(2.2.3.48) 
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we  have  from  Eq.  (2.2.3.46)  that 


g{j^2j  ^ 1.44E{E}X  + .125 
“ 1 - (2.88X)^ 


(2.2.3.49) 


Note  here  again  that  as  X -*■  .347  then  I 


2 

I 


■*  «. 


2.2.4  E{ delay}  versus  Arrival  Rate 

In  Section  2.2.2  we  derived  upper  and  lower  bound  to  E{6}  in  terms  of 
2 

£{£}  and  E{Jl  },  and  In  Section  2.2.3  we  determined  upper  and  lower  bounds  to 
2 

E{f.}  and  E{il  } In  terms  of  X.  In  this  section  we  will  combine  these  results 

to  obtain  upper  and  lower  bound  to  E{6}  In  terms  of  X.  We  begin  with  the 

upper  bound  to  £{6}.  Eqiiatlon  (2.2.2.20)  which  Is  the  upper  bound  to  E{6} 

2 

in  terms  of  E{E}  and  E{£  } is  rewritten  below. 

E{6}  < 1.05  Ulp  + .321  (2. 2. 4.1) 

Now  from  Eqs.  (2. 2. 4. 2),  (2.2.3.34)  and  (2.2.3.44)  we  have 

E{6}  <_ — - +.321  (2. 2. 4. 2) 

1 - (2.88X)  _ (2.88X)^)Ej^(X) 


where  c and  ^^^(X)  are  given  In  Eqs.  (2.2.3.41)  and  (2.2.3.34). 

Equation  (2. 2. 4. 2)  Is  the  desired  lower  upper  bound  to  E{6}.  This 
upper  bound  was  computed  and  the  result  Is  presented  In  Fig.  2. 2.4.1. 

The  derivation  of  the  lower  bound  to  E{6}  Is  similar  to  that  of  the 
upper  bound.  From  Eqs.  (2.2.2.25)  and  (2.2.3.49)  we  have 
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1 


E{6}  > 


.72A 


1 - (2.88X)^  32E{£}(1  - (2.88X)^) 


+ i E{)1}  (2.2.4. 3) 


> — (2. 2. 4. 4) 

1 - iZ.SQXy  ^ 

Equation  (2. 2. 4. 4)  is  one  lover  bound  to  E{6},  another  follows  from 
Eq.  (2.2.2.26):  that  is 


E{6}  > )lj^(X) 


(2. 2. 4. 5) 


The  best  lower  bound  Is  the  maximum  of  the  two  bounds  given  in 
Eqs.  (2. 2. 4. 4)  and  (2. 2. 4. 5).  This  is  given  below. 


E{6} 


.72X 

^ max  

.1  - (2.88X) 


(2. 2. 4. 6) 


Equation  (2. 2.4. 6)  has  been  computed  and  it  is  plotted  in  Fig.  2.2.4. 1 along 
with  the  upper  bound. 
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Figure  2.2.4. 1 


X (packet/slot) 

Upper  end  Lower  Bounds  to  the  Average  Delay  versus 
Average  Arrival  Rate  for  the  Binary  Tree/Poisson 
Source  System 

Mote:  An  algr.  step  equals  one  round  trip  delay 
plus  two  slots 


2.3  Average  Throughput 


The  average  throughput  of  the  system  is  defined  to  be  the  fraction  of 
the  time  that  the  channel  contains  valid  data,  l.e.,  exactly  one  packet  per 
slot.  Now  If  the  average  delay  Is  finite.  It  follows  that  over  a long 
Interval  of  time,  all  the  packets  that  arrive  will  be  successfully  trans- 
mitted. Therefore,  the  average  throughput  equals  the  average  arrival  rate 
If  the  average  delay  Is  finite.  From  the  above  discussion.  It  follows  that 
the  E{ delay}  vs  X results  of  Section  2.2  may  be  Interpreted  as  E{ delay}  vs 
E{ throughput} , (see  Fig.  2. 2. 4.1),  and  the  maximum  average  throughput  Is 
.347  packets/slot. 

In  the  preceding  paragraph  we  defined  and  determined  the  maximum 
average  throughput.  It  turns  out,  however,  that  the  system  may  be  operated 
at  a throughput  up  to  .43  packcts/slot  but  only  for  a limited  time. 

In  Fig.  2. 3. 0.1  we  have  E{l|y}  versus  y.  Ve  know  from  the  results  of 
Section  2.2  that  the  average  delay  Is  finite  for  X < .347.  However,  If 
.347  < X < .430  then  ^ Intersects  E{l|y}  In  two  places,  at  y^  and  y^  which 
corresponds  to  say,  I and  I . Mow,  If  .347  < X < .430  and  1.  < I then 

8 C — * J C 

snd 


the  system  will  operate  around  However,  statistical  perturbations  will 

eventually  cause  1.  to  become  greater  than  1 , In  which  case 
3 c 


Here,  the  system  Is  expected  to  overflow;  this  Is  accompanied  by  an  In- 
crease In  delay  and  a drop  In  throughput. 
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E[l//i]  (algr.  steps) 


/I  (packets) 


Figure  2. 3.0.1  Illuatratlng  the  Attalmaent  of  the  Maxlmtim  Through- 
put of  .430  pck/alot,  for  the  Binary  Tree/Poiaaon 
Source  Syatea 
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} 

f 

2.4  System  Stability 

Heurestlcally,  an  unstable  multi-access  system  Is  one  where  the 
following  scenario  Is  possible.  Originally,  seyeral  channel  packet 
collisions  reduce  the  number  of  packets  being  successfully  tranraltted. 

This  results  In  a packet  backlog  which  further  Increases  the  number  of 
channel  collisions,  this  In  turn  Increases  the  backlog,  etc.,  until  a total 
breakdown  occurs  when  essentially  everybody  la  trying  to  get  through  with 
very  few  actually  succeeding. 

A more  precise  definition  of  stability  of  a multi-access  system  Is 
the  following.  The  system  Is  defined  to  be  k'th  order  stable  If  the  k'th 
moment  of  the  delay  Is  finite.  In  Section  2.2,  we  showed  that  If  X < .347, 
then  the  average  delay  Is  finite,  therefore,  our  system  Is  at  least  first 
order  stable.  In  the  remainder  of  this  section,  we  will  show  that  If 

X < 1/3,  then  all  the  moments  of  the  delay  are  finite.  As  has  been  pointed  I 

out  previously,  this  stability  result  Is  overly  conservative,  and  Indications 
are  that  all  moments  are  finite  for  X < .347  packets/slot.  In  Appendix  A2.8, 
a way  Is  suggested  for  determining  a larger  lower  bound  to  X max. 

From  the  results  of  Section  2.3,  we  know  that  If  E{iL  } < od  then 

1 

E{6^  < “.  Therefore,  It  Is  sufficient  to  show  that  E{E^}  < <»  for  X < 1/3 

and  for  all  k.  This  la  accomplished  by  showing  that  11m  E{e  ^}  = E{e  } 

J-*- 

exists  for  0 £ s ^ s^,  X < 1/3,  a^  > 0.  The  analysis  Is  carried  out  by 

sE 

first  obtaining  an  upper  bound  to  E{e  } under  a general  set  of  conditions 
(this  Is  accomplished  In  Theorem  2. 4. 0.1),  and  then  relating  these  conditions 
to  the  multi-access  system.  Note  that  upperboundlng  a moment  generating 


Theorem  2. 4, 0.1;  Let  be  the  state  of  a Markov  chain  after  the  j'th 
transition.  Furthermore,  let 


E{e 


8(aA. 

|Zj}<e  ^ 


+ b) 


(2. 4. 0.1) 


where  a and  b are  constant,  0 a < 1,  and  the  initial  state,  is  finite 
Then 


llm  E{e  ^}<e^“* 


(2. 4. 0.2) 


Proof:  Multiply  both  sides  of  Eq.  (2. 4. 0.1)  by  p(f-j)  and  then  sum  over 
to  obtain 


E{e“^J'^h 


saE. 

E{e  ht 


sb 


Solving  the  above  equation  recursively,  we  have 


nl.  - J-1.  . 

E{e  J}  < E{e  0}  .“bCaJ  +...+1) 


(2. 4. 0.3) 


Since  ^ “ “nd  0 ^ a < 1,  Eq.  (2.4.0. 2)  follows  by  letting  j • in 
Eq.  (2. 4. 0.3). 

QED 
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Next  we  will  show  that  the  conditional  generating  function  of  can  be 
upper  bounded  as  In  Eq.  (2. 4. 0.1)  where  a < 3X.  In  Appendix  A2.8  we  proved 
that 


E{e*®’lp}  < (1+V)e®"*^  - (A/B^  + Au/B)e"^  + (A/B^)e^®“^^^  (2. 4. 0.4) 


where 

A - — S , B - 0 < s < 8.,  s.  > 0 (2. 4. 0.5) 

2-e®  x-e®  “*  00 

and  X Is  an  arbitrary  parameter  1 < x < 3.  Since  Ue~'^  and  e ^ 1 for  u ^ 0, 
we  have  from  Eq.  (2. 4. 0.4)  that 

E{e**’ly}  < 2e®  - A/B^  - A/B  + (A/B^)e^®“^^^  (2. 4. 0.6) 


Now  for 


. 3B 
*0  38 


8-0 


X 

x-1 


Sj^  > 0 exists  such  that 


B“1  < a^s  for  0 ^ 

Q 
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(2. 4. 0.7) 


1 


d 


I 

1 


Therefore,  from  Eq.  (2.4.0. 6)  ve  have 


E{e®*'|li}  < 2e®  - A/B^  - A/B  + (A/B^)  e 


(2. 4. 0.8) 


8(aQy  + b) 

1 e for  b > 0 < s < 82,  S2  > 0 (2. 4. 0.9) 


Equation  (2. 4. 0.9)  follows  from  Eq.  (2.4.0. 8,  because  at  s - 0,  both  equal  1, 

g 

both  are  continuous  and  of  Eq.  (2. 4. 0.9)  is  greater  than  that  of 
Eq.  (2. 4. 0.8). 

Finally  by  substituting  i • and  y - 2X1^  into  Eq.  (2. 4. 0.9)  we 

obtain 


si, s(2a.XE  b) 
E{e  ° J 


(2.4.0.10) 


si. 


From  Theorem  2. 4. 0.1,  we  know  that  for  convergence  of  lim  E{e  we  need 

J— 


X < 


2a, 


and  from  Eq.  (2. 4. 0.7) 


X 


< 


x-1 

2x 


(2.4.0.11) 


-71- 


since  I < X < 3,  we  have 


X < sup  [^] 
l<x<3 


1 

3 


(2. A. 0.12) 


This  concludes  this  section. 
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APPENDIX  A2 


PROPERTIES  OF  THE  BINARY  TREE/POISSON  SOURCE  SYSTEM 
A2.1  Derivation  of  E{fc|u} 

If  ve  let  the  number  of  packets  that  are  to  be  processed  In  an  epoch  be  a 
Poisson  random  variable  with  mean  y and  let  t be  the  number  of  algorithmic 
steps  required  to  process  these  packets  by  the  binary  tree  algorithm,  then 
In  this  appendix  we  will  show  that 

E{ily}  - 1 + 2D(y/2)  (A2.1.1) 


where 


D(y)  S I 2H(y/2^) 
1-0 


(A2.1.2) 


and 


C(y)  i 1 - e''^  - ye"'*  (A2.1.3) 

Before  proving  this  result  several  definitions  will  be  given.  (Fig.  A. 2. 1.1 
might  be  helpful  In  visualizing  these  definitions) 

n^j  - the  J'th  node  of  level  1.  There  are  2^  nodes  at  level  1.  (A2.1.4) 

T^^  - subtree  whose  root  node  Is  n^^  (A2.1.5) 

- a random  variable  that  equals  one  if  n^^  Is  visited  by  the 
algorithm  and  zero  otherwise.  Note  that  n^^  Is  always  visited 
by  the  algorithm,  therefore,  x^  • 1.  (A2.1.6) 
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As  has  been  pointed  out  an  algorithmic  step  Is  equivalent  to  visiting 
one  node.  Therefore,  the  total  number  of  algorithmic  steps  equals  the 
total  number  of  nodes  visited.  Motivated  by  this  observation,  we  prove  the 
following  lemma  which  Is  useful  In  determining  the  probability  that  a 
particular  node  Is  visited. 


Lemma  A2.1.1.  A node  n^  1^0  Is  transversed  In  the  binary  tree  algorithm  if 
and  only  If  there  are  at  least  two  active  sources  In 


Proof:  Let  [n  m - 0,  1,...,  1]  be  the  set  of  all  the  nodes  that  lie  on 
nir 

the  path  from  n^^  to  n^^j.  Since  Is  Included  in  : m - 0,  1,...,  1], 

It  follows  that  If  contains  at  least  two  active  sources  so  does 

[T  _ : m < 1].  Mow  If  the  protocol  Is  on  node  [n  : m < 1]  and  asks  for  the 
mr  — mr 

nioaber  of  active  sources  In  T the  answer  will  be  greater  than  one,  and 

iBTi,r 

It  will  move  to  node  n . Since  this  holds  for  at  least  m - 0,  1,...,  1-1 

iik^Xtr 

we  conclude  that  node  n^^  will  be  transversed  by  the  algorithm.  If,  on  the 

other  hand,  there  Is  at  most  one  source  In  and  the  algorithm  got  at  least 

as  far  as  n.  ^ and  asked  for  the  number  of  active  sources  In  T. .,  the  answer 
i-i,r  ij 

will  be  0 or  1 and  It  will  not  continue  to  n^^ . 

QED 


Theorem  A2.1.1;  Let  there  be  v active  sources  (where  v Is  a Poisson  random 
variable  with  mean  y)  and  let  £ be  the  nuniber  of  nodes  transversed  by  the 
binary  tree  algorltlai  In  resolving  the  conflict,  l.e. , the  process  of 
separating  all  the  sources  Into  subsets  such  that  each  subset  contains  at 
most  one  active  source.  Then, 
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i 

I 


E{«,|y}  - 1 + 2D(y/2) 


(A2.1.7) 


where 


D(y)  - I ihivnh 

i-0 


(A2.1.8) 


and 


5(y)  - 1 - e~^  -ye"'^  (A2.1.9) 

Proof:  From  the  definitions  In  Eqs.  (A2.1.4)  to  (A2.1.6)  we  have  that 


4 - Z z (A2.1.10) 

1-0  J-0 


From  Leomia  A2.1.1  we  have  that 


p(x^j-l|y) 


p(At  least  2 active  sources  In  T^^jy)  for 
1 If  1-0  (A2.1.11) 


But  since  the  number  of  active  sources  Is  Poisson  distributed,  and  the  sources 
are  Independent  of  each  other.  It  follows  that  the  number  of  active  sources 
In  Is  also  Poisson  distributed  with  parameter  y/2^.  Therefore,  from 
Eq.  (A2.1.11)  follows  that 


(A2.1.12) 


P(*lj*l|w)  - 1 - - ii/2^  i ,«  0 

i C(y/2S 

and  froa  Eq.  (A2.1.10)  we  have  that 

« 2^-1 

eU|w}  • I I p(x^j-i|v) 

i-0  j-0 

and  from  Eq.  (A2.1.11)  and  (A2.1.13) 


(A2.1.13) 


(A2.1.14) 


E{2|y}  - 1 + I 2h(v/2^) 
1-1 


(A2.1.15) 

QED 
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A2.2  Properties  of  E{i!,|u} 

Here  we  will  derive  several  properties  of  E{jl|y}.  Some  of  these 
properties  are  trivial,  whereas,  others  are  more  Involved.  However,  all 
of  them  are  presented  here  for  completeness.  The  quantity  E{l|y}  which 
was  derl/ed  In  Appendix  A2.1,  Is  rewritten  below  and  plotted  In  Fig.  A2.2.1. 


E{Jl|y}  - 1 + 2D(y/2)  (A2.2.1) 

00 

D(y)  - I iHivUh  (A2.2.2) 

1-0 

iiU)  - 1 - e"^  -ye‘^  (A2.2.3) 


Since  E{£|y}  depends  trivially  on  D(y)  and  since  D(y)  Is  a basic  quantity 
that  appears  In  other  expressions,  we  will  develop  the  properties  of  D(y). 

The  corresponding  properties  of  E{£|y}  follow  from  those  of  D(y)  In  a 
straight  forward  manner.  The  properties  of  D(y)  that  are  derived  here  are 
listed  on  page  84.  In  order  to  motivate  these  properties  we  calculated 
D(y),  D'(y)  and  D"(y)  for  ye [0,16].  The  results  are  shown  In  Figs.  A2.2.2, 
A2.2.3  and  A2.2.4. 

Before  beginning  the  derivations  we  will  comment  on  one  of  the  techniques 
used  In  this  appendix  as  well  as  several  other  places  In  the  thesis.  At  times 
we  need  to  show  that  a function  Is  positive  over  some  finite  range  of  Its 
argument.  This  la  accomplished  here  simply  by  calculating  the  function  at 
several  points.  An  exaisple  of  the  use  of  this  technique  can  be  found  In 
Property  4 where  It  Is  shown  that  — r > 0 for  0 < y < 4.  We  are  awere 
of  the  fact  that  calculating  the  function  at  a finite  number  of  points  does 
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iji&ist. 


5 


0 12  3 4 

ft  (packets) 


Figure  A2.2.1  B{i|y}  versus  u for  the  Blnery  Tree/Poieeon  Source 
) Syetea 
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Flgur*  A2.2.2 


12  3 4 

/I  (packets) 


D(v)  for  ClM  Binary  Traa/Polsson  Sourca  SysCaa 
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D'(/i) 
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0"(m) 


not  constitute  a proof  of  Its  properties  over  a range.  However,  given  the 
smoothness  of  the  functions  Involved  (Properties  2A,  2B  and  2C),  It  Is  felt 
that  this  technique  Is  justified.  A more  rigorous  approach  would  be  to  find 
a lower  bound  that  Is  valid  over  a small  finite  region  and  then  show  by 
computation  that  this  lower  bound  Is  positive  over  the  desired  range.  For 
example,  by  taking  the  second  derivative  of  D(vi)  we  have 


It  follows  easily  that 


D"(y)>  I -f V-e  ° 

2^  2*- 

1-0  ^ ^ 


for  Pq  1 U 1 Wq 


So  to  show  rigorously  that  D* ' (y)  ^ 0 for  0 ^ y ^ 4 one  need  simply  show  by 
computation  that  the  above  lower  bound  Is  positive  at 


Pg  ■ KAy  for  some  Ay  and  K - 0,1,2, . . . , |^| 

As  was  pointed  out  above.  It  la  felt  that  this  mors  rigorous  technique  Is 
not  necessary  for  our  purpose,  therefore,  when  needed  to  show  properties, 
such  as  positiveness,  of  a single  variable  function  over  some  small  range 
It  will  be  accomplished  by  computing  the  function. 

The  properties  that  will  be  proved  are  listed  below. 
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i 


00 

Property  lA  C(u)  - Z 


k+1 


k^l 


IB 


2 3 2 

^ ^ £ 5(y)  ± ^ for  0 £ y £ 1. 2 


1C 


k-l 


2'‘-l 


11) 


y^  - |y^  < D(y)  < 0<m<i.2 


2A  D(2y)  - 2D(y)  + C(2y) 


2B  D'(2y)  - D'(y)  + e’(2y) 


2C  D'’(2y)  - I D"(u)  C"(2y) 


2D  D(y)  - I C(y/2^)  + I D(y/2S 

1-0  1-1 


3A  Bln  D'(y)  - aln  D'(y)  = B.(y-) 


y>yji 


2y^>y^^ 


Bj(4)  - 1.440 
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3B  D'  (v)  2^  0 for  y ^ 0 

3C  max  D' (y)  i 8 - 1.443 

y>0  “ 

3D  D(y)  ^ S^y 

D(y)  0 < y < 4 

4 Ut  f - “ 

D'(4)(y-4)  + D(4)  y>4 

where  D'(4)  * 1.44  and  D(4)  ■ 4.77,  then  £ la  a convex 

Increasing,  positive  lower  bound  to  D(y). 

For  the  main  text,  we  will  use  values  of  8£(4)  and  8^  that  are  rounded 

off  to  three  significant  figures.  Therefore  we  will  assume  that  6j^(4)  - 8^  ■ 

1.44.  The  error  caused  by  this  approximation  Is  Insignificant. 


Here  we  will  show  that 


k-1 


(A2.2.4) 


2 3 2 

C(y)  for  0 < y < 1.2  (A2.2.5) 

Proof:  Substitute  the  Taylor  series  expansion  of  e ^ Into  Eq.  (A2.2.3)  to 
obtain 

t 

J 
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(A2.2.6) 


5(W) 


«» 

1-  I 

k-0 


ki 


k-0 


Equation  (A2.2.6)  Is  Chen  rearranged  to  obtain  (A2.2.4).  Equation  (A2.2.S) 
follows  froB  (A2.2.4)  becuase  for  u ^ 1.2  Che  sum  In  Eq.  (A2.2.4)  Is  Che 
SUB  of  an  alternating  sequence  whose  magnitude  Is  decreasing. 


Properties  1C  and  ID.  Here  we  will  prove  that 


D(y) 


m 

I 

k-l 


k 

(kfl)! 


(-U) 


k+1  2' 


2*^-1 


(A2.2.7) 


and 


^ < D(y)  for  0 < y £ 1.2  (A2.2.8) 

Proof:  To  obtain  Eq.  (A2.2.7),  first  substitute  Eq.  (A2.2.4)  Into  Eq.  (A2.2.2). 

D(y)  - I I 2^  (A2.2.9) 

1-0  k-l 

Equation  (A2.2.7)  follows  from  (A2.2.9)  by  first  Interchanging  the  order  of 
suBaatlon  and  then  susBlng  over  Che  Index  1. 

Equation  (A2.2.8)  follows  froB  Eq.  (A2.2.7)  because  the  sub  la 
Eq.  (A2.2.7)  Is  that  of  a decreasing  alternating  sequence. 

QED 
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Propertlea  2A.  2B.  and  2C.  Here  it  will  be  proved  that 


D(2y)  - 2D(y)  + C(2y) 

(A2.2.10) 

D'(2y)  - D’(y)  + 5'(2y) 

(A2.2.11) 

D"(2y)  - 1 D"(u)  + C'’(2v) 

(A2.2.12) 

Proof:  Fron  Eq.  (A2.2.2)  followa  that 

00 

D(2y)  ~ J,  2^  ^(u/2^~h 

1-0 

(A2.2.13) 

00 

• 2 1 2^”^  5(u/2^'^) 

1-0 

- 2(|  C(2y)  + D(y)] 

therefore 

D(2y)  - 2D(u)  ♦ C(2y) 

(A2.2.14) 

Equatlona  (A2.2.11)  and  (A2.2.12)  follow  from  Eq.  (A2.2.14)  by  differen- 
tiating with  reapect  to  y. 
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Propert’ 

1 2D.  Here,  we  will  prove  that 

00  00 

D(y)  - 1 C(y/2^)  + 1 D(y/2^) 

(A2.2.15) 

Proof : 

1-0  1-1 

Fron  Eq.  (A2.2.2)  we  have  Chat 

00  00 

D(y)  - 1 C(y/2^)  + 1 (2^-i)5(y/2h 

(A2.2.16) 

Now, 

1-0  1-0 

1-1 

2^-1  - I 2^^ 

80 

1^0 

. - 

1 (2^-l)5(y/2^)  • 1 1 

(A2.2.17) 

1-0  1-1  k-0 

Interchange  order  of  eoBMClon 

ml  1 

k-O  1-fcfl 

• 1 1 C(y/2^'^) 

k-O  1-1 
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IncerchAngc  order  of  suBBUtlon  once  again  to  obtain 


- I I (A2.2.18) 

1-1  k-0 


But 


I 2*^C(u/2^'^^  - D(M/2h  (A2.2.19) 

k-0 


Equation  (A2.2.1S)  follows  by  substituting  Eq.  (A2.2.19)  Into  Eq.  (A2.2.18) 
and  Eq*  (A2*2*18)  Into  Eq*  (A2*2*16)* 


Property  3A*  Let 


8,(y,)  - nln  D*(y)  (A2.2*20) 

Then  Property  3A  states  that 


D* (w)  > 8£(Mt)  tor  y > 


(A2*2.21) 


and  Property  3B  Is 


O' (y)  > 0 for  y > 0 


(A2.2.22) 
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Proof:  From  Eqs.  (A2.2.3)  and  (A2.2.11)  we  have 


D'(2u)  - D'(w)  + 2ye'^ 


(A2.2.23) 


since 


2ve“^^  > 0 


It  follows  that 


D'(2v)  > D’(w) 


therefore 


Min  D'(p)  £ nln  D'(y) 

Vj<V<2Vj  2yj^W<4Uj^ 


(A2.2.24) 


Equation  (A2.2.21)  follows  by  Induction  from  Eq.  (A2.2.24). 


QED 


The  quantity  6£(4)  was  calculatad  and  It  Is  given  by 

B£(4)  - 1.440  (A2.2.25) 
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Property  3B.  Prove  that 


D'(u)  > 0 for  w > 0 


(A2.2.26) 


Proof:  Differentiating  Eq.  (A2.2.2)  we  have 


D'(li)  - I (u/2^)e"‘'^^^ 
1-0 


(A2.2.27) 


which  la  positive  for  u > 0. 


Property  3C.  Let 


then 


OB  ^ 

0'(U)  - I (y/2S«"^'^ 


1-0 


D’,(c)  - I \ 

1—  2*^ 


U > 0 


(A2.2.28) 


sup  D'(v)  - uea  D'^(c) 
U>0  l<c<2 


(A2.2.29) 


(A2.2.30) 


-91- 


First  we  will  prove  Eq.  (A2.2.29)  and  then  we  will  calculate  6 . 

u 

Proof:  Substituting  Into  Eq.  (A2.2.28} 

U - c2^  for  l<c<2  and  k-0,1,2,... 
and  then  rearranging  we  have 

D'(c2*^)  - I (A2.2.31) 

1— k 2^ 

Since  the  sinaand  In  Eq.  (A2.2.31)  Is  positive  we  have 

llB  D'(c2S  - D'  (c) 
kr*« 

> D'(c2^)  for  llc<2,  J-0,1,2,...  (A2.2.32) 

QED 

The  Moclausi  of  D^(c)  over  lfc^2  was  calculated  and  It  equals  1.443. 

Therefore 

®u  ■ (A2.2.33) 

Property  3D  D(u)  < 6^  W (A2.2.34) 

This  follows  froa  Property  X and  tha  fact  that  D(0)  > 0. 
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Property  4.  A postlve, convex  non-decreasing  lower  bound  to  D(y)  Is  f(u) 
where 


f(y)  - 


D(vi)  for  0 £ y £ 4 

D'(4)[y-4]  + D(4)  for  y > 4 


(A2.2.3S) 


Proof:  From  Property  3B  we  have  that  D'(v)  ^0  and  since  D(0)  • 0 it 
follows  that  f(y)  is  positive  and  non-decreasing.  That  D(y)  is  convex 
for  0 ^ U 4 follows  from  Fig.  A2.2.4  and  from  a more  detailed  computer 
printout  of  — 7.  To  show  that  f is  convex  for  y > 0,  note  that 

~ - 0 for  y > 4 (A2.2.36) 

3y^ 

and  that  f(y)  is  continuous  and  has  a continuous  first  derivative  at  y ■ 4. 
Therefore  f(y)  is  convex. 

Next  we  will  show  that 

f (y)  < D(y)  for  y ^ 0 (A2.2.37) 

Since  f(y)  ■ D(y)  for  y ^ 4,  Eq.  (A2.2.37)  would  be  true  if  f'(y)  £D'(y) 
for  y ^ 4.  But 

f*(y)  - D’(4)  for  y > 4 (A2.2.38) 
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so  we  need  to  show  that 


D'(4)  < D'(p)  for  p > 4 


(A2.2.39) 


From  Property  3a 


min  D'  (p)  ^ D'  (p)  for  p ^ 4 
4<p<8 


(A2.2.40) 


D'(p)  was  computed  for  4 ^ p ^ 8 (see  Fig.  A2.2.2).  The  minimum 
did  occur  at  4 where 


D’(4)  - 1.440 


QED 
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( 


A2.3  Derivation  of  E{£^lu} 

•y 

In  this  section  ve  will  derive  an  expression  for  Ed^jy}.  This  Is 
accomplished  In  Theorem  A2.3.1.  However,  before  this  theorem  Is  proved 
It  will  be  necessary  to  prove  the  following  Lemma. 

Lemma  A2.3.1.  Let  be  a random  variable  as  defined  in  Eq.  (A2.1.6) 
with  the  probability  density  given  In  Eq.  (A2.1.11).  Then  for  l^i^l 

C(u/2^)  If  n^^j  and  n^  lie  on  the  same  path 
e(v/2^)C(y/2”)  otherwise  (A. 2. 3.1) 

Note,  two  nodes  (n^^  and  n^^  : 1 ^ m)  lie  on  the  same  path  If  n^  lies 
on  the  line  connecting  n^^  and  n^^. 

Proof:  If  n^j  and  n^  lie  on  the  same  path,  then 

P(x^-l|x^j-l)  - 1 {A2.3.2) 


and 


p(x^j-l|y)  - C(U/2S 


(A2.3.3) 


therefore 


p(x^j-l  and  x^-l|y)  - 5(y/2^) 
Next  note  that 

1 Iff  X..-1  and  X .-1 
IJ  mk 

0 othervlae 

Therefore « from  Eqs.  (A2.3.4)  and  (A2.3.5) 


(A2.3.4) 


(A2.3.5) 


“mk  ^ “1  * ?;(u/2^) 


(A2.3.6) 


If,  on  the  other  hand,  n^^  and  n^^^  do  not  lie  on  the  same  path,  then  x^^  and 
x^  are  Independent  because  the  sources  of  and  are  Independent. 


Therefore 


E[Xlj*mkl'^*“lj  “mk  °“  " 5(u/2^)5(u/2")  (A2.3.7) 

QES 

A 

Now  we  are  ready  to  proceed  with  Theorem  A2.3.1.  where  E[1  |y]  is 
derived. 

f 
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4. 


Theorca  A2.3.1.  Let  V be  a Poisson  random  variable  with  mean  y corresponding 
to  the  number  of  contending  sources.  Also  let  I (given  In  Eq.  (A2.2.10)  and 
rewritten  below)  be  the  number  of  nodes  transversed  by  the  binary  tree 
algorithm. 

« 2^-1 

t - 1 + Z I (A2.3.8) 

1-1  j-0 

Where  x^  Is  defined  be  Eqs.  (A2.1.6)  and  (A2.1.11),  then 

% 

E[t^|y]  - 1 + 2D(y/2)  + [2D(y/2)]^  + F(y)  (A2.3.9) 


where 


D(y)  - I 2^C(y/2^) 

1-0 


(A2.3.10) 


and 


P(y) 


2 I 2S(y/2^)(l  - C(y/2Sl  + Z 2^  C^(y/2S 

1-1  1-1 


(A2.3.11) 


Proof:  Square  both  sides  of  Eq.  (A2.3.8)  and  then  take  the  expectation  to 
obtain 


2^- 

E(£^|y]  - 1 + 2E  I I 
a-1  j-0 


1 

8 

‘ljl“ 

+ E 

n.  i iw 

• 

L\i-i  J-0  / J 

(A2.3.12) 
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From  Theorem  A2.1.1  and  Eq.  (A2.3.10),  the  above  expression  can  be  simplified 
to 


1 + 4D(v/2)  + E 


2^-1 


I I 

^-1  j-0 


(A2.3.13) 


Next  note  that 


« 2^-^  \2  » 2^'^ 
I I - I I 


.1-1  j-0 


i-1  j-0 


2’"-l 


2 I I 

m>l  n-n(l,m) 


X..X  - X.. 

Ij  mn  ij 


(A2.3.14) 


where 


n(i,m) 


j for  m-1 
0 otherwise 


Now  the  sums  over  m and  n will  be  decomposed  into  three  sums  and 

^ij  * ^ij  over  all  (m,n)  e >2^^  is  the  sum  over  all 

(m,n)  t and  m > i,  and  is  the  sum  over  all  (m,n)  i such  that 
m-1,  n>j.  Figure  A2.3.1  might  be  helpful  in  visualizing  the  preceding. 

In  terms  of  the  above  definitions  we  have 

^Ij  " ^ *lj*mn  (A2.3.15) 

(m,n)eT^j 

hi"  I I *ij*»n  (A2.3.16) 

kex  ■ 

I 

i 
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Where  kEK  If  T.  ..  . fl  T.,»0;  note  that  there  are  2^^^-’’  such  T. . . 

IJ  1+1  ,k 


2^-1 

”ij  “ ^ *lj*ln 

n-j+1 


(A2.3.17) 


Finally  from  Eqs.  (A2.3.14)  through  (A. 2. 3. 17)  we  have 


/ » 2 -1  \ 2 00  2^-1 

( I I *lj  1 “ ^ ^ (A2.3.18) 

\i-l  j-0  / 1-1  J-0 

Next  we  will  determine  E[Y^j|ii],  E[Z^j|y],  E(W^j|w]  and  E[Xy|vi]. 

First  note  that  from  Lemma  A2.3.1  we  have 

E[Xy^|yl  - C(v/2^)  (A2.3.19) 

Now  from  Leomia  A2.3.1  and  Eq.  (A2.3.15) 

EtYyly]  - I (2“"^)C(y/2“) 
m-1 

and  from  Eq.  (A2.3.10) 

E(Y^jlyl  - D(y/2^)  (A2.3.20) 
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Next  apply  Lenma  A2.3.1  to  Eq.  (A2.3.16)  and  use  symnetxy  to  obtain. 


E[Z^j|y]  - (2^'^^-2)  I 2“-^-H(w/2S5(y/2") 

OPl+l 


Simplifying 


E[Zy(y]  - 2(2^-l)C(y/2^)D(y/2^‘^b 


Similarly 


EfW^jIyJ  - (2^-l-J)C^(y/2S 

Next  take  expectations  of  both  sides  of  Eq.  (A2.3.18),  substitute 
Eqs.  (A2.3.19),  (A2.3.20),  (A2.3.22)  and  (A2.3.23)  into  (A2.3.18)  sum 
and  then  substitute  this  into  Eq.  (A2.3.13)  to  obtain, 

- 00 

E(Z^|yl  - l+2D(y/2)  + I (2(2^)D(y/2^)+4(2^)(2^-l)5(y/2^)D(y/2^‘^^) 

i-1 

+2^2^-l)e^(y/2*)) 

Now  from  Property  2A  we  have 

D(y/2^‘^^)  - I D(u/2^)  - | C(y/2^) 


(A2.3.21) 


(A2.3.22) 


(A2.3.23) 


over  j 


(A2.3.24) 


(A2.3.25) 
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Substituting  this  into  Eq.  (A2.3.24)  and  then  rearranging  we  have. 


E[a^|y]  - l+2D(y/2)  + I 2^[2D(y/2 V2(2^-l)D(y/2^)C(y/2S-(2^-l)C^(y/2^)] 

(A2.3.26) 

Next  note  that 

00  2 

(2D(y/2))^  - [ I 2^  5(y/2^)]  iA2.3.27) 

1-1 


A2.4  Upper  and  Lower  Bounds  to  Effc'^lul 
In  Appendix  A2.3  we  showed  that 

Etil^ly]  - 1 + 2D(y/2)  + [2D(y/2)]^  + F(y) 


where 

CO 

D(y)  - I 2^  C(y/2S 
1-0 


00  OP 

F(y)  • Z I D(y/2S[l-C(y/2^)]  + I 2V(y/2S 

i-1  i-0 


The  quantity  /Ell^ly]  was  calculated  and  It  Is  plotted  In  Fig.  A2.4. 
this  Appendix  we  will  derive  upper  and  lower  bounds  to  E[t^|y].  More 
specifically  we  will  show  that 

E(£^|yl  1 (1  + 1.44y)^ 

Et4^|yl  > (.25  + 1.44y)^ 


We  begin  with  the  proof  of  Eq.  (A2.4.4).  Two  other  quantities  will  be 
needed  In  the  analysis  that  follows.  These  are 


(A2.4.1) 

(A2.4.2) 

(A2.4.3) 

In 

(A2.4.4) 

(A2.4.5) 
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(A2,4.6) 


F’  (16)  - max  F’  (y)  - 4.369 
y>16 

D(8)  - 10.54  (A2.4.7) 

The  compuCaClons  of  F| _ (16)  Is  somewhat  Involved  and  It  Is  performed 

at  the  end  of  this  Appendix.  The  value  for  F'  (16)  given  In  Eq.  (A2.4.6) 

tnx 

is  calculated  In  Eqs.  (A2.4.11)  through  (A2.4.22).  The  computation  of 
D(8)  is  straightforward  and  no  special  explanation  for  Its  computation 
will  be  given. 

From  Fig.  A2.4.1  and  from  a more  detailed  computer  printout.  It  can  be 
seen  that  Eq.  (A2.4.4)  Is  true  for  ()^y^l6.  The  validity  of  this  equation 
for  y>16  will  follow  If  we  prove  that 

1^  E[Jl^|y]  (1  + 1.44y)^  for  y > 16  (A2.4.8) 

We  begin  by  differentiating  Eq.  (A2.4.1) 

1^  ElE^ly]  - D*(y/2)  + 4D(y/2)D’(y/2)  + P’(y)  (A2.4.9) 

Now  by  applying  Property  3C  we  can  upper  bound  the  right  side  of 
Eq.  (A2.4.9)  as  follows: 
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E(£^ly]  1 By  + 46yt^(y-16)+D(8)]  + F^(16)  for  y > 16 


or 

E(A^|y]  < 4.147y  + .164  for  y > 16  (A2.4.10) 

Equation  (A2.4.8)  follows  by  carrying  out  the  Indicated  differentiation 
on  Its  right  side  and  then  comparing  the  result  to  Eq.  (A2.4.10). 

The  above  argument  was  baaed  on  Eq.  (A2.4.6).  Next,  we  will  prove 
this  equation.  Now  from  Eq.  CA2.4.3)  we  have 


F(y)  - I 2^  f(y/2^) 
1-1 


(A2.4.11) 


where 


f(y)  = 2D(y)(i-C(y))  + TCy) 

Differentiating  Eq.  (A2.4.11)  we  have 

OD 

F'(y)  - I f(y/2^) 

1-1 

Next  any  y 2^  16  can  be  expressed  as 


y - x2  where  8<x<16  and  k - 1,2,3,, 


(A2.4.12) 


(A2.4.13)  t 


(A2.4.14) 
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Substituting  Eq.  (A2.4.14)  into  Eq.  (A2<4.13)  and  then  rearranging  we  have 


k 

F'(x2*^)  - F'(x)  + I f'(2^x) 

1-1 

Next  we  will  show  that 
k 

^ f'(2^x)  ^ 10  ^ for  1^1  and  8^x<16 

1-1 


(A2.4.15) 


(A2.4.16) 


Differentiating  Eq.  (A2.4.12)  and  then  substituting  D(vi)  - 2D(y/2)+C(U) 
(Property  2A)  we  have 

f'(y)  - 2D'(y)(l-^(y))-4D(u/2)5'(w) 

- 2e"^lD’ (y)(l-hj)  - 2D(y/2)y)  (A2.4.17) 

A slnple  upper  to  f'(y)  Is 


f’(y)  < 2e"‘*D'(y)(l+u) 


(A2.4.18) 


(A2.4.18)  can  bet  further t upper  bounded  by  applylt^^'  Property  3C. 
Doing  this  we  have 


f’(y)  < 2B^(14v)e“'‘ 


(A2.4.19) 
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therefore 


I f'(2^x)  < I 2ey(l+2^x)e"^  * < 10“^  for  xe[8,16]  (A2.4.20) 

1-1  1-1 


Finally  from  Eq.  (A2.4.15)  we  have  that 

max  F'(y)_<  max  F'(y)  + 10"^  (A2.4.21) 

VI^16  8£y£l6 

The  right  side  of  Eq.  CA2.4.21)  was  calculated  and  Is  given  by 

- ^-369  (A2.4.22) 

This  concludes  the  proof  of  Eq.  (A2.4.4). 

The  proof  of  Eq.  (A2.4.5)  Is  similar  to  that  of  Eq.  (A2.4.4)  and  will  not 
be  given  here.  That  Eq.  (A2.4.S)  Is  true  should  be  obvious  from  Fig.  A2.4.1. 


In  this  Appendix  we  will  prove  that  ^ E[Aj^]E[ll2]  where  and 

^2  are  the  lengths  of  two  consecutive  epochs.  This  will  be  accomplished  In 
two  steps;  first.  In  Lemma  A5.1  It  Is  shown  that  E[£r2l^]^l  a non-decreasing 
function  of  and  then  In  Theorem  A5.1,  the  main  result  Is  proved. 

Le— a A5.1.  Let 

£[^21^1  “ 1 2D(y/2)  (See  Eq.  (A2.1.7))  (A2.5.1) 


where 


y - 2Xf,^  (See.  Eq.  (2. 2. 1.1)) 
and  X > 0 Is  the  packet  arrival  rate.  Then 


(A2.5.2) 


(A2.5.3) 


Proof:  From  Eqs.  (A2.5.1)  and  (A2.5.2)  we  have 


Fro*  Proporty  3B,  D'(ij)  ^ 0 and  since  X ^ 0 it  follows  that  i ®* 


Now  wa  are  ready  to  prove  the  main  result. 


Theorem  A2.5.1.  Let  £[£.2!^]^]  non-decreasing  in  Then,  assuming  all 
expectations  to  exist. 


1 E(£2]E(£2] 


(A2.5.5) 


Proof:  An  equivalent  statement  to  Eq.  (A2.5.5)  is 


E[(i2-A^)(42-^2^^  - ° 


(A2.5.6) 


Et(4^-T^)a2-1^2>  “ I (V^i)[E[42Ui)-^2^P^^1^ 


l_  ^ (E2-Il)tE(E2|E2)-lf2^P^^l^ 


4 <1 
1—1 


SL  >X 
1 1 


(A2.5.7) 


i I <Ti-ti)lT2-E<l2|ti)Jp(li)  + J (tl-Tl)IE(t2|T,)-T2)p(t,) 

*1-^1  ‘l>Tl 


(A2.S.8) 


(A2.5.9) 
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Equation  (A2.5.8)  follows  from  (A2.5.7)  because  la  non-decreasing. 

Carrying  out  the  sunmatlon  In  Eq.  (A2.5.8)  results  In  Eq.  (A2.5.9), 


This  concludes  Appendix  A2.5. 


QQ) 
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A2.6  Derivation  of  E[dlul 

Assume  that  there  are  V ^ 1 packets  to  be  processed  by  the  binary  tree 
algorithm  where  v Is  a Poisson  random  variable  with  mean  u,  and  let  d 
correspond  to  the  number  of  nodes  visited  by  the  algorithm  before  a randomly 
selected  packet  from  the  set  of  the  v active  ones  Is  successfully  transmitted. 
In  this  Appendix  we  will  determine  an  expression  for  E[d|y].  Note  that 
E[^|yl>  as  derived  In  Appendix  A2.1,  Is  the  expected  number  of  nodes  visited 
given  y,  whereas  E[d|y]  equals  the  expected  number  of  nodes  visited  before  a 
randomly  selected  packet  Is  transmitted.  The  relationship  between  E[Jl|y]  and 
E[d|y]  Is  examined  In  Appendix  A2.7. 

The  binary  tree  Is  shown  In  Fig.  A2.6.1.  As  has  been  pointed  out 
previously,  the  sources  correspond  to  the  leaves  of  this  tree.  Another 
representation  of  the  sources  which  will  be  more  convenient  for  the  work  In 
this  Appendix  Is  to  represent  each  source  by  the  binary  number  S where, 

S ■ SqS^S2*«**  (A2.6.1) 

The  equivalence  between  the  two  representations  Is  established  through  the 
following  convention. 

If  >2^*0  take  the  upper  branch  emanating  from  the  node  of 
level-1,  and  If  e^*l  take  the  lower  branch.  This  Is 

performed  consecutively  beginning  with  1-0,1 etc. 

In  Fig.  A2.6.1,  for  example,  the  circled  nodes  correspond  to  S-0110... 


-112- 


Below  are  some  definitions  that  we  will  need. 

S ■ SqS^S2...  the  blt-by-blt  complement  of  S.  (A2.6.2) 

E(d^|u]  > the  expected  delay  source  S undergoes,  given  u*  (A2.6.3) 
E[d|u]  - the  expected  delay,  l.e. , E[d^|y]  averaged  over  S.  (A2.6.4) 

* the  subtree  of  level-1  whose  root  node  Is  one  branch 

away  from  S,  (see  Fig.  A2.6.1)  (A2.6.5) 

m 

D(y/2S  - I 2^  C(y/2^'^^)  (A2.6.6) 

j-0 

This  Is  the  expected  number  of  nodes  that  must  be 
visited  by  the  algorithm  in  order  to  process  all  the 
active  sources  of  T^.  (See  Appendix  A2.1) 

■ expected  number  of  nodes  lying  on  path-S  that  are  used 
to  process  source  S given  y.  (A2.6.7) 

- the  expected  number  of  nodes  lying  above  S that  are 

visited  before  source  S transmits,  given  y.  (A2.6.8) 
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(A2.6.9) 


e(y/2^) 


. €(y/2^) 


This  Is  the  probability  that  there  are  at  least  two 
active  sources  In  given  that  there  Is  at  least 
one  active  source. 


In  terms  of  the  above  definitions  we  have 


E[d  ly]  - X + Y (A2.6.10) 

8 3 8 

Next  we  will  calculate  E[d^|y]  and  E[d^|yl  and  show  that  E[d^|y]  + E[d^|y]  is 
Independent  of  S.  This  result  will  lead  us  to  conclude  that 


Eldly]  - i [E[d^ly]  + Etdyly]] 


(A2.6.11) 


First  we  will  calculate  X . A node  of  level  1 will  be  visited  by  the 

s 

algorithm  if  there  are  at  least  two  active  sources  In  the  subtree  that  emanate 
from  It  (see  Leimna  A2.1.1).  Now  any  subtree  whose  root  node  lies  on  S contains 
at  least  one  active  source,  l.e.,  source  S.  Therefore,  the  probability  that  a 
node  of  level-i,  that  Is  on  S will  be  used  Is  C(y/2^).  Now  as  in  Theorem 
A2.1.1  use  the  fact  that  expectation  of  a sta  equals  the  sum  of  the  expecta- 
tion to  obtain. 


X.  - 1 + I e(y/2^) 

1-1 


(A2.6.12) 
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Now  we  will  make  several  observations  before  deriving  Y . Figure  A2.6.1 

8 

should  be  helpful  In  conceptualizing  the  following  observations. 

A)  T®  lies  above  S Iff  ■ !• 

B)  The  set  of  all  the  that  lie  above  S equals  all  the  nodes  that 

lie  above  S. 

C)  The  average  number  of  nodes  In  T^  that  are  used  (not  necessarily 
before  S Is  transmitted)  Is  D(y/2^).  This  follows  from  Theorem  A2.1.1. 

D)  A tree  T^  lying  above  S will  be  processed  before  S If  there  are  more 
than  one  active  source  In  the  subtree  of  level-1  that  contains  S. 

This  Is  so  because  If  S Is  the  only  active  source  In  T^,  then  that 
source  would  have  been  transmitted  at  some  node  at  level  less  than  1. 
Therefore,  the  probability  that  T^  that  lies  above  S will  be  processed 
before  S Is  6(y/2^). 

Next  we  will  combine  the  above  observations  so  as  to  derive  an  expression 

for  Y . 

s 

From  observations  A and  D we  have 

lies  above  S and  it  Is  processed  before  S]  > s^CCy/Z^***^)  (A2.6.13) 


From  Eq.  (A2.6.13)  and  observation  C we  have 


E[nuaber  of  nodes  in  . that  are  processed  before  S] 


- S^0(y/2^'^^)D(v/2^‘^^) 


From  Eq.  (A2.6.14)  and  observation  B ue  have 


00 

Yg  - I 8^0(y/2^'^^)D(y/2^‘^^) 
i-0 

Therefore  substituting  Eqs.  (A2.6.12)  and  (A2.6.15)  Into  Eq.  (A2.6.10) 
and  then  rearranging  we  have 

00 

E[d^|y]  - 1 + Z (l+8^_^D(y/2^))0(y/2^) 

1-1 

similarly 

00 

Eld_|y]  - 1 ^ (l-*^.  ,D(y/2^))0(y/2^) 

* 1-1 


^lEld^lyl+Eldglhll  - 1 + I ^(y/2S)e(y/2S 

1-1 

Since  8.  + 7.  - 1,  4tEtd  1+E[d  I]  Is  Independent  of  S,  therefore, 

11  A 9 9 


(A2.6.14) 


(A2.6.15) 


(A2.6.16) 


(A2.6.17) 


(A2.6.18) 
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(A2.6.19) 


E[d|y]  - 1 + Z 0(y/2S[l4D(p/2^)] 
1-1 


This  concludes  this  Appendix.  In  Appendix  A2.7  that  follows,  we  develop 
the  relationship  between  E[2.|y]  and  E[d|yl. 


A2.7  Relatlonahlp  Between  Efdlul  and  Effclul 


In  Appendices  A2.1  and  A2.6  we  derived  expressions  for  Etlly]  and 


E[d|vi].  These  expressions  are  retnrltten  below. 

E(8.1wl  - 1 + 2D(y/2) 

00 

El<i|y]  - 1 + I e(u/2^)[l  D(u/2^)] 

1-1 


(A2.7.1) 

(A2.7.2) 


where 


D(y)  - I 2H(u/2h 
1-0 


(A2.7.3) 


C(U)  - 1 - e"'^  -ye"^' 


(A2.7.4) 


e(y) 


(A2.7.5) 


In  this  Appendix  we  will  prove  In  Theorems  A2.7.1  and  A2.7.2,  respectively 
that 


Eldjy]  ^ Y E[l|yl  for  y 2;^  0 


(A2.7.6) 


and 
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E[<i|y]  < a EtJ!.|y]  + b for  y ^ 0 


where  (a,b)  * (.55,. 321) 

Equation  (A2.7.7)  la  specifically  proved  for  the  constants  (a,b)  > (.55, .321). 
However,  It  should  be  noted  that  there  are  other  (a,b)  pairs,  that  satisfy 
Eq.  (A2.7.7). 

Now  we  are  ready  to  prove  Eqs.  (A2.7.6)  and  (A2.7.7).  In  the 
Theorems  A2.7.1  and  A2.7.2  that  follow,  we  will  refer  extensively  to  the 
properties  of  D(y)  that  were  developed  In  Appendix  A2.2.  We  will  also  need 
the  following  Lemma. 

Lemma  A2.7.1.  Let  E[d|y]  be  as  given  in  Eq.  (A2.7.2).  Then 

•|^Etd|y]>0  (A2.7.8) 


Proof:  Differentiating  Eq.  (A2.7.2)  and  then  rearranging  we  have 


Now,  It  can  easily  be  shown  that  -1  + y + e~^  > 0 for  y > 0.  And 
since  D(y)  ^ 0 and  D'(y)  > 0 (See  Property  3B),  It  follows  that  all  the 
terms  In  Eq.  (A2.7.9)  are  positive. 

QED 

Theorem  A2.7.1.  Let  E[d|y]  and  E[l|y]  be  as  given  In  Eqs.  (A2.7.1)  and 
(A2.7.2)  then, 

E[d|y]  > Y E[£|y]  for  y > 0 (A2.7.10) 


Proof:  By  computation  It  can  be  shown  that  the  maximum  of  E[£|y]  > 

E(d|y]  In  the  range  of  ye[0,41  occurs  at  y-0  and  It  equals  0.  Therefore, 

El£|y]  - E(d|y]  < 0 for  0 < y £ 4 (A2.7.11) 

For  y ^ 4 we  will  first  prove  that  for  any  y > 4 there  exists  a 

ff 

y e[2,4]  such  that 

E[4|y]-E(d|y]  < | E[l|y*}-Eld|y*]  for  y > 4 (A2.7.12) 

The  conclusion  to  the  proof  will  follow  by  using  Eq.  (A2.7.11)  to  upper 
bound  the  right  side  of  Eq.  (A2.7.12). 
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4.. 


From  Eqs.  (A2.7.1)  and  (A2.7.2)  and  Properties  2A  and  2D  we  have 


i E[E|u]-E[dIu]  - - I + Y I (A2.7.13) 

1-1 


where 

" j \yj  (A2.7.14) 


Now  for  any  y > 4 there  exists  a k such  that  2 < \i/2  £ 4.  For  that  k 
Eq.  (A2.7.13)  can  be  rewritten  as  follows. 

k-1 

YlE[a|u]-E[d|yl]  - |lE[E|u2'VEtd|u2"’"l]  + I f^iu/2^)  (A2.7.15) 

1-1 


From  the  definition  of  k It  follows  that  y/2^  > 4 for  1 < k.  Next,  therefore, 
we  will  show  that  f^(x)  £ 0 for  x > 4,  our  ultimate  objective  being  to  show 
that  the  sum  In  Eq.  (A2.7.1S)  Is  negative  for  y ^ 4. 

From  Property  3D  we  have 


fl(x)  £ 


C(x)  + (3^x) 


(A2.7.16) 


It  Is  a straightforward  exercise  to  show  that  the  right  side  of  Eq.  (A2.7.16) 

Is  negative  for  x ^ 4.  So 

fj^(x)  <0  for  X > 4 (A2.7.17) 

And  from  Eq.  (A2.7.1S)  we  have 


i[E[E|y]-E[d|ii]]  < ^ElE|y2"*‘]-E[dlv2’’']] 


(A2.7.18) 


where 


y > 4 and  2 < 


y2"*‘  < 4 


QED 

Theorem  A2.7.2.  Let  E[£|y]  and  E(d|y]  be  as  given  In  Eqs.  (A2.7.1)  and 
(A2.7.2).  Then 

Eldly]  < .55  E[£|y]  + .321  (A2.7.19) 
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Proof:  The  structure  of  this  proof  is  similar  to  that  of  Theorem  A2.7.1. 
By  computation  It  can  be  shown  that 


Etdly]  - .55  E[ll|y]  < .321  for  0 < y < 8 


(A2.7.20) 


To  prove  Eq.  (A2.7.19)  for  y > 8 we  have  from  Eqs.  (A2.7.1)  and  (A2.7.2) 
and  Properties  2A  and  20  that 


E[d/y]  - .55  Etl/y]  - .45  + I 

1-1 


where 


(A2.7.21) 


(A2.7.22) 


As  In  Theorem  A2.7.1  we  conclude  this  proof  by  showing  that  f2(M)  ^ 0 for 
y ^ 8.  This  follows  by  substituting  the  following  lower  bound  for  D(y) 

Into  Eq.  (A2.7.22). 

D(y)  > 64(A)ty-4]  + 0(4) 

- 1.44(y-4)  + 4.77  (A2.7.23) 

QEO 


This  concludes  Appendix  A2.7. 
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fil. 

A2.8  An  Upper  Bound  to  E{e 

8li 

Let  G(s,p)  • E{e  (y}  be  the  generating  function  for  the  number  of 
steps  required  to  resolve  a conflict  by  the  binary  tree  algorithm,  given 
that  the  number  of  contending  packets  is  a Poisson  random  variable  with  mean 
y.  In  this  appendix  we  will  upper  bound  G(8,y)  and  9G(8,3)/3s.  More  speci- 
fically if  we  will  let 


G^(8,y.x)  - (I4y)e®"^  - (A/B^  + Ay/B)e“*'  + (A/B^)  (A2.8.1) 

then 

G<«.y)  < G^(8,y,x)  (A2.8.2) 

and 

h - h (A2.8.3) 

where 


A.-i 


, B - 1 < , < J, 


8 depends  on  x and  is  positive  for  1 < x < 3;  at  x-2,  s -ln2.  x is  a 
“ o 

variable  which  is  chosen  so  as  to  minimize  the  right  side  of  Eq.  (A2.8.1). 
Note  that  the  right  side  of  Bq.  (A2.8.1)  Is  decreasing  in  x.  therefore,  its 
greatest  lower  bound  occurs  at  xp3. 

Now  we  are  ready  to  begin  with  the  proof  of  Eq.  (A2.8.2).  Let, 

^j<«)  ■ packets  to  resolve}  (A2.8.S) 
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\ 

\ 


(A2.8.6) 


j-1 

- I eVJ(J)  fj^(s)fj_^(8)  + 

1-1 

We  have  divided  j Into  1 packetu  In  the  upper  subtree  and  j-1  In  the  lower. 
Note  that  e la  the  contribution  of  the  root  node,  2 ''(^}  Is  the  probability 
of  the  assumed  division,  and  f^(s)  and  f^_^(s)  are  the  generating  functions 
for  the  number  of  steps  to  resolve  each  subtree.  Furthermore,  note  that 
f^(s)  - f^(s)  - 1;  this  Is  because  exactly  one  step  Is  required  when  there 
are  either  zero  or  one  packet. 

Now  solving  Eq.  (A2.8.6)  for  f ^ (s)  we  have, 

j-1 

I eVJ(J)f^(s)fj^^(8) 

f,(8)  - ^ for  J > 1 (A2.8.7) 

j 1_2 

Solving  the  above  equations  recursively  we  have 


f^Cs)  - 1 


f2<->  " 

^ 2-e® 


f,(8)  - (-S^)(^) 

■*  2-e*  4-e* 


Motivated  by  the  above  expressions,  we  will  show  that 


f (g)  < .(x-l)e®  ^ ^ 

J'^  — _ s^  g ) l<x<3  and  s < s 

z-e  x_g  — c 


(A2.8.8) 


where  s^  > 0 for  1 < x < 3. 
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Assuming  that  the  above  Is  true  for  j-1,  then 


fj(.)  i 


8 2 . 8 j-4 

I e*2‘^d)(-S— ) ) + 

1-2 


2.’32-l  «* 


s j-3 


2-e*  x-e® 


l-2"J''‘^e® 


O s 8 

2-e  x-e 


2-e" 


x-e® 


l-2-:‘+le® 


(A2.8.9) 


(A2.8.10) 


for  j>2  l<x<3  and  0 < s < s^, 

— o 

where  8^(x)  > 0 for  1 < x < 3. 

(A2.8.11) 

Equation  (A2.8.11)  follows  from  Eq.  (A2.8.10),  by  noting  that  at  s«0;  both 
equal  one,  both  are  continuous,  and  of  Eq.  (A2.8.9)  is  greater  than  that 
of  Eq.  (A2.8.8). 

Finally 


2-e® 


, -.8  j-2 

(jjEzlle-) 


x-e 


E{e®*’|y)  - e®p(l-0,l|y)  + I f j (8)p(t-J |y)  (A2.8.12) 

J-2 

< e*(l-hi)e"'*  + i AB^"^  (A2.8.13) 

3-2 

where  A and  B are  given  in  Eq.  (A2.8.4).  Carrying  out  the  indicated  sum- 
mation in  Eq.  (A2.8.13),  we  obtain  Eq.  (A2.8.2). 
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•1 

rT 

X 


A tighter  bound  to  G(s,y)  should  be  possible  If  we  use  for  an  upper 


bound  to  f^(s)  the  following 


fj(s)  < (■ 


:)( 


^ 8 , 8 
2-e  4-e 


3e”®  ^^(x-l)e^^^  ^ 


s 


x-e 


for  j^3,  O^s^s^. 

(A2.8.1A) 


The  variable  x Is  chosen  so  that  the  Induction  step  Is  satisfied  for  all 
j ^ 3 and  for  Sq  > 0.  Although  we  have  not  carried  out  the  derivation 
of  this  upper  bound.  It  Is  felt  that  the  resulting  lower  bound  to  the  maxi- 
mum arrival  rate  at  which  the  moments  of  the  delay  are  finite,  will  be 
greater  than  that  obtained  from  the  bound  given  In  Eq.  (A2.8.2)  (see  Section 
2.4). 

Now  we  will  prove  Eq.  (A2.8.3).  From  Eq.  (A2.8.12),  we  see  that  It 
Is  sufficient  to  show  that 


3f^(8) 

38 


J-1 


2-e“ 


s j-1 
2-e® 


for  j > 1 


(A2.8.15) 


< 


(A2.8.17) 


CHAPTER  3 


OPTIMUM  DYNAMIC  TREE  ALGORITHMS  WITH  POISSON  SOURCE  MODEL 

3.1  Introduction 

I 

In  Chapter  2 we  presented  and  analyzed  the  static  binary  tree  algorithm, 
a protocol  where  the  tree  Is  fixed  to  be  binary  and  Independent  of  traffic  \ 

condition.  In  this  chapter  we  will  consider  a dynamic  algorithm  where  the 

1 

tree  Is  allowed  to  vary  from  epoch  to  epoch  depending  on  traffic  conditions. 

Even  though  the  tree  may  vary  from  one  epoch  to  another.  It  Is  held  fixed 
within  any  one  epoch.  The  tree  search  Is  the  same  as  that  of  the  static 
algorithms  that  were  considered  In  Chapter  2 and  may  be  carried  out  serially 
or  In  parallel,  deterministically  or  randomly.  The  source  model  that  will 
be  assumed  In  this  chapter  Is  the  Poisson  Source  Model. 

Now  we  will  state  the  Issues  of  this  chapter  more  precisely.  Three  main 
problems  will  be  considered  here.  The  first  Is  the  determination  of  the 
optimum  tree  to  be  used  In  optimum  tree  for  is  defined  to 

be  that  tree  which  minimizes  the  expected  number  of  slots  used  In 
given  the  observation  of  the  transmission  process  up  to  the  end  of  . The 
second  problem  will  be  the  analysis  of  this  optimal  dynamic  protocol,  l.e., 
the  determination  of  the  delay,  throughput  and  stability  properties.  The 
third  Is  the  optimization  and  analysis  of  a suboptimum  algorithm  where  the 
tree  Is  restricted  to  have  binary  nodes  everywhere  except  for  the  root  node; 
the  degree  of  the  root  node  Is  constrained  to  be  a power  of  two.  As  will 
be  seen,  this  algorithm  does  have  certain  Implementation  advantages. 
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4 


Given  chat  the  number  of  arrivals  In  any  one  slot  Is  a Poisson  random 
variable  vhlch  Is  Independent  of  the  arrivals  In  any  other  slot.  It  follows 
that  the  only  quantity  from  the  transmission  process  that  Is  needed  to  fully 


If  we  let  be  Che  degree  of  the  optimum  tree  at  depth-1,  Chen  In 
Section  3.2  It  Is  proved  that 


1 y < 1.70 

n 1.70  + 1.15(n-2)  < y < 1.70  + 1.15(n-l) 


g^  ■ 2 for  i ^ 1;  all  y 


(3. 1.0. 2) 


In  terms  of  Eqs.  (3. 1.0.1)  and  (3. 1.0. 2),  the  optimum  dynamic  algorithm  may 
be  stated  as  follows: 


characterize  the  packets  to  be  processed  during  is  hj  - the  number  of 

slots  In  Cj.  This  Is  because,  given  hj,the  number  of  packets  that  arrived  In 
Cj  Is  a Poisson  random  variable  with  mean  y where. 


y ■ Ah 


J 


(3. 1.0.1) 


t \ 


4.^ 


1.  Observe  h,  Che  number  of  slots  In  Che  previous  epoch 

2.  Calculate  y from  Eq.  (3. 1.0.1) 

3.  Determine  the  optimum  tree  to  be  used  In  the  following  epoch 
from  Eq.  (3. 1.0. 2) 

4.  Execute  Che  search  In  Che  following  epoch  using  one  of  Che  tree 
search  algorithms  of  Chapter  2. 
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In  Section  3.3  we  will  analyze  the  slightly  suboptlmum  but  easier  to 

Implement  tree,  where  Is  restricted  to  be  an  even  positive  Integer.  There, 

we  will  calculate  upper  and  lower  bounds  to  the  average  delay  as  a function 

of  the  arrival  rate  (see  Fig.  3. 3. 2. 2),  prove  that  the  maximum  average 

throughput  Is  .430  packets/slot,  and  show  that  system  Is  stable  for  X < .430, 

in  the  sense  that  all  the  moments  of  the  delay  exist  for  X < .430. 

The  third  problem  that  was  posed  above  will  be  considered  In  Section  3.4. 

There  we  will  show  that  If  the  root  node  degree  of  a binary  tree  Is  constrained 
K * 

to  be  2 ; then  K , the  optimum  K,  Is  given  by 


1 for  p _<  3.40 

K 3.40  (2*^"^^  < y < 3.40(2^"^)  (3. 1.0. 3) 

We  will  also  determine  upper  and  lower  bounds  to  the  E[ delay]  (see  Fig. 

3. 4. 2.1)  and  we  will  show  that  for  this  dynamic  algorithm  the  maximum  E[thr] 
Is  less  than  .430  but  greater  than  .420  packets/slot. 

In  Appendix  A3  that  accompanies  this  chapter,  we  prove  several  of  the 
theorems.  This  appendix  Is  consulted  frequently  as  we  proceed  through 
this  chapter. 
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3.2  The  Optimum  Tree 


In  this  section  we  will  prove  that  If  the  number  of  contending  packets 
Is  a Poisson  random  variable  with  mean  y,  then  the  symmetric  tree  that 
minimizes  the  expected  number  of  slots  required  to  process  the  contending 
packets  Is  given  by  Eq.  (3. 1.0. 2).  Althougn  we  do  not  prove  that  symmetric 
trees  are  optimum,  that  this  Is  so  should  be  obvious  from  the  symmetry  of 
the  problem. 

It  can  be  shown,  by  a procedure  analogous  to  that  of  Appendix  A2.1.1, 
that  the  expected  number  of  slots  required  to  process  the  packets  Is  given 
by 


E{h|y} 


1-0 


(3. 2. 0.1) 


where 


and 


U 

*0*1*"*1 


(3. 2. 0.2) 


CCM)  - 1 - e"**  - Me"** 


(3. 2. 0.3) 


In  terms  of  the  above  three  equations,  our  goal  In  this  section  Is  to 
determine  that  {g^;  1 - 0,1,2,...}  which  minimizes  E{h|y}  subject  to  the 
constraint  that  g^  Is  an  Integer  and 


i 


1 If  1 - 0 

2 otherwise 


(3. 2. 0.4) 
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The  solution  will  be  carried  out  In  three  steps.  First,  we  will  show  that 
g*,  the  optimum  g^,  equals  2 for  sufficiently  large  1.  This  Is  accomplished 
In  Section  3.2.1.  Secondly,  In  Section  3.2.2,  we  will  prove  that  the  minimum 

of  E{hlu}  over  and  g^,  subject  to  the  condition  that  g^^  ■ 2 for  1 > k, 

* 

occurs  at  then  g^^  « 2.  This  result  applied  recursively  with  the  result  of 

* 

Section  3.2.1  as  the  boundary  condition  proves  that  g^  * 2 for  1^1.  Finally, 
* 

In  Section  3.2.3,  g^  Is  calculated  as  a function  of  y;  this  result  Is  given 
by  Eq.  (3. 1.0. 2). 


3.2.1  Optimum  Tall  End  of  the  Tree 

We  begin  this  section  with  the  following  definitions: 


^ ^®k’  *lri-l’ ^ ° 


(3. 2. 1.1) 


= the  G^  which  minimizes  E{h|y}  subject  to  (3. 2. 1.2) 

g^  (1  ^ k)  being  an  Integer  and  greater  than  1 


2 i [2,2,2 ] 


(3. 2. 1.3) 


(3. 2. 1.4) 


««k) 


. 1-1  _ -1 
i ^ *k+l^  " ^ ^ 

1-0  j-0  j-0 


(3. 2. 1.5) 


3 * 

In  this  section  we  will  prove  that  If  k > 0 and  If  ^ 2g  then  G^^  - 

It  should  be  observed  that  since  » y/gQg^...g^  and  8^  ^ 2 for  1 > 0, 
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In  other  words,  the 


3 

then  for  any  finite  y,  a k exists  such  that  ^ ^ 

main  result  of  this  section  states  that  beyond  some  finite  depth  the  optimum 
tree  has  binary  nodes. 

The  solution  will  be  carried  out  as  follows: 


a) 

b) 

c) 


Show  that  If  minimizes  It  will  also  minimize 

E{h|y}.  This  Is  Theorem  3. 2. 1.1. 

In  terms  of 

fi(G^).  This  Is  accomplished  In  Theorem  3. 2. 1.2. 

Show  that  the  upper  bound  of  K evaluated  at  G^^  >■  2^  Is  less  than  the 
lower  bound  at  G^  ^ 2.  This  Is  accomplished  In  Theorems  3. 2. 1.3 
and  3. 2. 1.4  and  Corollary  3. 2. 1.1. 


Determine  upper  and  lower  bounds  to  H(G^,Mj^_ 


I 


Theorem  3. 2. 1.1:  If,  Gj^  minimizes  K(Gj^,M^_j^)  for  all  then  it  will 

also  minimize  E{h|y}. 

Proof:  By  combining  Eqs.  (3. 2.0.1)  and  (3. 2. 1.4),  E{h|y}  may  be  written  as 
follows: 

g 

E{h|y}  - gQ  + y I eCM^)  + (3.2. 1.6) 

1-0  ^ 

Since  of  the  three  terms  on  the  right  side  of  Eq.  (3. 2. 1.6),  only  H depends 
on  G^  and  since  fi  and  are  greater  than  zero.  It  follows  that  If 

G^  minimizes  R then  It  will  also  minimize  E{h|y}. 

QED 
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Theorem  3. 2. 1.2:  Let  and  Gj^  be  as  defined  in 

Eqs.  (3. 2. 1.1)  through  (3. 2. 1.5)  and  (3. 2. 0.2).  Furthermore,  let 
0 £ ^ 1 8^  'J®  Integer  and  greater  than  one.  Then, 

(1  - |M^_^)fi(Gj^)  < (3.2.1. 7) 

-^1 

Proof : First  substitute  the  Taylor  series  expansion  of  e Into 
Eq.  (3.2.0. 3)  and  then  substitute  Eq.  (3. 2. 0.3)  Into  Eq.  (3. 2. 1.4)  and 
rearrange  to  obtain 


fi(G. 


00  00 

k’“k-l^  “ *k+i”k-l+l  ” *k+l^^\-l+i^  (3. 2. 1.8) 


where 


■ I (snr 

n“2 


(3. 2. 1.9) 


If  0 £ ^ *^^k-l+l^  alternating  sequence  of 

decreasing  magnitude  and  It  follows  that. 


« 1 i i 


(3.2.1.10) 


From  Eqs  (3. 2. 1.8)  and  (3.2.1.10)  we  have 


00 

fi(Gk.Hj^_l)  < 2m  ^ *kH\-l+l 
1-0 


(3.2.1.11) 
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and 


fi(G, 


...  C3.2.1. 


1 r 2 

^ at  ^ ^ A ~ 

K A 1»0 


Since  0 < for  ill.  It  follows  that 


OD  *® 

1 V „ Vl+i-  ^ i y g M. 

3 i *\i+i  - 3 ^ ®V.+i^-l+i 


i-0 


i-0 


(3.2.1.13) 


Substituting  Eq.  (3.2.1.13)  into  Eq.  (3.2.1.12)  we  have 


OB 

fi(G^,M^_l)  > (1  - I \_i)  aT— - I *lcfiVl+l 

*■1  1-0 


(3.2.1.14) 


Finally  from  Eqs.  (3.2.0.2)  and  (3.2.1.5)  we  have 


S(G 


OB 

\}  " aTT  ^ *fcfi\-i+i 

1-0 


(3.2.1.15) 


and  by  substituting  Bq.  (3.2.1.15)  into  Eqs.  (3.2.1.11)  and  (3.2.1.14).  we 
obtain  the  desired  result. 


QED 


12) 
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Theorem  3. 2. 1.3;  Let  l>«  *8  given  in  Eq.  (3. 2. 1.5)  and  let  be 

integer  and  greater  than  one.  Then  the  minimum  of  is  2 and  it  occurs 

•t  0^  - 2. 

Proof:  For  n > k we  have  from  Eq.  (3. 2. 1.5) 


S(\) 


n-k-1 

1 I 

1-0 


-1 

h(G 


(3.2.1.16) 


From  Eq.  (3.2.1.16)  it  follows  that  if  g^  minimizes  ^(Gj^)>  it  will  also 

. * 
minimize  n(G^).  Therefore  we  will  determine  g^. 

From  Eq.  (3. 2. 1.5)  we  have 


(3.2.1.17) 


Now  assume  that  has  been  determined  in  Eq.  (3.2.1.17)  and  we  would  like 

to  know  g*.  But  h(G  ) is  convex  in  g , therefore  the  optimum  g*  is  the 
n n n n 

smallest  Integer  greater  than  one  such  that 


fi<*n*Cl^  i Cl> 


(3.2.1.18) 


Substituting  Eq.  (3.2.1.17)  into  Eq.  (3.2.1.18)  and  rearranging,  we  have 

* 

the  following  equivalent  expression  that  g must  satisfy, 


*n<*n  * i 


(3.2.1.19) 
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But  since 


fi(G^l)  lfi(2)  - 2 


It  follows  that  the  smallest  Integer  g greater  than  one  that  satisfies 

n 


Eq.  (3.2.1.19)  Is 


* 

8n-2 


Corollary  3. 2.1.1;  Let  be  such  that  8^  ^ 3.  Then  ^ 13/6. 


Proof:  From  Eq.  (3.2.1.17)  we  have 


- r ^ -17- 


(3.2.1.20) 


and  from  Theorem  3. 2. 1.3  we  have 


fi(G,^l)  > 2 


(3.2.1.21) 


therefore 


2 

fi(G.  ) > 

k “2  g^ 


(3.2.1.22) 


Since  the  derivative  of  the  right  side  of  Eq.  (3.2.1.22)  Is  positive  for 


g^  > 2,  It  follows  that  for  i 3 
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or 

«<S>  i r 


(3.2.1.23) 


QED 

Now  we  are  ready  to  prove  the  main  result  of  this  section  In  the  following 
theorem. 

Theorem  3. 2.1. At  If  "Jg*  then  ^ minimizes 


Proof:  Take  any  G,  2 and  let  g for  n > k be  the  first  element  of  G,  that 
ic  tl  ^ 

is  greater  than  2.  Then  from  Eq.  (3. 2.1.4)  we  have. 


n-k-1 


vui 


*lcfl 


M 


n-1 


(3.2.1.24) 


First,  note  that  the  G„  that  will  minimize  fi(G_,M  ,)  will  also  minimize 

n n n—i 

R(Gk.Mj^_l).  Then  from  Theorems  2 and  3 and  Corollary  1,  we  have 


fi(2,  M^j^)  < 2 (3.2.1.25) 

# 

i - 3 "n-l’F  (3.2.1.26) 
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BuC  for  n ^ k,  0 ^ < 3/26,  therefore  from  Eq.  (3.2.1.26), 

fi(Gn,Mn_i)  > 2 (3.2.1.27) 

and  it  follows  from  Eqs.  (3.2.1.25)  and  (3.2.1.27)  that, 

fi(2,M„  -)  < fi(G^,M^  .)  (3.2.1.28) 

n— 1 n n^x 

Finally  since  n was  chosen  so  that  " 2 for  It  £ 1 £ n-1,  we  have  that 

M . 

fi(G^.\_l)  - fi(2,M^_^)  - f*'<Gn’“n-l^  " (3.2.1.29) 

and  from  Eq.  (3.2.1.28)  we  conclude  that 

fi(G^,K^_j^)  - fi(2,M^_j^)  > 0 for  (3.2.1.30) 

QED 

Theorems  3. 2. 1.1  and  3. 2. 1.4  taken  together  prove  the  main  result  of 
this  section,  i.a.,  tha  tall  end  of  the  optimum  tree  is  binary.  In  the 
next  section  we  will  be  concerned  with  the  degrees  of  the  nodes  of  the 
optimum  tree  from  the  tell  end  to  the  root  node. 

3.2.2  A Recursive  Optimisation  Technique 

Vfhereas  in  Section  3.2.1  we  were  concerned  with  the  tail  end  of 
the  optlmtmi  tree,  in  this  section  we  will  develop  a recursive  relationship 
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4 


which  In  conjuctlon  with  the  results  of  the  previous  section  proves  that 

* * 

> 2 for  1 > 0.  The  optimum  degree  g^  of  the  root  node  does  depend  on 

y and  this  will  be  considered  In  Section  3.2.3. 

The  recursive  relationship  that  will  be  proved  here  is  as  follows. 
Let  and  Integer  and  greater  than  one  (for  k«0,l...)  and  let 

g^  * 2 for  1 > k + It  then  for  any  y ^ 0,  the  point  that  minimizes  E{h|y} 
satisfies. 


- 2 (3. 2. 2.1) 

The  minimization  problem  has  a slightly  different  objective  for  k>0 
then  It  does  for  k > 0.  Therefore,  the  analysis  of  this  section  Is 
organized  Into  two  subsections.  In  Subsection  1 the  case  where  k > 0 Is 
considered  and  In  Subsection  11  the  case  where  k^O  Is  considered. 


1.  Recursive  Relationship  for  k > 0 

Now  we  begin  the  analysis  for  k > 0.  Setting  g^  ■ 2 for  1 ^ k 2 
In  Eq.  (3. 2. 0.1)  we  have 


E{h|y} 


*0 


k-2 


1-0 


CMj^)  + 


^^*k’*fc+l’*Hc-l^ 


(3. 2. 2. 2) 


where 


w M 

L(gk»gk+i*\_i)  = *k*k+l^^  g^^  ^ ^*k*k+l® ^*k*fcf  1^  (3. 2. 2. 3) 
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and 


D(x)  - I iHix/lh  (3. 2. 2. 4) 

1-0 

As  can  be  seen  from  Eq.  (3. 2. 2. 2),  the  that  minimizes  L will  also 

minimize  E{h|y}.  Therefore,  we  will  be  concerned  with  the  minimization  of 
over  la  order  to  simplify  the  notation  somewhat 

we  will  drop  the  subscript  of  M and  also  use  and  g2  Instead  of  and 
g^^.  In  what  follows,  therefore,  we  will  prove  that  for  any  M ^ 0,  the 
minimum  of  L(g^,g2,M)  occurs  at  g2  - 2.  This  will  be  accomplished  In  two 
parts.  In  Part  1 It  will  be  proved  for  M £ 8 and  In  Part  2 for  M > 8. 

Part  1.  Minimization  of  L(g^,g2,M)  for  M £ 8 

In  this  part  we  will  prove  that  If  M ^ 8,  and  If  g^  and  g2  are 

Integers  and  greater  than  one,  then  the  min  L(g. ,g.,M)  occurs  at  g.  - 2. 

®1*®2  ^ ^ ^ 

The  results  of  Section  3.2.1  prove  this  statement  for  M £ 3/26;  therefore, 
we  will  be  concerned  with  Its  proof  for  3/26  ^ M £ 8. 

First  divide  the  (g^  x g2)  space  Into  three  regions,  R^,R2,  and  R^  as 
follows. 


Rl-I(gl.g2)s  2 < g^  < 9,  2 < g2  < 5] 
^2  ■ t(8j^»82)5  8^  1 10,  g2  1 3] 

83  ■ 1(83,82)5  2 < gj^  £ 9,  g2  i 61 


(3. 2. 2. 5) 


(3. 2. 2. 6) 


(3. 2. 2. 7) 


(See  Fig.  3.2. 2.1) 
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2 4 6 8 10 


Figure  3.2. 


.1  The  Subdivision  of  the  (g.  x g.)  Space  Into  R. , R, 
and  Rj 
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j 

I 


4! 


In  the  Appendix  to  this  chapter  we  prove  In  Theorem  A3.1  that 


Lj^(gl.g2.M)  < L(gj^.g2,M)  for  M/(gj^g2)  < 1.2  (3, 2. 2. 8) 

where 


4(gl.g2.M)  = giCCM)  + gig,?(?-)  + 2M(  ” 


_ 

* Q 


g2^g2  ^ ^1®2 


y)  (3. 2. 2. 9) 


Note  that  the  above  lower  bound  Is  valid  In  R2  and  R^. 

In  Theorem  A3. 2 It  Is  proved  that  for  (g^tg2)  In  R2  and  .1  ^ M ^ 8, 
then  Lj^(gj^,g2,M)  Is  minimum  at  (gj^,g2)  ■ (10,3). 

In  Theorem  A3. 3 It  la  proved  that  If  (gj^,g2)  in  R^  and  .1  ^ m ^ 8, 
then  L2^(gj^,g2.M)  Is  minimum  at  g2  ■ 6. 

Now  for  M ■ .Ik,  k « 0,1 80,  the  quantity  L(g^g2,M)  was  minimized 

over  the  32  (g^,g2)  points  In  R^.  The  minimum  did  occur  at  g2  * 2 and  It  Is 
plotted  In  Fig.  3. 2. 2. 2 as  a function  of  M. 

In  R2>  because  of  Theorems  A3.1  and  A3. 2,  we  calculated  (10, 3,. IK) 

and  In  R^,  because  of  Theorems  A3.1  and  A3. 3,  we  calculated 


min  L, (g. ,6,.1K)  for  k - 0,1,..., 80 

218i19  ^ ^ 

These  results  are  also  Included  In  Fig.  3. 2. 2. 2.  As  can  be  seen  from  this 
figure. 
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^Min  L(g,,g2,m) 
9i92 

^Min  Lji{g„6,m) 

9i 


L«(10,3,m) 


(3.2.2.10) 


Min  L(g  ,2,M)  < L.(10,3,M) 
2<g^<9 


Min  L(g  ,2,M)  < L(g,,6,M)  ; 2 < g.  < 9 (3.2.2.11) 

2<g^<9 

for  .1  < M < 8 

Therefore,  we  conclude  Chat  L(g^,g2,M)  attains  Its  Integer  minimum  at  g2*2, 
if  gj^  ^ 2,  ^2  — ^ .1  £ m ^ 8. 

Part  2.  Minimization  of  L(g^,g2,M)  for  M ^ 8 

In  this  part  we  will  prove  that  if  g^  and  g2  are  integers  and 
greater  than  one,  then  for  any  M ^ 8,  the  minimum  of  L(g^,g2,M)  lies  on  the 
line  g2*2. 

In  this  range  of  M it  is  more  convenient  to  work  with  H(g^,g2,M)  Instead 
of  with  L(g^,g2,M),  where 

H(gj,g2.M)  H i L(gj^,g2,M)  (3.2.2.12) 

Since,  for  any  M > 0,  if  (g^,g2)*  minimizes  H(g^,g2,M)  it  will  also  minimize 
L(gj^,g2tM),  it  is  sufficient  to  prove  that  the  minimum  of  H(g^,g2,M)  for 

*1  — *2  — ^ *2“^* 

Here  again,  the  way  this  will  be  proved  is  to  show  that  an  upper  bound 
to  H(g^,g2,M)  evaluated  at  82*2  and  minimized  over  g^  ^ 2 is  less  than  a 
lower  bound  evaluated  at  {gj^  ^ 2,  82  i 3} . 
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since  ^(M)  is  an  increasing  function  bounded  by  1,  we  have  from  Eqs.  (3. 2. 2. 3) 
and  (3.2.2.12)  that. 


H^(gl,g2.M)  1 H(gj^,g2,M)  < H^(g^,g2,M)  ; for  M > 8 (3.2.2.13) 


where 


*1  . *1*2  „ r ^ *1*2...  M 


H^(*1.*2.M)  - 5(8)jp  + C(^)  + 2 [ — ^ (3.2.2.14) 

1 1-0  2 *1*2 


*1  , *1*2  , . „ V ^ *1®2,,  M 


H^(gl.g2.M)  - iT  + + 2 I ) (3.2.2.15) 


1-0 


2 gj_g2 


Since 


Si  * M 

H„(gi.2.M)  -5^+2  I 


(3.2.2.16) 


1-0 


*1^ 


we  have  from  Theorem  A3. 4 that 


Min  H^(g^,2,M)  < 2.34 


*1>2 


(3.2.2.17) 


In  Theorem  A3. 5 It  Is  proved  that  If  M/g^  f.  82  — M ^ 8,  then 


H^(gj^,g2,M)  >2.498 


(3.2.2.18) 


-148- 


And  In  Theorem  A3. 6 It  Is  proved  that 


Hj^(gl,g2.M)  > 2.34  for  .1  £ M/g^^  £ 3 and  g2  > 3 (3.2.2.19) 

Therefore,  from  Theorems  A3. 4,  A3. 5 and  A3. 6 we  conclude  that  If 
M/gj^  < 3,  M > 8,  gj^  > 2,  g2  i 2 then 
the  minimum  of  H(g^,g2,M)  occurs  on  Che  line 

®2  - 2 

In  Theorem  A3. 7 It  Is  proved  that  the  minimum  of  H(g^,g2,M)  for  M/g^  > 3 also 
occurs  on 

g2  * 2 

This  concludes  Part  2. 


11.  Determination  of  Optimum  g^ 

What  we  proved  up  to  this  point  Is  chat  the  minimum  of  E{h|w}  over 
gj^  and  g^^  (subject  Co  Che  constraints  Chat  g^  and  g^^  are  Integer  and 
greater  chan  one;  g^  * 2 for  1 ^ k^2,  and  k > 0)  occurs  at  gj^^*2.  Ua  still 
have  to  prove  the  above  statement  for  k>0  and  this  la  what  we  will  do  next. 


-149- 


Setting  ■ 2 for  1 2^  2 In  Eq.  (3. 2.0.1)  we  have 


E{h|lj} 


(3.2.2.20) 


We  are  going  to  show  that  the  minimum  of  the  above  expression  subject  to 
gQ  and  gj^  being  Integer  and  Eq  ^ 1.  ^ 2 occurs  at  g^^  - 2.  By  comparing 

Eq.  (3.2.2.20)  to  Eq.  (3. 2. 2. 3)  It  can  be  seen  that  this  problem  Is  very 
similar  to  that  solved  In  Subsection  1.  Therefore,  one  may  go  through  a 
similar  procedure  In  solving  the  above  programming  problem.  Here,  however,  we 

•k 

are  going  to  be  less  rigorous  and  simply  show  that  g^  - 2 by  computation. 

The  procedure  to  be  used  here  Is  as  follows. 

Divide  Eq.  (3.2.2.20)  by  u and  let  x - W/Eq  to  obtain  f(x,gj^). 


f (x.gj^) 


i(l  + g^5(x)  + 2g^D(|-)) 


(3.2.2.21) 


Given  the  preceding  definition  of  f(x,g^).  It  follows  that  If  the 
minimum  of  f(x,gj^)  occurs  at  gj^‘2  then  the  minimum  of  E{h|M}  will  also  occur 
at  8}^"2.  We  computed  f(x,g^)  for  various  values  of  (x,g^)  and  the  results 
are  shown  In  Fig.  3. 2. 2. 3. 

Now  If  y < 3.08  then  x ■ ^ < 1.54,  and  as  can  be  seen  from  this 

*0 

figure 


f(x,2)  < f(x,gj^)  for  X < 1.54;  > 2 (3.2.2.22) 

If,  on  the  other  hand,  y ^ 3.08  then  It  can  be  shoini  that  a g^  ^ 7 exists 
such  that  .89  £ y/g^^  £ 1.54.  But  here  again  as  can  be  seen  from  Fig.  3. 2. 2. 3 
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(3.2.2.23) 


f(x,2)  ^ f(x,g^)  for  .89  X ^ 1.54;  8^  > 2 

Therefore,  we  will  assume  that 

g*  - 2 (3.2.2.24) 

This  concludes  Section  3.2.2. 

3.2.3  The  Root  Node  of  the  Optimum  Tree 

In  this  section  we  will  determine  the  optimum  g^,  l.e.,  that  g^ 
which  minimizes  EChlu.g^}  subject  to  gg  being  Integer,  gQ  ^ 1 and  g^  > 2 
for  1^1. 

We  begin  by  setting  In  Eq.  (3. 2. 2. 2)  to  obtain 

E{h|y,  gg)  - ggd  + 2D(u/gQ)}  (3. 2. 3.1) 

This  equation  was  calculated  and  the  results  are  Illustrated  In  Fig.  3. 2. 3.1. 
The  problem  posed  In  this  section  Is  solved  by  determining  0(n),  a quantity 
defined  by, 

g*(u)  - n for  0(n)  < U £ U(n+1)  (3. 2. 3. 2) 

As  can  be  seen  from  Fig.  3. 2. 3.1,  0(gQ)  la  that  u which  satisfies 

E{h|y,gQ-l}  - E{h|y,  gjj}  (3. 2. 3. 3) 
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E[h//x.,go]  (slots) 


Equation  (3. 2. 3. 3)  was  solved,  and  U(n)  to  within  two  significant  places  to 
the  right  of  the  decimal  point  Is  given  by. 


0(n)  - 1.70  + 1.15(n-2) 
This  concludes  Section  3.2.3. 


(3. 2. 3. 4) 


I 


3.3  Analysis  of  the  Dynamic  Algorithm 


In  this  section  we  will  analyze  (for  the  Poisson  Source  Model)  the 
dynamic  tree  algorithm  where  all  the  nodes  except  for  the  root  node  are 
binary  - the  degree  of  the  root  Is  restricted  to  be  even  and  It  Is  chosen 
so  as  to  minimize  the  expected  number  of  slots  used,  given  the  expected 
number  of  contending  sources.  This  algorithm  should  be  ognlzed  as  being 
very  similar  to  the  optimum  algorithm  of  Section  3.2,  the  only  difference  being 
chat  here  Che  root  node  Is  restricted  to  taking  on  an  even  degree  where  In 
Che  optimum  algorithm  It  Is  not.  There  are  two  reasons  for  considering 
this  suboptimum  algorithm.  First,  the  even  degree  root  node  tree  Is  easier  to 
Implement  since  here  all  nodes  are  even,  and  secondly  the  analysis  Is  neater, 
even  though  It  Is  not  essentially  less  complex. 

This  algorithm  will  be  executed  serially  by  assigning  two  consecutive 
slots  to  It  for  each  rotmd  trip  Interval.  If,  for  example,  Che  Initial  node 
has  degree  2r,  Chen  the  tree  search  Is  Identical  to  r consecutive  serial 
searches  (of  Che  type  described  la  Chapter  2)  idiere  each  of  the  r searches 
Is  over  1/r  of  the  sources. 

As  In  Chapter  2,  an  algorithmic  step  consists  of  the  actions  taken  In 

two  consecutive  slots.  Therefore,  the  root  node  correponds  to  r steps  but 

all  ocher  nodes  correspond  to  a single  step  each.  It  also  follows  that  If 

we  let  d be  the  number  of  steps  between  the  arrival  and  the  successful 

transmission  of  a packet,  and  T and  T the  round  trip  delay  and  length 

r s 

of  one  slot,  Chen  the  delay  (In  seconds)  experienced  by  that  packet  Is 
given  by 

I 

i 
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(3. 3.0.1) 


packet  delay  ■ 

Since  in  the  above  equation  and  are  constant.  It  follows  that  In  order 
to  characterize  the  packet  delay,  one  need  only  obtain  the  statistics  of 
5.  In  this  section,  therefore,  we  will  obtain  upper  and  lower  botinds  to 
E{5}  as  a function  of  the  packet  arrival  rate  X.  A by-product  of  this 
analysis  Is  the  determination  of  the  maximum  throughput  and  the  characteri- 
zation of  stability. 

Since  the  expected  delay  with  the  parallel  execution  of  the  algorithm 
Is  less  than  that  with  the  series  execution,  the  delay  results  obtained  here 
may  be  considered  to  be  upper  bounds  to  those  of  the  parallel  execution.  The 
maximum  throughput,  however,  is  the  same  for  both  schemes. 

The  analysis  is  organized  as  follows.  In  Section  3.3.1,  expressions 
2 

for  E{j.|y,r},  E{1  |vi,r}  and  E{d|y,r}  are  derived.  The  quantities  1,  y and  d 
have  the  same  definition  as  they  did  In  Chapter  2;  that  Is,  y Is  the 
expected  number  of  packet  arrivals  In  the  previous  epoch,  I Is  the  number 
of  algorithmic  steps  In  the  present  epoch,  and  d Is  the  number  of  algorithmic 
steps  until  a randomly  chosen  packet  from  the  set  of  conflicting  packets  is 
successfully  transmitted.  Furthermore,  note  that  the  above  conditional 
moments  with  r>l  are  Identical  to  the  corresponding  ones  of  Chapter  2.  In 
Section  3.3.2,  we  determine  the  optimum  r and  obtain  upper  and  lower  bound 
to  E{5}  (see  Eqs.  (3. 3. 2. 3),  (3.3.2.13),  and  (3.2.2.14)).  In  Sections  3.3.3 
and  3.3.4  the  stability  and  throughput  are  considered,  respectively.  It  should 
be  noted  that  the  analysis  of  the  dynamic  tree  is  very  similar  to  that  of  the 
static  binary  tree  which  was  considered  In  Section  2.2.  Therefore,  since  many 
of  the  results  of  Section  2.2  carry  over  to  the  dynamic  tree  analysis,  we  will 
not  be  as  detailed  here. 
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3.3.1.  Derivation  of  E{£.|ii.r}.  E{)l^lu.r}  and  E{dlu.r} 

Consider  a tree  where  the  root  node  has  degree  2r  and  where  all  other 
nodes  have  degree  two.  Furthermore,  the  slots  are  used  In  pairs,  and  the 
algorithm  Is  executed  serially  as  described  above.  In  this  section,  expressions 
of  E{2,|u»r},  E{ll^|v,r}  and  E{d(vi,r}  will  be  derived  In  terms  of  E{f,|^,l}  and 
E{)tr^  1-^,1}  (the  corresponding  quantities  of  the  binary  tree  that  was  analyzed 
In  Section  2.2.  These  derivations  are  developed  In  Theorems  3. 3. 3.1  and 
3. 3. 3. 2 below. 

Theorem  3. 3. 3.1;  Let  2.  and  y be  as  defined  above  and  let  the  source  model 
be  Poisson.  Then 

E{Ely,r}  - rE{)ll^,l}  (3.3. 1.1) 

E{4^1u,r}  - r(r-l)E^{£|^,l}  + rE{4^|^,l}  (3.3. 1.2) 


or 


E{E^ly,r}  - rV{4|^,l}  + r Var{i|^,l} 


(3.3.1. 3) 


Proof:  Divide  the  2r  subtrees  that  emanate  from  the  root  node  Into  r pairs 
and  let  be  the  number  of  nodes  that  are  visited  In  the  I'th  pair.  Then 
since  the  root  node  corresponds  to  r steps  and  all  other  nodes  correspond  to 
a single  step,  we  have 


t - I (1  + 1^)  (3.3.1.4) 

1-1 
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Taking  expectations  of  the  above  equation  and  then  noting  that  1 + equals 
the  number  of  nodes  in  a binary  tree  with  mean  u/r  we  have  Eq.  (3. 3. 1.1). 

To  prove  Eq.  (3. 3. 1.2),  square  both  sides  of  Eq.  (3. 3. 1.4)  to  obtain. 


r-1  r r 

~ 2 I I (1  + )l^)(l  + £j)  + I (1  + (3. 3.1. 5) 

i«l  j-i+1  1-1 


Taking  expectation  of  Eq.  (3. 3. 1.5)  and  noting  the  Independence  of 
for  1 j we  have  Eq.  (3. 3. 1.2). 


and 

1 


QED 


Theorem  3. 3. 1.2;  Let  the  source  model  be  Poisson,  and  E,  d,  u be  as  defined 
above,  (Note  that  d Is  the  same  as  d^).  Then 


E{d|y,r}  < .55E{]lly,r}  + .321 


(3. 3. 1.6) 


and 


E{d(y,r}  ^ j E{il|M,r} 


(3. 3. 1.7) 


Proof:  Here,  again,  divide  the  tree  into  r pairs  of  subtrees  and  observe 
that  each  pair  Is  equivalent  to  a binary  tree  where  the  expected  number  of 
contending  packets  Is 

Assuming  that  each  pair  Is  processed  In  sequence,  we  have 

E{d|vi,r,  test  packet  Is  In  the  I'th  pair}  - (l-l)E{E|^,l}  + ECdl-^,!} 

(3. 3. 1.8) 
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Next  multiply  by  — and  sum  over  1 to  obtain 


E{d|n,r}  - ~ E{Jt|^,l}  -f  E{d|^,l} 


(3. 3. 1.9) 


But  from  Appendix  A2.4.2 


E{d|y,l}  < .55E{2ly,l}  + .321 


(3.3.1.10) 


and 


E{d|y,l}  E{2iy,l} 


(3.3.1.11) 


Substituting  Eq.  (3.3.1.10)  Into  Eq.  (3. 3. 1.9)  we  have 


E{d|^,l}  r E{1|^,1}  + .321 


(3.3.1.12) 


But  for  r > 1, 


< 55 

2r  - •” 


(3.3.1 


.13) 


Therefore,  substituting  Eq.  (3. 3.1.1)  Into  Eq.  (3.3.1.12)  and  then  applying 
Theorem  3. 3. 3.1,  we  get  Eq.  (3. 3. 1.6).  To  obtain  Eq.  (3. 3. 1.7),  substitute 
Eq.  (3.3.1.11)  Into  Eq.  (3. 3. 1.9)  and  apply  Theorem  3. 3.1.1. 


3.3.2  Upper  and  Lower  Bounds  to  Average  Dela'^ 


Here,  upper  and  lower  bounds  to  E{6},  that  are  functions  of  X,  will 
be  obtained  for  the  dynamic  tree.  This  will  be  accomplished  as  follows. 
First,  bounds  to  E{5}  that  depend  on  E{Jl}  and  E{£  } are  obtained.  Secondly 

it 

r (u),  the  optimum  relationship  between  y and  r,  is  obtained.  (Note:  2r  is 

the  degree  of  the  root  node.)  And  finally,  r*(u)  is  substituted  into 
2 

E{£.|y,r},  E{2  |y,r}  and  E{d|y,r}  so  that  techniques  similar  to  those  used 

in  Section  2. 2. 1.3  may  be  applied  here  to  obtain  bounds  to  e{2.}  and 
2 

E{£  } that  are  functions  of  X. 


E{6}  vs  E{il}  and  E{£^} 

Since  Eqs.  (3. 3.1. 6)  and  (3. 3. 1.7)  are  identical  to  Eqs.  (2.2.2.16) 
and  (2.2.2.17),  the  results  of  Section  2.2.2  (Eqs.  (2.2.2.20)  and  (2.2.2.25)) 
are  applicable  here.  Therefore, 


E{6}  < 1.05 


E{Jl^} 

“eTIT 


+ .321 


(3. 3. 2.1) 


E{6}  > + (3. 3.2. 2) 

Determination  of  r (y)  for  Even  Degree  Node  Tree 

* 

The  quantity  r (y)  is  the  optimum  relationship  between  the  degree  of 

the  root  node  and  the  expected  number  of  contending  packets  for  a tree  where 

the  root  node  is  restricted  to  be  even  and  all  others  are  binary.  It  was 

calculated  by  a procedure  analogous  to  that  used  in  Section  3.2.3  to  obtain 
* 

g^;  the  result  is. 


(3. 3. 2. 3) 


. I 1 for  y < 3.40 

r (w)  ■ I 

I r for  3.40  + (r-2)2.30  < V < 3.40  + (r-l)2.30 

Determination  of  E{6}  vs.  X 

Next  r*(v)  is  substituted  Into  Eqs.  (3. 3. 1.1)  and  (3. 3. 1.2)  to  obtain 

E{f,|y,r*(M)}  and  / E{E^ly,r*(y)}.  These  two  quantities  are  plotted  In 

2 A 

Pig.  3. 3. 2.1.  The  discontinuities  of  E{Jl  |y,r  },  which  are  evident  In  this 
figure,  are  due  to  changes  In  the  degree  of  the  root  node. 

As  In  Chapter  2,  It  can  be  shown  that 

E{ilj^j^|)lj,r*}  < 2.325X(jlj-l)  + E{£j_^j^|ilj-l,l}  (3. 3. 2. 4) 

/E{£j^j^l4^,r*}  < 2.325X(Ej-l)  + (3.3. 2. 5) 


} > 2.325XEj  + .6 


(3. 3.2. 6) 


E{J.|y,l}  for  y £ 2 

2.325(y-2)  + E{t|y-2,1}  y > 2 

Where  f^  la  a convex.  Increasing  lower  bound  to  E{l|y,r  (y)}. 

Continuing,  as  In  Chapter  2,  It  can,  further  be  shotm  that 


f/f) 


(3. 3.2. 7) 


v/  0 \ 


(3.3. 2.8) 


E{4}  > I* 
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(3. 3.2. 9) 


where 


* * 

a - ) 


and  f (A)  - f (2X1) 

c y 


^ 4.65cXE{a}  -f 
1 - (2.325A)^ 


(3.3.2.10) 


where 


c 


- 2X}  - 2.325X 


(3.3.2.11) 


and 


j^2j  ^ 2.79XE{t}  -I-  .36 
“ 1 - (2.325X)^ 


(3.3.2.12) 


Finally,  from  Eqs.  (3. 3. 2.1),  (3. 3. 7. 2)  and  (3. 3. 2. 7)  through  (3.3.2.12)  we 
have 


E{6}  < ^r?8cX 


1.05c 


1 - (2.325X)^  (1  - (2.325X)^]Jlj^(X) 


+ .321  (3.3.2.13) 


E{6}  > max 


r/l^395X ,18 (X)\;  1 (X)] 

L\  1-(2.325X)^  tl-(2.325X)^U^(X)  / J 


(3.3.2.14) 
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Equations  (3.3.2.13)  and  (3.3.2.14)  have  been  calculated  and  they  are  plotted 
In  Fig.  3. 3. 2. 2.  Note  that  the  maxlnum  arrival  rate  Is  .430  packets/slot. 

3.3.3  Average  Throughput 

Based  on  the  results  of  the  preceding  section  and  Section  2.3,  It 
follows  that  the  maximum  average  throughput  for  the  optimum  dynamic  algorithm 
Is  .430  packets/slot.  Furthermore,  the  average  delay  vs.  average  arrival 
rate  results  of  Section  3.3.2  may  be  Interpreted  as  average  delay  vs. 
average  throughput  (see  Fig.  3. 3. 2. 2). 


3.3.4  System  Stability 

In  this  section,  we  will  consider  the  stability  of  the  optimum 

* 

dynamic  tree  where  the  root  node  has  degree  2r  and  all  other  nodes  are 

binary.  The  definition  of  stability  of  the  multi-access  system  with  the 

dynamic  algorithm  Is  the  same  as  It  Is  with  the  static  algorithm.  That  Is, 

the  system  Is  k'th  order  stable  If  the  k'th  moment  of  the  delay  Is  finite. 

In  this  section  we  will  prove  that  If  E{)^}  exists,  then  the  generating 

function  of  I also  exists  (£.  corresponds  to  the  number  of  algorithm  steps). 

Since  In  Section  3.3.2,  we  showed  that  £{£}  Is  finite  Iff  X < .430  and 
k**l  k 

since  E{ [delay]  } < » If  E{E  } < <>,  the  preceding  statement  allows  ua 
to  conclude  that  all  the  moments  of  the  delay  are  finite  as  long  as  X < .430 
packets /slot. 

Let  G(8,y)  and  0^(8, u)  be  the  generating  functions  of  I given  y for  the 
binary  and  for  the  dynamic  tree  algorithms.  The  analysis  of  this  section  Is 
organized  as  follows. 
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X (packets/ slot) 


Figure  3.3. 


.2  Upper  end  Lower  Bounds  to  the  E{ delay}  versus  the 
Arrival  Rate  for  the  Optlnum  Dynamic  Tree/Poisson 
Source  System 

Note:  An  step  equals  one  round  trip  delay 

plus  two  slots 


i.  Show  that  G(s,|i)  ^ fgy  0 f.  s £ 8q»  0 5.  W £ 3.4, 

Sq  > 0,  and  e arbitrary. 

8(E,{Jlly}  + e] 

li.  Show  that  G^(8,u)  £ e for  all  y and  0 ^8  £ 8q, 

Sq  > 0.  Note  that  E^{2.|y}  la  the  Indicated  expectation  for  the 
optimum  dynamic  tree  algorithm. 

fi£. 

ill.  Argue  that  E{e  } exlats  if  E{i|y}  exlata  or  equivalently  if 
A < .430  packeta/slot. 

1.  Upperbound  to  G(a.y) 

Let 

<!„  <..u)  ■ - (^))  + y.-"  (.•-!)  + ^ 

for  0 £ 8 < ln2 

Then  in  Appendix  A2.8  we  show  that 

G(8,y)  £G^(8,y)  for  s < ln2  (3. 3. 4. 2) 


• ■] 


(3. 3.4.1) 


and 


3G 

■51 


% 

(a,y)  1 j^u^a.y)  for  a < Ini 


(3. 3.4. 3) 
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It  can  be  shown  that  If  a generating  function  Is  bounded,  then  It 

sg 

Is  continuous  In  s.  The  same  Is  true  for  Therefore,  It  follows  from 

a 

Eqs.  (3. 3. 4. 2)  and  (3. 3. 4. 3)  that  G(s,u)  and  -^^(3,^1)  are  continuous  In  s 

g 

for  s < ln2.  Note  that  G(s,y)  and  -r— G(s,m)  are  also  continuous  In  y.  This 

® k -y 

follows  from  the  fact  that  P{j  packets  are  active}  ■ — . 

g g 

Since  G(s,y)  and  i^^(s,y)  are  continuous  In  s and  y,  ■^lnG(s,y)  Is 

continuous  In  s and  y for  s < ln2  and  y.  Then,  In  any  closed  bounded  region, 
1 3 

say  0 ^ s ^ ln2,  0 £ y ^ 3.4,  -j^lnG(s,y)  Is  uniformly  continuous.  Thus 
for  any  C > 0 we  can  choose  s^  > 0 so  that 


•g^lnG(s,y)  £ -g^lnGCs.y)!  + e for  0 ^ s £ Sq,  0 ^ y ^ 3.4  (3. 3. 4. 4) 

® ® I s-0 

Therefore,  since  x— ln(G(s,y)  I ■ E{ll|y}, 

I s-0 

G(s,y)  £ for  a ^ a^,  a £ y £ 3.4  (3. 3. 4. 5) 

11.  Upperbound  to  G^(s,y) 

j 

Let  2r  (y)  be  the  degree  of  the  optimum  root  node.  Note  that  this 
* 

tree  corresponds  to  r (y)  Independent  binary  trees  where  the  mean  number 
of  contending  packets  at  each  subtree  Is  y/r*.  Therefore,  from  Eq.  (3.3.4. 5) 
we  have 

8r*[E{E|V  + e] 

Gj(8,y)  ^ e ^ ® i *0 
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I 


But  rE{J.|^}  - E^{£|p}  and £ 3.4  (see  Section  3.3.2),  therefore 

r 

8[E  {^lu)  + 3^  e] 

Gj(s,u)  < e • ® 1 ®0  (3. 3.4. 6) 

111.  Exlstance  of  E{e^^} 

In  Section  3.3.2  we  showed  that  E^{E1u}  ;i  1 + 1.163vi.  Substituting 
this  along  with  and  W“2Xf,j  Into  Eq.  (3. 3. 4. 6)  we  have 

* 11  (3.3.4.7) 

Finally,  from  Theorem  2. 4. 0.1  we  have 

E{e®^}  exists  for  X < ifa.iftK  e/3.4)  ® 1 ®o 

but  e Is  arbitrary  so  X > .430  packets/ slot . 

in&x 

sit  Ic 

It  can  be  shown  that  If  E(e”  } exists  for  s £ s^  then  so  does  E{£'‘} 

for  any  k. 

This  concludes  this  section. 
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3.4  An  Efficient  Suboptimum  Dynamic  Tree  Algorithm 


In  this  section  we  will  consider  a tree  algorithm  whose  root  node 
degree  Is  constrained  to  be  a power  of  2,  but  whose  all  other  nodes  are 
binary.  The  advantage  of  such  an  algorithm  Is  that  It  Is  relatively  easy 
to  Implement.  The  reason  for  this  Is  that  as  the  root  node  degree  is  varied 
dynamically.  It  Is  not  necessary  to  choose  a different  tree  and  a different 
addressing  scheme  each  time;  variations  In  the  root  node  degree  may  be 
realized  In  a single  binary  tree  simply  by  varying  the  depth  of  the  nodes 
where  the  algorithm  originates. 

In  Subsection  3.4.1,  we  determined  the  optimum  tree  given  the  above 
constraints.  That  Is,  If  we  let, 

(2^;  K - 1,2,...  for  1 - 0 

(3. 4. 0.1) 

2 1 > 0 

Then  we  will  show  that  K , the  K which  minimizes  E{h/vi,K},  Is  given  by 


^ ( 1 for  u ^ 3.40 

K - (3. 4. 0.2) 

I K for  3.40(2’^"^)  < y < 3.40(2*^“^),  r > 1 

In  Subsection  3.4.2,  where  the  above  algorithm  Is  analyzed,  we  obtain 
upper  and  lower  bounds  to  the  E{ delay);  these  are  Illustrated  In  Fig.  3. 4. 2.1. 
In  this  section  we  also  show  that  the  maximum  throughput  Is  less  than  .430 
but  greater  than  .420  packets/slot. 
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3.4.1  Optimizing  the  Suboptlmum  Tree 


Let  a dyanmlc  algorithm  have  a tree  with  node  degrees  given  by 


g 


1 


2^,  K ■ 1,2,3,...  for  1 ■ 0 
2 1 > 0 


(3. 4. 1.1) 


Then,  In  this  section,  we  will  determine  the  K which  minimizes  E{£./y,K}. 
Note  that  since  h ■ 22.,  minimizing  E{2/vi,K}  Is  equivalent  to  minimizing 
E{h/u,K}. 

V 

By  making  the  association  that  2r  - 2 , we  have  from  Eq.  (3. 3. 1.1) 

that 


E{2|y,K}  - E{2|-^,1} 


(3.4.1. 2) 


As  In  Section  3.2.3,  K may  be  determined  by  setting 


E{2|y,K-l}  - E{2|y,K} 


(3. 4. 1.3) 


and  then  solving  for  0(K)  for  K-  2,3,4 Where  (1(K)  Is  defined  by 


K - K for  0(K)  < y < 0(K+1) 


(3. 4. 1.4) 


This  problem  Is  considerably  simplified  If  both  sides  of  Eq.  (3.4. 1.3) 


are  divided  by  y before  It  Is  solved.  (Note  that  this  operation  does  not 
affect  the  solution  0(k).)  The  simplification  arises  from  the  fact  that 


K K 

E{£|y,2  }/p  Is  simply  a function  of  u/2  , therefore,  as  will  be  shown  below, 
0(K)  ■ 0(2)2^  and  It  follows  that  Eq.  (3. 4. 1.3)  need  be  solved  only  for 
K - 2.  Let 


H(y/2‘^)  S E{Jl|^,l} 
then  If  0(2)  satisfies 


H(y/2)  - H(u/4) 


(3. 4. 1.5) 


(3. 4. 1.6) 


It  follows  that 

0(K)  - 0(2)2*^’^ 


(3. 4. 1.7) 


will  satisfy 

H(y/2^"^)  - H(y/2^) 


(3. 4. 1.8) 


Equation  (3. 4. 1.6)  was  solved  numerically;  and  the  answer  Is 

y(2)  - 3.40  (3. 4. 1.9) 

The  final  result  to  this  section  follows  fron  Eqs.  (3. 4. 1.4),  (3. 4. 1.7)  and 
(3. 4. 1.9).  It  Is  given  by 
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1 for  u < 3.40 


(3.4.1.10) 


* 

K 


K for  3.40(2^"^)  < u < 3.40(2^“^),  K > 1 


3.4.2  Analysis  of  the  Subopt Imum  Algorithm 

The  first  step  of  the  analysis  Is  to  determine  0 and  0,,  the  slopes 

U X 

of  the  straight  lines  that  respectively  upper  and  lower  bound  E{£/y,R  }. 
These  two  quantities  are  given  by, 


max 

y>3.4 


E{Jl/y.K  } 

y 


(3. 4. 2.1) 


0,  - min  ^ (3. 4. 2. 2) 

^ y>3.4  ^ 

From  Eq.  (3.4.1.10)  and  the  fact  that  E{il|y,K}/y  Is  simply  a function  of 

j( 

y/2  , It  can  be  shown  that 


max 

y>3.4 


E{l|y.K*> 

y 


max 

3.4£y<6.8 


E{E|y.K*} 

y 


and 


min 

y^3.4 


E{t|y.K*} , 

y 


min 

3.43i<6.8 


E{&|y.K*} 

y 


(3. 4. 2. 3) 


(3. 4. 2. 4) 


-172- 


Equations  (3.4. 2. 3)  and  (3. 4. 2. 4)  were  solved  numerically  and  the  answers 
are 


6 - 1.189 
u 


(3. 4. 2. 5) 


- 1.164 


(3. 4. 2. 6) 


The  corresponding  B's  of  equal  those  of  E{il|y}.  To 

ir 

this,  first  note  that  by  setting  2r  > 2 In  Eq.  (3. 3. 1.2)  we  have 


see 


* * * * * 

E{i^|y,K*}  - 2‘^"^(2^"^-l)E^{«.|y/2‘^“^,l}  + 2‘^"^E{&^ |y/2^"^,l}  (3. 4. 2. 7) 


Next  take  the  square  root  of  this  expression,  divide  by  y,  let  y 
and  use  Eq.  (3.4.1.10)  to  obtain. 


/ E{£^|y.K^}  _ E{&|y.K*} 


< 1.189 
> 1.164 


(3. 4. 2. 8) 


Using  the  above  results  and  numerical  calculation  of  E{2.|y,K  } and 

, 2i  *1 

EIE  |y,K  } It  can  be  shown  as  In  Section  3.3.2  that 

E{£}  11+  1 i 2.37SX  (3. 4. 2. 9) 

E{a}  > I*  (3.4.2.10) 

where  £*  ■ f (I*) 
c 
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(U)  > f ,(2A)l) 
c U 


f 


y 


E{iL|  u,K-l}  for  y £ 2 

2.328(y-2)  + E{Jl|u-2,K-l}  y > 2 


Furthermore  It  can  be  shown  that 


g{j^2j  ^ 4.756cAE{Jl}  + 

1 - (2.378A)^ 


(3.4.2.11) 


where 


c 


/ E{ilj^j^|y-2A}  - 


2.378A 


and 


j;{£2j  ^ 2.328E{11}  + .25 
“ 1 - (2.328A)^ 


(3.4.2.12) 


Now  substitute  Eqs.  (3.4.2.9-3.4.2.12)  into  Eqs.  (3. 3. 2.1)  and  (3. 3. 2. 2)  to 
obtain 


E{6}  ^ ^ ^ 

1 - (2.378A)^ 


+ 


l.OSc^ 

1-  j:2.378A)5^(A) 


+ .321 


(3.4.2.13) 
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and 


E{6}  > 


'/I- 164  ^ 

.\1  - (2. 


.125 


328X)' 


1 - (2.328X)^  2^(1) 


+ i H^Cl) 


) = ‘t»>] 

(3.4.2.14) 


Equations  (3.4.2.13)  and  (3.4.2.14)  are  Che  desired  bounds;  they  are 
Illustrated  In  Fig.  3.4. 2.1.  Note  that  Che  upper  bound  approaches  Infinity 
at  X ■ .420,  whereas,  the  lower  bound  approaches  Infinity  at  X ■ .430. 
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APPENDIX  A3 


SELECTED  THEOREMS  OF  CHAPTER  3 


Theorem  A3.1;  Let  .1  < y < 8;  — = — < 1.2; 

- 


OD 

L(gi.g2.y)  - giC(y)  + gigjCC^)  + 2 I 2Sig25(-^) 


^ 1-0  *1*2^ 


and 


A,  U ^2^ 


Chen; 


L(8^fg2fy)  ^ gg2»y^ 


Proof:  This  follows  from  Property  ID  In  Appendix  A. 2. 2. 

Theorem  A3. 2:  Let  L2^(g^,g2tU}  be  as  In  Eq.  (A. 3. 2),  furthermore 
gl  ^ 10,  ^2  ^ ^ .1  £ y £ 8.  Then  Che  minimum  of  L2(g]^,g2gy) 
at  g^  ■ 10,  g2  - 3. 


Proof : 

aij 


This  theorem  will  be  proved  by  showing  Chat 

> 0.  Differentiating  Eq.  (A3. 2), 

■3 


^^1 

•» — > 0 and 


(A. 3.1) 


(A. 3. 2) 


(A.  3. 3) 

QED 

let 

occurs 
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(A3. 4) 


3L„ 


agj  " 


iC(^)-2(^)' 

*1  *1*2 


since  g.  >10,  — > 0 (A3. 5) 

~ 8182 


^•|(^)^  - Property  IB  in 

*1  ^ *1  *1*2  Appendix  A2.2 


(A3. 6) 


> — (^)^  - — since  a > 3 

- 18 'g/  18  *2  - ■* 


(A3. 7) 


> 0 


since  ^ < T 

gl-6 


(A3. 8) 


For  Che  second  part  of  this  proof  substitute  Che  expression  for  C(u) 
given  in  Property  lA  in  Appendix  A2.2  into  Eq.  (A3. 2)  to  obtain. 


3L, 


3g.  " ^^*2  ^ (k-1)! 

^ k-1  ^ *’ 


k-1 


^ (-^)^  + 2y(-^^^^^ — + — — 

'i-a-1 ' 2 2 9 3 2^ 

8^  82^82 


c(p)  - ^ (^)^  - 

^ *1  ■*  *1  *2  *1  9g2  *1  *1 


C(u)  - .t  .j  - 3 


(A3. 9) 
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5(u)  > 5(1.2)  - .337 


But 


^ < .312 

gj^  O gj^  81  — 


Therefore  from  Eq.  (A3. 9) 


3L- 

^>0  for  lJ>1.2 

*»1 


Eqs.  (A3. 11)  end  (A3. 12)  taken  together  conclude  the  proof. 


(A3. 12) 


QED 
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Theorem  A3. 3;  Let  Lj^(g^,g2,u)  be  as  given  in  Eq.  (A3. 2);  furthermore, 
let  2 £ gj^  £ 9,  ^2  — ^ .1  £ y £ 8.  Then  the  minimum  of  Lj^(gj^,g2,y) 

occurs  at  g2  ■■  6. 


Proof:  We  will  prove  this  by  shoving  that  when 


3Lj^ 

in  the  above  statement  then  -r— 


> 0. 


From  Eq, 


g^,g2,y  are  as  specified 
(A3. 5)  we  have  that 


3La 


"82  ®1®2 


i *2  i ‘ 


(f-)^I5(%^  5(¥-)  - tIi 


g 


y 


gl'  18^ 


Since  ^ > 0,  we  need  to  show  that  (7^)^^(^)  > t4. 
gj^  — y gj^  — 18 


Now  since  2 £ gj^  ^ 9 and  .1  £ y ^ 8 it  follows  that 


90  - gj  - “ 


d 2 

And  since  ^(^(x)/x  ) £ 0,  we  have  that 


> (T)^e(y)  for  0 < ^ < 4 


'gj  - '4 


*1 


- .05678 
> 1/18 


(A3. 13) 

(A3. 14) 


(A3. 15) 


(A3. 16) 

(^.17) 
(A3. 18) 
QED 
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Theorem  A3. 4;  Define, 


■ 5 + 2 I ^ 

1^  2^8 

And  let  w ^ 8 and  g be  Integer  and  greater  than  one.  Then 

min  H^(g,u)  £ 2.34 
8 

Proof;  This  will  be  proved  by  first  finding  an  Interval  X such  that 

U ^ 8,  an  Integer  y > exists  such  that  y/g  e X and  then  showing  that 

^u^a^  — 2-34.  First,  we  prove  that  X - ; 1.1  < < 1.257} 

(y/g)eX  * g - g - 

satisfies  the  above  criteria.  It  Is  sufficient  to  show  that 

sup  C min  (^)']  - 1.257 
W>8|_g3j/l.l 

But  m|n  ^ ^ for  8 < u < 8(1.1).  Therefore, 

sup  [#]  - (1.1)|  - 1.257 

8£m<8.8  ' 


(A3. 19) 

(A3. 20) 

for  any 

(A3. 21) 

(A3. 22) 
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similarly  for  1.1k  £ P ^ l.l(fcfl)  ; y ^ 8 

EquaClons  (A3. 22)  and  (A3. 23)  prove  Eq.  (A3. 21) 
Now  leC 


X - y/g 


(A3, 24) 


then  from  Eq.  (A3. 19)  we  have 


H (x)  - - + 

u X 


I 1 5(7) 
1-0 


(A3. 25) 


Finally,  H^(x)  was  maximized  over  xcX  and  the  result  Is 


max  H (x)  - 2.34  (A3. 26) 

xeX  “ 


QED 


Theorem  A3. 5;  Define 


Hj^(x,g) 


^(8)  + g^  + 2 


I 

1-0 


^(^) 

2^g 


(A3. 27) 
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and  let  0 < X £ .5  and  g ^ 3.  Then 

Hj^(x.g)  ^ 2.498  (A3. 28) 

Proof:  Since  the  sum  in  Eq.  (A3. 27)  is  positive  we  have  that 
H^(x,g)  1^(8)  + 

^ ^(8)  + 3^^-^  for  g > 3 

l^(8)+|x-x^  (A3. 31) 

The  last  equation  follows  from  Property  IB  of  Appendix  A2.1.2.  Now  for 
0 X £ 'S*  the  right  side  of  Eq.  (A3. 31)  la  decreasing  in  x.  Therefore, 
evaluating  it  at  x - .5  we  have 

H^(x,g)  > 2.498 

QED 

Theorem  A3.6t  Let  Hj^(x,g)  be  as  in  Eq.  (A3. 27).  Also  let  .5  < x < 3 
and  g ^ 3.  Then 

Hj^(x,g)  ^ 2.34  (A3. 32) 
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I 


(A3. 29) 


(A3. 30) 


I 


Proof:  First  we  prove  this  for  g > 12.  It  can  easily  be  shown  that. 


. (1  + X + x^)e'*]  _<  0 for  X £ 3 

dx*^*-  * -I  x"* 


therefore 


min  ] 

.5<x<3  * 


min  [ 


iLii . mil 

.5  ’ 3 ^ 


mil 

.5 


(A3. 33) 


But  from  Eq.  (A3. 29) 


Hj^(x,g)  > 3C(8)  + 12^^!^  for  .5  £ X £ 3 and  g > 12 


(A3. 34) 


- 2.497  for  g > 12  (A3. 35) 

Finally,  Hj^(x,g)  la  minimized  over  g ■ 3, 4,..., 12  and  .5  ^ x ^ 3.  The 
result  Is 


Hj^(x,g)  ^ 2.34  for  *5  £ x £ 3 and  g • 2,3, ...  ,12  (A3. 36) 


Theorem  A3. 7:  Define  H(g^,g2,y)  as 


81^2  u ? 

H(g^,g2,M)  - ^(U)  + + 2^2 

^ 1-0 


.1  *1*2, 


(A3. 37) 


y 8182 
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and  let  g^,g2  be  Integer  and  greater  than  one.  Furthermore,  let  y ^ 8 

and  ^ > 3. 

*1 

then  the  min  H(g. ,g,,y)  lies  on  g-  - 2 (A3. 38) 

0-1  O*}  i A Z 


Proof:  Since  ^(y)  is  an  increasing  function  we  have  for  y > 8 that 


i-0  *1*2 


(A3. 39) 


H(g 


8i  I,  II 

,.g,.u)  > ^tKCS)  + * 2 I r)  (A3. 40) 

1 1.0  *1*2^ 


“ 2 g.g. 


for  y > 8 


y/gi 


Now  for  2 ^ g2  ^ ■ the  following  holds 


5(8)  + g25(y/gi)  1 82 


(A3. 41) 


therefore. 


g^go  ® 2 g'lg^  YJ 

H(8i.g2.^')i^+  2 I — ^(_L) 


(A3. 42) 


i-O 


*1*2^ 


and 
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(A3. 43) 


ao  ^ 

H(g,,g2,u)  ^ mln[^  + 2 


But  the  right  side  of  Eq.  (A3. 43)  Is  the  upper  bound  (Eq.  (A3. 39))  evaluated 

£(8)e^^®l 

at  g2  * 2 and  minimized  over  gj^.  So  we  have  proved  that  if  g2  ^ — , 

then  the  minimum  of  H(gj^,g2,y)  lies  on  g2  ■ 2. 

The  proof  is  concluded  by  demonstrating  that  Eq.  (A3. 43)  also  holds 
r(8)e^^*l 

if  g2  > + y/g  • Letting  x ■ y/g  in  Eq.  (A3. 40)  we  have  from  Theorem  A3.1 


H(x,g2)  > ^(8)  + - I (x/g2)^ 


(A3. 44) 


> ^(8)  4.  ^(x)  + g - f(ife  e-)2  for  g2  > ^ (A3.45) 


Now  let  y be  defined  by 


y - g 


C(x)  ^ 

2 X g- 


(A3.46) 


Then  it  can  easily  be  shown  that  y is  convex  in  g2  for  x ^ 0 with  the 
minimum  occurring  at 


‘VtliT 


*2  " C(x) 


(A3. 47) 


-186- 


4 


But  for  X > 3 


* ^ ® 

82 -1?;^ 


(A3. 48) 


and  since  y Is  convex  It  follows  that 


y > fr-  + 2x(~)  for  X > 3 

^ — 1+X  X X — 

e 


Therefore  from  the  above  equation  and  Eq.  (A3. 45)  we  have  that 


’®2  — X 1+x  X X 

e 


8/  X Xn2 

9'T5r®  > 


> 2.55 


and  from  Theorem  A3. 4 we  conclude  that 


^ min  [J+2  \ ^:^(-^)] 

8l"2  i.o  " 2S 


QED 
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CHAPTER  4 


THE  TREE  ALGORITHMS  WITH  A FINITE  SOURCE  MODEL 
4.1  Introduction 

In  this  chapter  we  will  consider  a multiple  access  system  with  2^ 
Independent  sources.  The  transmissions  process  Is  essentially  the  same  as 
that  which  Is  described  In  Chapter  2.  That  Is,  packets  that  arrive  In  one 
epoch  are  processed  In  the  following  epoch  by  a tree  algorithm  (both  static 
and  dynamic  algorithms  will  be  considered  In  this  chapter).  Here,  as  In 
Chapter  2,  a group  of  sources  chat  has  undergone  a conflict  Is  divided  In 
half  and  the  two  subgroups  transmit  their  packets  In  two  consecutive  slots. 

The  decision  to  divide  a group  In  half,  the  transmission  of  the  two  halves 
and  the  observation  of  Che  results  of  those  two  transmissions  constitute 
an  algorithmic  step.  The  number  of  algorithmic  steps  In  one  epoch  will  be 
designated  by  H. 

It  will  be  assumed  chat  a source  may  receive  at  most  one  new  packet  per 
epoch.  The  probability  Chat  a source  will  receive  a packet  In  Che  next  step, 
given  Chat  It  has  not  yet  received  one  In  the  present  epoch.  Is  a constant  and 
will  be  designated  by  p.  It  follows  Chen,  that  q,  Che  probability  that  a 
source  will  receive  a packet  In  an  epoch  of  length  H,  Is  given  by 

q - 1 - (l-p)*".  (4. 1.0.1) 

It  should  be  noted  Chat  In  Che  above  source  model,  a source  can  have 
at  most  two  packets  at  any  one  time;  one  Chat  arrived  In  the  previous  epoch 
and  which  Is  In  Che  process  of  being  transmitted,  and  one  Chat  arrived  In 
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the  present  epoch  and  which  will  be  processed  In  the  following  epoch. 

The  channel  model  la  the  same  as  that  which  was  considered  In  Chap- 
ters 2 and  3.  That  Is,  It  Is  slotted  and  the  sources  have  Che  means 
to  determine  whether  there  are  0,  1 or  more  chan  1 packets  In  any  one  slot. 
If  a slot  contains  more  than  one  packet,  then  It  la  assumed  that  no  one 
gets  through. 

The  tree  searches  of  the  algorithms  that  we  will  consider  here  will 
be  carried  out  serially  and  Che  source  addresses  are  assigned  deterministi- 
cally. The  parallel  search  and  the  random  address  assignment  will  not  be 
considered  here.  It  Is  Important  to  point  out  however  that,  even  chough 
for  Che  Poisson  source  model  the  random  and  deterministic  source  assign- 
ments had  the  same  delay  and  throughput  properties.  In  the  finite  source 
model  It  can  be  shown  that  the  random  address  assignment  has  larger  average 
delay  chan  that  of  the  deterministic. 

The  rest  of  this  chapter  is  organized  Into  two  sections.  In  Section 
4.2  we  consider  the  static  binary  tree.  Here  we  obtain  an  upper  bound  to 
Che  average  delay  and  a lower  bound  to  the  average  throughput  In  tern  of 
p.  We  also  combine  these  two  bounds  to  obtain  an  E{ delay}  vs  E{ throughput) 
performance  curve.  This  Is  Illustrated  in  Fig.  4.2. 2.2  for  M*6.  In  Sec- 
tion 4.3  we  consider  the  optimum  dynasdc  tree  algorithm.  Here  we  restrict 
all  nodes  to  be  binary  except  for  the  root  node  which  Is  allowed  to  have 

II 

a degree  chat  Is  a power  of  2 but  less  than  or  equal  to  2 . Subject  to 
the  preceding  restrictions,  g^  Is  chosen  so  as  to  minimize  the  expected  num- 
ber of  slots  needed  to  process  the  contending  packets,  given  q.  The  optimum 
gg  Is  given  by  Eqs.  (4. 3. 1.1)  end  (4. 3. 1.6).  Following  the  determination 
of  the  optimum  tree,  first,  wa  obtain  upper  and  lower  bounds  to  the  average 
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delay  and  average  throughput  respectively,  and  then  obtain  an  E{ delay}  vs 
E{ throughput)  curve.  The  E{ delay)  vs  E{ throughput)  curve  Is  shown  In 
A. 3. 2. 5.  Section  4.3  Is  concluded  with  a theorem  proving  that  the 
E{delay}  for  the  optimum  dynamic  tree  protocol  Is  less  than  or  equal  to  the 
E{ delay)  of  the  TOMA  protocol. 
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I 

i 

^•2  Static  Biiurv  Tre«  Algorithm  with  Finite  Source  Modal 

section  contains  two  subsections:  4.2.1  where  an  upper  bound  to 
the  E{ delay}  Is  obtained  and  4.2.2  where  a lower  bound  to  E{ throughput)  Is 
obtained.  Both  of  these  bounds  are  functions  of  p,  the  traffic  parameter. 

They  have  been  computed  for  1^6  and  they  are  plotted  In  Elga.  4. 2. 1.4  and 

4. 2. 2.1.  The  E{delay}  va.  E{ throughput)  curve  la  given  In  4. 2. 2.1. 

4.2.1  Average  Delay 

We  begin  this  section  by  presenting  several  definitions.  We  will  be 
using  the  following  quantities:  e^,  1^,  d^  6,  p,  and  q.  The  first  five 
of  these  are  defined  In  Section  2.2  and  the  last  2 are  defined  in  Section 

4.1.  Furthermore,  let 

'l>(q.n)  - 1 - (l-q)2“  - 2*q(l-q)^2““^>  (4.2.1.1) 

(4.2. 1.2) 

(4. 2.1.3) 

t|>(q,a)  is  the  probability  that  there  are  at  least  two  active  sources  In  a 
branch  with  2*  leaves.  4(q»m)  Is  the  probability  that  there  Is  at  least 
one  other  active  source  la  a branch  with  2*  leaves  given  that  one  particular 
source  of  that  branch  Is  active.  D(q,m)  Is  the  expected  nuid>er  of  nodes 
visited  In  a branch  with  2*  leaves.  These  quantltltes  are  considered 
In  more  detail  In  Appendix  A4. 

The  analysla  Is  carried  out  aa  follows:  | 

1.  Calculate  E{d|lj^,p}.  j 
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♦(qfW)  - 1 - (1-q)^^ 

1W"1  . 

D(q,n)  - ^ 2^|f(q,m-l) 

1-0 


I 


II.  Show  that  E{6|£^,p}  la  concava  and  Incraaalng  In 

III.  Calculate  1^,  an  upper  bound  to  £{!}. 

Iv.  Apply  Jenssen's  Inequality  to  the  above  three  steps  to  prove 
that  E{6lp}  < E{6|I^,p}  5 6^. 

1.  Derivation  of  E{d|Ej^,p} 


The  delay  that  a packet  uxtdergoes  can  be  decomposed  Into  d^  and 
62  where  d^  Is  the  time  spent  by  the  packet  In  the  epoch  of  arrival  and 
d2  Is  the  time  spent  In  the  following  epoch,  l.e. , the  epoch  where  It  Is 
successfully  transmitted.  As  a consequence  of  this  observation  we  have 
that 

R{6|l^,p}  - E{d^|4j^,p}  + Efd^lij^.p}  (4.2. 1.4) 

The  subscript  of  I will  not  be  used  where  aad>lgultles  do  not  arise, 
note  that  refers  to  length  of  the  epoch  In  which  the  packet  arrived 
and  whan  a random  variable  Is  conditioned  on  1 we  will  mean  1^^. 

Next  an  expression  for  E{dj^|i,p}  will  be  determined;  this  Is  accom- 
plished by  noting  that,  given  that  has  length  1,  the  delay  of  a packet 
arriving  In  the  step  of  la, 

d^  - 1 - J - 1 (4.2.1.5) 


The  probability  of  a given  packet  arriving  In  the  step  Is 

o 
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P(jll,p) 


for 


A-1. 


(4. 2. 1.6) 


1 - (1-P) 


Jfaltlplylng  Eqs.  (4.2. 1.5)  «nd  (4.2.1. 6)  and  sunning  over  j we  have 

E{dJl.p}  - 1 - 1 - 1-P  - (4.2.1.7) 

^ p(l  - (l-p)**) 

This  la  the  desired  expression  for  E{dj^(l,p}.  E{d2|ltp}  Is  derived  In 
Appendix  A4  and  It  Is  rewritten  below, 

N-1 

E{dJl,p}-l+  I [1  +4  D(q,»-1)]  ♦(q.H-l)  (4.2.1.8) 

^ 1-1  ^ 

where  q,  and  D are  given  by  Eqs.  (4. 1.0.1),  (4. 2. 1.2)  and  (4. 2. 1.3).  The 
•xpresslon  for  E{5|l,p}  follows  fron  Eqs.  (4. 2. 1.4),  (4. 2. 1.7)  and  (4. 2. 1.8). 

II.  Properties  of  ECdIt.p} 

For  M-6  and  selected  values  of  p.  the  quantity  E{8|i,p}  was  cal- 
culated for  1-1, 2, 3,... 63.  The  results  are  shown  la  Fig.  4. 2. 1.1.  From  the 
cowputer  printout,  as  well  as  fron  this  figure.  It  la  evident  that  E{6|l,p} 
Is  concave  and  Increasing  in  1.  These  properties  will  be  uaed  In  Subsection 
Iv. 

III.  Petemlnatlon  of  Upper  Bound  to  E{t} 

In  this  subsection,  first,  ws  will  derive  7^,  an  upper  bound  to 
E{1},  and  then  compute  that  upper  bound  for  1h6.  The  upper  bound  follows 
from  the  followlag  theorem. 
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Theorem  4. 2.1. It  Let  ij  be  Che  outcome  of  a Markov  chain  after  the 
CranaiCion  and  let  f^(£j}  be  a soadecreaalng  concave  upper  bound  Co 
Furthermore,  let  I*  be  defined  by 

(4. 2. 1.9) 

and  asaume  that  4^  > fy(Ej)  for  Ij  > 4*  then 

EUj}  < 4*  = 4^  (4.2.1.10) 

Proof:  This  proof  ia  similar  to  that  of  Theorem  2. 2. 3. 2.  Therefore,  it 
will  not  be  given. 

QED 

Next  we  will  calculate  for  »«6.  In  Appendix  A4  the  following  ex- 
pression is  derived. 

E{42|q}  - 1 + 2D(q,N-l)  (4.2.1.11) 

Equation  (4.2.1.11)  is  plotted  in  Fig.  4. 2. 1.2  for  M*6,  where  for  comparison 
we  also  plot  on  the  same  figure  E{42|q}  for  the  optimum  dynamic  tree  — 
thla  ia  derived  in  Section  4.3.1.  E{42|42^,p}  follows  from  Eqs.  (4.2.1.11) 

and  (4. 1.0.1).  This  quantity  la  llluatrated  in  Fig.  4. 2. 1.3.  Aa  can  be 
seen  from  this  figure,  for  p > .016,  f^(4*)  - E{42|4*}.  Equation  (4. 2. 1.9) 
was  solved  for  this  f^  and  for  p ^ .016.  The  resulta  are  listed  In  Table 
^•2.1.1.  For  p ■ .004,  .008  and  .012,  the  following  three  linear  expressions 
upperbound  E(42|42}  rsspeetivsly. 
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l 

r 


E (alg>’.  steps) 


p 


I 5 th^  th. 

u u 1 2 


004 

1.04 

1.24 

.125 

001 

1.2 

1.6 

.232 

012 

1.9 

2.6 

.314 

.016 

19 

22.5 

.372 

.331 

.020 

32 

36 

.410 

.370 

.030 

50 

39.5 

.436 

.436 

.040 

57 

70 

.436 

.470 

050 

60 

77 

.437 

.489 

.100 

63 

89 

.439 

.507 

TABLE  4.2.1. 

.1 

AVERAGE  DELAY  AND  AVERAGE  THROUGHPUT 

FOR 

BINARY  TREE  ALGORITai  WITH  64  SOURCES 
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E{i2lii.-004} 

<,  .2241j^  + .804 

(4.2.1.12) 

Ed^l^i.^OOS} 

< .500£j^  + .605 

(4.2.1.13) 

ECijUi.-OW} 

< .7501 + .473  . 

(4.2.1.14) 

These  f^'s  were  used  in  Eq.  (4. 2. 1.9)  to  obtain  for  the  Indicated 
p's.  The  results  are  listed  In  Table  4. 2. 1.1.  On  that  same  table  we  also 
present  5^,  th^  and  tb^.  5^  Is  the  upper  bound  to  the  E{ delay}  which  Is  de- 
termined In  the  following  subsection  and  th^  and  th2  are  lower  bounds  to 
the  expected  throughput  which  are  derived  In  Section  4.2.2. 

Iv.  Determination  of  Upper  Bound  to  E{ delay} 

In  Subsection  11,  It  was  shown  E{5|£,p}  Is  concave  and  nondecreaa- 
Ing  In  1.  Therefore,  from  Jenssen's  Inequality  and  the  concavity  of 
E{d|l,p}  It  follows  that 

E{«lp}  < E{6ll,p}  (4.2.1.15) 

and  from  the  nondecreasing  property  of  E{5|E,p}  we  have  the  following  upper 
bound. 

E{6|p}  < E{«|Tj^,p}  5 5^  (4.2.1.16) 

Dote  that  the  right  side  of  the  above  equation  equals  Eq.  (4. 2. 1.4)  with 

i 
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replaced  by  The  valusa  of  A^(p)  that  were  calculated  in  the  pre- 

ceding subsection  were  substituted  into  Eq.  (4.2.1.16)  and  the  results 
are  presented  in  Table  4. 2. 1.1  and  in  Fig.  4. 2. 1.4. 

One  more  result  will  be  developed  before  comcluding  this  section;  it 
is  an  exact  expression  for  £{ delay}  at  p»l.  This  is  a useful  quantity 
since  it  is  the  maximum  E{ delay}  over  p.  When  p"l,  a packet  arrives  at 
each  source  in  the  first  step  of  each  epoch  with  probability  one.  There- 
fore, it  can  be  shown  that 

E{dj^|p-1}  - 2^-2  (4.2.1.17) 

and  E{d^|p-1}  (q. 2 1.18) 

and  we  have 

E{6|p-1}  - ~ ■ (4.2.1.19) 

This  should  be  compared  with  the  maximum  average  delay  for  the  optimum 
dynamic  tree  which,  as  we  will  see,  is  given  by 

Edy„{5lp-1}  - (3(2“‘h-l)/2  (4.2.1.20) 

This  concludes  Section  4.2.1.  Next  we  consider  the  expected  throughput. 
4.2.2  Average  Throughput 

This  section  has  two  objectives.  The  first  is  the  determination  of 
a lower  bound  on  the  E{ throughput)  as  a function  of  p,  and  the  second  is 
the  determination  of  an  E{ delay)  vs  e{ throughput)  curve.  These  two  results 
are  displayed  in  Figs.  4. 2. 2.1  and  4. 2. 2. 2. 
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The  definition  of  throughput  for  the  finite  source  model  Is  the  same 
as  It  was  for  the  Poisson  source  model.  That  Is,  the  average  throughput 


is  the  fraction  of  time  that  the  channel  contains  valid  data,  l.e.,  exactly 
one  packst/slot.  Even  though  the  definition  of  E( throughput)  Is  the  same 
for  both  cases,  the  computations  are  more  Involved  In  the  finite  source 
model  because  of  the  assumptions  that  a source  can  accept  only  one  packet 
par  epoch. 

We  begin  by  deriving  an  expression  for  E{ throughput)  In  terms  of  the 
system  parameters.  By  using  the  law  of  large  numbers,  one  sees  that 
In  an  interval  of  length  2kE{ll)  slots  the  number  of  packets  that  Is  succes- 
sfully  transmitted  approaches  k[2  E{q)  + o(k)]  for  large  k,  where  1,  q, 
and  N are  as  defined  previously  and  o(k)*0.  It  follows  then  that. 


E{ throughput) 


E{q)2” 

2EU) 


(4.2. 2.1) 


Mote  that  In  the  above  equation  both  E{q)  and  E{E)  are  functions  of  the 
traffic  parameter  p. 

Next  we  will  develop  two  lower  boxinds  to  E{ throughput) ; they  will  be 
designated  by  thj^(p)  and  th2(p).  Since,  as  will  be  shown  shortly,  neither 
of  these  bounds  Is  tightest  over  all  p,  we  will  take  the  lower  bound  to  the 
E{ throughput)  to  be. 


th^  ■ Max{th2,  th2) 


(4. 2.2. 2) 


Nov  we  will  derive  th^.  Equation  (4. 2. 2.1)  can  be  rewritten  as  follows 
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E{ throughput}  - 2 


N-1  I 


. .»-i  1.  "’IV  -p 
i ^ ij  Ba'jiy  ■ “i 


(4. 2. 2. 3) 


Vlhere  E{q|ij^}  and  are  given  by  Eqa.  (4.1.0.1)  and  (4.2.1.11). 

Next  we  will  derive  th^.  Flrat  note  that  E{ilq}  (see  Eq.  4. 1.0.1) 
can  be  lower  bounded  as  follows: 


E{q|E}  > ai  - a + p for  1-1,2, ... ,2*^-1 


where 


2»-i 

1 - (1-p)^  p 

2«-2 


Taking  expectations  of  both  aides  of  Eq.  (4. 2. 2. 4)  we  have 


E{q}  ^ a E{4}  - a + p 


Substituting  this  Into  Eq.  (4. 2. 2.1)  we  have 


E{th}  > (•  + 


It  can  be  shown  that  p - a > 0,  and  therefore 
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(4. 2. 2. 4) 


(4.2. 2.5) 


(4. 2. 2. 6) 


(4.2.2. 7) 


4 


(4. 2. 2. 8) 


E{th}  > [a  + ^^]2^  ^ = th, 
u 

Thla  is  tha  daalred  fora  for  th2*  ^ calculated  for 

N*6  for  sevaral  values  of  p.  The  results  are  presented  in  Table  4.2.1«1 
and  In  71e«  4«2«2«1« 

The  performance  curve  (fi  vs  th, ) can  be  obtained  directly  from 

u *> 

Table  4. 2. 1.1.  This  Is  plotted  In  Fig.  4. 2. 2. 2. 
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ur^at 


thj 

fh, 


p 


Flgur*  4. 2. 2.1  Lo«*r  Bounds  to  ths  Avsrags  Throughput  versus  p for 
the  Binary  Tree  Algorltha  with  64  Sources 


I 


1 

* 


E[thr]  (packets/slot) 

Figure  4. 2.2.2  The  Upper  Bound  to  E{deley}  versue  the  Lower  Bound 
the  B{ throughput)  for  the  Binary  Tree  Algoritha 
with  64  Sourcea 

Note:  An  Algr.  etep  equele  one  round  trip  delay 
plue  two  alota 


4.3  Optimum  Dynamic  Tree  Algorithm  with  Finite  Source  Model 

There  ere  two  major  objectives  Co  this  section;  the  determination  of 
Che  optimum  finite  tree  algorithm  and  the  analysis  of  that  algorithm.  These 
two  problems  are  considered  In  Subsection  4.3.1  and  4.3.2,  respectively.  In 
Subsection  4.3.2  we  also  prove  that  the  optimum  d3maad.c  tree  algorithm  Is 
superior  to  the  TDMA  protocol. 

4.3.1  Optimum  Tree 

The  criterion  of  optimality  to  be  used  here  Is  the  same  as  that  of 

Chapter  3.  That  is,  the  optimum  tree  Is  that  which  minimizes  the  expected 

number  of  slots  needed  to  process  the  2 sources  given  that  the  probability 

that  any  one  of  them  has  a packet  to  transmit  Is  q.  Since  the  number  of 

sources  Is  finite,  the  optimization  Is  carried  out  over  a smaller  set 

of  trees  then  It  was  In  the  Poisson  source  model.  More  specifically,  we 

are  going  to  restrict  all  nodes  to  be  binary  except  for  the  root  node  which 

r 

can  have  a degree  that  Is  a power  of  2.  In  other  words,  if  we  let  g^"!  , 
then  the  problem  Is  to  choose  X-1,2,...,M  so  that  2E{)l|q,M,K}  Is  mlnisnim. 

The  reason  for  restricting  the  initial  degree  to  be  a power  of  two  Is 
easier  Implementation.  Variations  in  the  degree  of  the  root  node  under  this 
restriction  are  equivalent  to  starting  the  binary  tree  algorithm  at  different 
levels,  thus  not  requiring  a different  tree  each  time  the  degree  of  the  root 

n 

node  Is  changed.  For  example,  an  algorithm  whose  tree  has  2 leaves  and 
root  node  degree  eqtial  to  2^  Is  equivalent  to  a binary  tree  algorithm  that 
starts  with  the  nodes  of  depth  K. 

Mow  we  will  determine  K*,  the  optimum  K,  as  a function  of  q.  More 

A 

specifically  we  will  obtain  an  equation  whose  solution  is  q(K,N),  where 

A 

q(K,N)  Is  defined  by 


i 
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(4. 3. 1.1) 


K*  - K for  q(K,H)  < q £ q(K+l,N) 

A A 

As  ve  will  shortly  sss,  q(K,II)  ■ q(»-K).  It  can  be  shown  by  a procedure 

n 

analogous  to  that  of  Section  3.3.1  that  for  a tree  with  2 leaves  and  a 
root  node  degree  2^,  that 

EdjIq.NfK)  - 2*^"^  E{l2hiN-K+l.l}  (4. 3. 1.2) 

Substituting  Eq.  (4.2.1.11)  Into  this  expression  results  In, 

EdjU.H.K}  - 2^"^  {1  + 2D(q,N-K)}  (4. 3. 1.3) 

A 

q(N,K)  Is  deteralned  by  setting 

EdjIq.H.K-l)  - B{)t2|q,!l,K}  (4. 3.1.4) 

and  then  solving  for  q.  Substituting  Eqs.  (4. 2. 1.3)  and  (4. 3. 1.3)  Into 
Eq.  (4. 3. 1.4)  and  then  rearranging  we  have  the  following  expression  which 

A 

defines  q. 


H-E+1)  - 1/2 


or 


11  + q(2**"‘^‘*’^^-l)l(l-q)^^"'*^  ^-1)  • 1/2 


(4.3.1. 5) 
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Equation  (4. 3.1.5)  was  solved  for  N-K*^,l,2,3t4,5.  The  results  are 


given  In  Table  4. 3. 1.1  below: 


.707 

.38 

.20 

.10 

.05 

.025 


H-K 


0 

1 

2 

3 

4 

5 


TABLE  4.3.1.1 
OPTIMUM  SOOT  NODE  DEGREE 


Froa  the  above  table  the  following  expression  for  q(M-R)  Is  evident 


q(»-K) 


for  N-K-0 
N-K-1 
N-K  > 2 


(4.3.1.6) 


The  optlania  K for  N«6,  which  follows  from  Eqs.  (4. 3. 1.1)  and  (4. 3. 1.6) 
was  substituted  into  Eq.  (4. 3. 1.2)  and  the  result  Is  plotted  In  Fig.  4. 2. 1.2. 
As  can  be  seen  from  that  figure,  E{ll|q}  is  the  same  for  both  binary  and 
optimum  tree  for  small  q.  As  q Increases,  however,  the  optimum  protocol  Is 
definitely  superior.  What  this  suggests  Is  that  the  optimum  algorithm 
should  be  used  when  the  traffic  Is  heavy. 

When  the  tree  algorithm  Is  equivalent  to  the  TDMA  protocol.  Since 
q(8-K*0)  ■ l/yr  ■ .707  we  see  that  TDMA  Is  optimum  If  q > 1//2,  What  Is 
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surprising  Is  chat  Chls  result  does  not  depend  on  N.  The  relationship 
between  the  tree  algucicbm  and  TDMA  Is  considered  In  more  detail  In  Theorem 
4. 3. 2.1  In  the  following  section. 

We  conclude  this  section  %flth  the  presentation  of  the  optimal  dynamic 
strategy.  This  Is  as  follows: 

1.  Observe  1^^,  the  length  of  the  previous  epoch. 

2.  Substitute  1^  Into  Eq.  (4. 1.0.1)  to  obtain  q. 

3.  Use  the  q,  determined  In  Step-2,  to  obtain  R from  Eq.  (4. 3. 1.6). 

4.  Use  Che  tree,  determined  In  Step-3,  to  resolve  any  conflicts 

Chat  may  exist. 

4.3.2  Analysis  of  the  Optimum  Dynamic  Tree  Algorithm 

In  Chls  section  we  will  analyze  the  optimum  dynamic  algorithm  ^en  It 
Is  used  In  conjunction  with  the  finite  source  isodel.  First,  we  will  develop 
an  upper  bound  to  the  S{delay}  In  terms  of  p,  secondly,  we  will  develop  a 
corresponding  lower  bound  to  the  E{ throughput) , and  finally,  combine  these 
two  bounds  to  obtain  ad  vs  ch,  perforswnce  curve.  These  three  results  are 
Uluacrated  In  Figs.  4. 3.2. 3,  4. 3. 2. 4 and  4. 3. 2.5.  In  this  section,  we  will 
also  prove  (In  Theorem  4. 3. 2.1)  that  Che  E{delay}  of  the  optimum  dynamic 
tree  algorithm  la  smaller  than  or  equal  to  that  of  the  TOMA  protocol.  This 
Interesting  result  should  be  evident  from  the  work  of  the  preceding  section. 

1.  Upper  Bound  to  B( delay) 

As  In  Section  4.2.1  the  delay  la  decomposed  Into  d^  and  d2.  There- 
fore, we  can  write 

E{6|l,p,H,IC)  - E{dj^|l,p,H,K)  + E{d2|l,p,N,K)  (4.3.2. 1) 
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E{d^|  j.,p,N,K}  In  the  preceding  equation  Is  the  same  as  that  of  the  static 
algorltha.  That  Is, 


E{d jA.p.N.K}  - 1-1- 

Pd-a-l)*’) 


(4. 3.2.2) 


E{d2|l,p,M,K}  can  be  shown  to  be 


E{d2ll,p,S,K}  2 E{l2lqdtP).M-E+l,l) 

, , , (4.3.2. 3) 

E{d2|q(l,p),  N-K+l.l} 

In  the  preceding  equation,  q.  B{d2fq.«,l},  and  Ea2lq.a.l>  «re  given  by 
Eqs.  (4.1.0.1),  (4. 2.1.8),  and  (4.2.1.11). 

The  conditional  delay  given  by  Eq.  (4. 3.2.1)  was  coaputed  for  N-6 
and  K"K  (from  Table  4. 3. 1.1).  The  results  are  shown  in  Fig.  4. 3. 2.1.  As 
can  be  seen  from  that  figure,  B{« |i,p,M,E*}  Is  Increasing  In  t but  It  Is 
not  concave.  The  nonconcavity  of  this  function  Is  especially  evident  for 
P-.4  around  1-10.  Hext,  we  proceed  by  obtaining  the  tightest  con- 

cave lower  bound  to  E{6|l,p,H,K*}  and  than  upper  bounding  E{6>  by 


E{6)  < f^i(I„). 


(4.3.2. 4) 


The  function  'vl  Is  determined  graphically  from  Pig.  4. 3. 2.1  and  from  a 
more  detailed  computer  printout. 
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r 


Fl|ur«  4. 3. 2.1  E{d«l«y|l,p}  v 

Algorlctaa  with 


The  next  step  Is  to  calculate  It  can  be  seen  that 

^ * 

E{A2Ui»P»N,K  } - 2*^  E{42Ui.P,N-K*+1.1}  (4.3.2.5) 

where  the  conditional  mean  on  the  right  of  Eq.  (4. 3.2. 5)  Is  given  by 
Eq.  (4. 2. 1.8).  Equation  (4. 3. 2. 5)  Is  plotted  In  Fig.  (4. 3. 2. 2)  for  11-6. 

By  a procedure  similar  to  that  of  Section  4.2.1,  r^(p)  Is  calculated.  The 
results  are  listed  In  Table  4. 3. 2.1.  Finally  Is  substituted  Into 
Eq.  (4. 3. 2. 4)  and  the  upper  bound  to  E{6},  thus  derived.  Is  tabulated  In 
Table  4. 3. 2.1  and  plotted  In  Fig.  4. 3. 2. 3. 

11.  Lower  Bound  to  Ef throughput) 

Here  as  In  Section  4.2.2,  we  will  determine  th^  and  ^2*  Ihese 
two  quantities  are  defined  by  Eqs.  (4. 2. 2. 3)  and  (4. 2. 2. 9).  They  have 
been  calculated  using  the  values  of  given  In  Table  4. 3. 2.1,  and  the  re- 
sults are  presented  on  that  seme  table  and  In  Fig.  4. 3. 2. 4.  Finally,  by 
using  the  lower  bound  to  B{ throughput}  the  max{th2,  th2}  we  obtain  from 
Table  4. 3. 2.1  the  7^  vs  th^^  performance  curve  shown  In  Fig.  4. 3. 2. 5. 

Ill*  On  the  Superiority  of  the  Optimum  Dynamic  Tree  Over  the  TDK4  Protocol 
The  following  theorem  Is  based  on  the  observation  that  E{5|q} 

^ 0.  That  this  Is  so  follows  by  showing  that 

If  E{6|1}  ||>  0. 

Theorem  4. 3.2. It  Let  q be  the  probability  that  a source  has  a packet  to 
transmit  and  1st  6 be  the  delay  that  a packet  undergoes  when  It  Is  processed 
by  the  optimum  tree  algorithm.  Furthermore,  assume  that. 
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Figure  4. 3. 2. 2 van 
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.314 

.016 

6.8 
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.372 
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14.9 

.410 

.381 

,03 

19 

30.5 

.454 

.443 

.04 

32 

35.4 

.462 

.481 

05 

32 

36.2 

.505 

1 

32 

39.6 

.549 

2 

32 

43.7 

.600 

.4 

32 

46.2 

.700 

6 

32 

47.0 

.812 

TABLE  4. 3. 2.1 

AVERAGE  DELAY  AND  AVERAGE  THROUGHPUT  FOR  OPTIMUM 
DYNAMIC  TREE  ALGORZTH<  WITH  64  SOURCES 
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Figure  4. 3. 2. 4 Lower  Bounde  to  B{ throughput}  versua  p for  the 

OptlauB  Dynaalc  Tree  Algorltha  with  64  Sources 
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E[Thr.]  (packets/slot) 


Vlgur*  4.3. 


.3  Th«  Upper  Bound  to  E{d«lny}  versus  tlie  Lower  Bound 
to  E( throughput)  for  the  Optlaun  Dynaalc  Tree 
Algorlthn  with  64  Sources 

Note:  An  Algr.  step  equels  one  round  trip  delep 
plus  two  slots 


It  °* 


Than  th«  E{ delay}  of  the  optlania  dynamic  tree  protocol  la  leas  than  or  equal 
to  the  E{delay}  of  the  TDHi^  protocol. 

Proof: 

1 

E{6}  - / E{6|q}  dP(q) 

0 

Integrating  by  parts  we  have, 

1/  1 a 

E{6}  - E{6|q}  F(q)l  - /(|-E{6lq}  F(q)dq. 

0 0 

But  E{d|q}  0,  F(q)  ^ 0,  F(0)  > 0,  and  F(l)  > 1;  therefore 
B{6}  < E{«|q-l} 

t 

But  from  Eq.  (4. 3. 1.6)  it  follows  that  at  q"l  the  optimum  tree  algorithm 
is  the  TMA  protocol. 

QED 
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APPENDIX  A4 


PROPERTIES  OF  THE  BINARY  TREE/FINITE  SOURCE  SYSTEM 
Let  there  be  2^  sources  In  a multiaccess  system  where  each  source  may 
be  active  with  probability  q,  and  assume  that  the  multi-access  protocol  Is 
the  static  binary  tree  algorithm.  Also  let  I and  d have  the  same  definitions 
as  they  did  In  Chapter  2.  Then  In  this  appendix  we  will  derive  expressions 
for  E(l|q}  and  E{d|q}.  Since  the  work  of  this  appendix  parallels  that  of 
Appendix  A2,  we  will  not  be  as  detailed  here. 

A4.1  Derivation  of  E{1|q} 

The  number  of  nodes  I visited  by  the  algorithm  may  be  written  as  follows: 
N-1  2^-1 

4-1  + 1 I *44  (A4.1.1) 

1-1  j-0 

where 


*lj 


1 If  node  n^^  Is  visited 
0 otherwise 


CA4.1.2) 


Since  a node  n^^ 
we  have  that 


Is  visited  If  there  are  at 


least  two  active  sources  In  T 


Ij 


P(Xjj-l/q)  - ^(q»N-i) 


(A4.1.3) 


where 


'l»<q 


,m)  5 1 - (l-q)*"  - 2*q(l-q)^"l 


(A4.1.4) 
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and  It  follo«rs  that 


E{4|q} 


1 + 


N-l 

I 

i-1 


2"-  iKq.N-l) 


(A4.1.5) 


or 


B{A|q)  - 1 + 2D(q,N-l) 


(A4.1.6) 


where 

a-l  . 

D<q,a)  • J 2^  <{>(q,a-l)  (A4.1.7) 

1-0 

Equations  (A4.1.5)  through  (A4.1.7)  are  the  desired  results. 

A4.2  Derivation  of  E{d|q} 

E(d|q}  Is  the  number  of  nodes  that  ere  visited  before  a randonly  se- 
lacted  packet  from  the  contending  set  Is  successfully  transadtted.  First 
we  will  calculate  E{d  |q},  where  d Is  the  niad>er  of  nodes  visited  before 
source-s  Is  successfully  transmitted.  E{d^|q}  may  be  decomposed  as 

B{d^|q}  - (A4.2.1) 

Where  and  have  the  same  definitions  as  la  Appendix  A2.6.  That  Is*  X^ 
la  the  average  number  of  nodes  lying  on  s,  and  Is  the  average  number  of 
nodes  above  s that  were  visited  before  the  successful  transmission  of  s. 

As  In  Appendix  A2.6  we  may  write  X.  and  T as 

c.) 
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N-1 

X - 1 + I 'Kq,N-i)  (A4.2.2) 

® i-1 

and 

H-l 

^8  " h-l  ♦(q.N-1) 

* 1-1  ^ ^ 

vhara 

<Kq,m)  - 1 - (1-q)^^ 

■ Pr  {at  least  2 active  sources  In  a branch  of  depth 
m,  given  that  a particular  source  in  that  branch 
is  active). 

To  obtain  the  final  result,  first  substitute  Eqs.  (A4.2.2)  and  (A4.2.3)  into 
Bq.  (AA.2.1)  and  then  add  y E{d^|q}  and  E{d^|q},  as  shown  in  Eq.  (A4.2.3). 
Note  that  a is  the  ones*  conplaamt  of  s. 

Ill  1 

j ECd^lq)  + -j  {dj|q>  - 1 + [1  + y D(q,H-i)l  ♦(q.H-i)  (A4.2.3) 

- X{dlq)  (A4.2.6) 

The  last  step  follows  because  E{d^|q}  + ECd^lq)  is  independent  of  s. 


(A4.2.3) 


BIBLIOGRAPHY 


1.  Abramson,  Norman,  Packet  Switching  With  Satellites,  National  Computer 
Conference,  Vol.  42,  1973. 

2.  Roberts,  Lawrence  G.,  Aloha  Packet  System  with  and  without  Slots  and 
Capture,  ASS,  Mote  8,  June  26,  1972. 

3.  Metcalf,  Robert  M.,  Packet  Communication,  M.I.T.,  MAC  TR-114,  December 
1973. 

4.  Klelnrock,  Leonard  and  Lam,  Simon  S.,  Packet  Switching  In  a Multiaccess 
Broadcast  Channel:  Performance  Evaluation,  IEEE  Trans,  on  Communications, 
Vol.  COM-23,  No.  4,  April  1975. 

5.  Lam,  Simon  S.,  and  Klelnrock,  Leonard,  Packet  Switching  In  a Multiaccess 
Broadcast  Channel:  Dyanmlc  Control  Procedures,  IEEE  Trans,  on 
Communications,  Vol.  COM-23,  No.  9,  September  1975. 

6.  Carlelal,  Aydano  B.  and  Heilman,  Martin  B. , Bistable  Behavior  of  Aloha- 
Type  Systems.  Same  as  [4]. 

7.  Crowther,  W. , el  al.  A system  for  Broadcast  Communication:  Reservation- 
Aloha,  Proceedings  of  the  Sixth  Hawaii  International  Conference  on  System 
Sciences,  January  1973. 

8.  Roberts,  Lawrence  G.,  Dynamic  Allocation  of  Satellite  Capacity  Through 
Packet  Reservation.  National  Computer  Conference,  Vol.  42,  1973. 

9.  Binder,  Richard,  A Dynamic  Packet-Switching  System  for  Satellite  Broadcast 
Channels.  University  of  Hawaii,  Tech.  Report  B74-5,  August  1974. 

10.  Klelnrock,  Leonard  and  Tobagl,  Fouad  A. , Packet  Switching  In  Radio  Channels 
Part  I.  IEEE  Trans,  on  Communications,  Vol.  COM-23,  Mo.  12, 

Decamber  1975. 

11.  Tobagl,  Fouad  A.  and  Klelnrock,  Leonard,  Packet  Switching  In  Radio  Channels 
Part  11.  Same  as  [10]. 

12.  Drake,  Alvin  W. , Fundamentals  of  Applied  Probability  Theory,  New  York: 
McGraw-Hill,  1967. 


-223 


BIOGRAPHICAL  NOTE 


John  Ippocratls  Capetanakls  was  born  in  Greece  on  December  S,  1944. 

He  attended  public  schools  in  New  Britain,  Connecticut  and  graduated  from 
New  Britain  High  School  in  June,  1964. 

Mr.  Capetanakls  received  his  B.S.  and  M.S.  degrees  in  electrical 
engineering  fromM.I.T.  in  June,  1970.  From  September,  1968,  to  June,  1972, 
he  was  supported  by  a teaching  asslstantship.  In  June,  1972,  he  joined  the 
M.I.T.  Lincoln  Laboratory  staff  where  he  was  employed  until  September,  1975. 
From  September,  1975,  to  the  present  time  he  has  been  a full-time  graduate 
student  at  M.I.T.  supported  by  the  Lincoln  Laboratory  staff  associate 
program. 

Mr.  Capetanakls  is  a member  of  Tau  Beta  Pi,  Eta  Kappa  Nu,  and  Sigma  Xi. 


9 


